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Abstract 



This Thesis is based on quantum electrodynamic (QED) bound-state calculations [1, 2, 3, 4, 5, 6, 7, 8, 9, 
10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22], as well as on investigations related to divergent series, 
convergence acceleration and applications of these concepts to physical problems [23, 24, 25, 26, 27, 28, 
29, 30, 31]. 

The subjects which are discussed in this Thesis include: the self energy of a bound electron and the 
spin-dependence of QED corrections in bound systems (Chs. 1-5), convergence acceleration techniques 
(Chs. 6 and 7), and resummation methods for divergent series with an emphasis on physical applications 
(Ch. 8). 

In Chapter 1, we present numerical results [] 2], accurate to the level of 1 Hz for atomic hydrogen, for the 
energy correction to the K and L shell states due to the one-loop self energy. We investigate hydrogcnlike 
systems with low nuclear charge number Z = 1-5. Calculations arc carried out using on-mass shell 
rcnormalization, which guarantees that the final result is written in terms of the physical electron charge. 
The purpose of the calculation is twofold: first, we provide accurate theoretical predictions for the one-loop 
self energy of K and L shell states. Second, the comparison of the analytic and numerical approaches to 
the Lamb shift calculations, which have followed separate paths for the past few decades, is demonstrated 
by comparing the numerical values with analytic data obtained using the Za-expansion [1, 2, 32]. Our 
calculation essentially eliminates the uncertainty due to the truncation of the Za-expansion, and it 
demonstrates the consistency of the numerical and analytic approaches which have attracted attention 
for more than five decades, beginning with the seminal paper [-V^]. The most important numerical results 
are summarized in Table 1.5. 

In Chapter 2, we investigate higher-order analytic calculations of the one-photon self energy for excited 
atomic states. These calculations rely on mathematical methods described in Sec. 2.2 which, in physical 
terms, lead to a separation of the calculation into a high- and a low-energy part for the virtual photon. 
This separation does not only permit adequate simplifications for the two energy regions [1, 2], but also 
an adequate treatment of the infrared divergences which plague all bound-state calculations (see also 
Ch. 7 of [34]). The investigation represents a continuation of previous work on the subject [1, 2, 32]. The 
calculations are relevant for transitions to highly excited states, which are relevant for the extraction 
of fundamental constants from the high- precision measurements in atomic hydrogen [35], and for the 
estimation of QED effects in more complex atomic systems where some of the electrons occupy highly 
excited states. 

Chapter 3. We investigate the two-loop self energy. The calculation is based on a generalization of the 
methods introduced in Sec. 2.2. Historically, the two-loop self energy for bound states has represented a 
major task for theoretical evaluations. We present an analytic calculation [16] of the fine-structure differ- 
ence of the two- loop self energy for P states in atomic hydrogen in the order (Zoif' ra^ (? . This energy 
difference can be parameterized by two analytic coefficients, known as i?6i and i3go (see Sec. 3.2). These 
coefficients describe the two- loop self energy in the sixth order in Za, with an additional enhancement 
due to a large logarithm ln[(ZQ!)~^] (in the case of the Bgi-coefficicnt). The calculations arc relevant for 
an improved theoretical understanding of the fine-structure in hydrogenlike systems. They are also in part 
relevant for current experiments on atomic helium [36, 37, 38, 39], whose motivation is the determination 
of the fine-structure constant with improved accuracy. Finally, it is hoped that numerical calculations of 
the two-loop effect will be carried out in the near future for which the current analytic evaluation will be 
an important consistency check. The calculation of the analytic corrections discussed in Sec. 3 represents 
a solution for the most problematic set of diagrams on the way to advance our understanding of the 
fine-structure in atomic hydrogen to the few-Hz level. Explicit results for the fine-structure difference of 
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the Bqi- and i?60-coefScients can be found in Eqs. (3.59) - (3.64). 

In Chapter 4, We investigate the spin-dependence of the Breit hamiltonian and quantum electrodynamic 
effects in general. Specifically, we consider a bound system of two spinless particles. The calculation is 
motivated in part by current experiments (DIRAC at CERN) whose aim is the experimental study of 
pionium, which is the bound system of two (spinless) oppositely charged pious. The evaluation of the 
two-body relativistic corrections of order {Za)'^ leads to a different result [17] than expected for a system 
of two spin-1/2 particles of equal mass, e.g. positronium. In particular, we conclude in Sec. 4.3 that 
the so-called zitterbcwegung term is absent for a system of two spinless particles (the absence of the 
zitterbewegung term in a bound system consisting of a spinless and a spin-l/2 particle was previously 
pointed out in [40]). Final results for the relativistic correction, the vacuum polarization, and the self 
energy of a system of two scalar particles can be found in Eqs. (4.15), (4.21) and (4.31). A summary of 
the QED calculations is provided in Chapter 5. 

In the second part of this Thesis, we discuss methods for accelerating the convergence of series, and for 
the resummation of divergent series. As discussed in Chapter 6 and Sec. 7.1.1, convergence acceleration 
is essentially based on the idea that information hidden in trailing digits of elements of the sequence 
can be used in order to make "educated guesses" regarding the remainder term, which can be used 
for the construction of powerful sequence transformations (see Sec. 7.1.5). In Chapter 7, we discuss 
convergence acceleration in detail. After a short overview of relevant mathematical methods (Sec. 7.1), 
we discuss applications in applied biophysics (Sec. 7.2.1), in experimental mathematics (Sec. 7.2.2), and 
other applications, mainly in the evaluation of special functions (Sec. 7.2.3). In particular. Sees. 7.2.1 - 
7.2.3 illustrate how the combined nonlinear-condensation transformation (CNCT) described in Sec. 7.1.6 
can be used for the accelerated evaluation of nonalternating series, with an emphasis to applications of 
practical significance (Sec. 7.2.1). 

The discussion on convergence acceleration in Ch. 7 is complemented in Chapter 8 by an overview of re- 
summation techniques and relevant applications. Using the Stark effect and the associated autoionization 
width as a paradigmatic example, we discuss basic concepts like the Borel resummation method and its 
generalizations (Sec. 8.2.2), and the conformal mapping of the complex plane (Sec. 8.4). We then proceed 
to discuss further applications of resummation methods like zero-dimensional model theories (Sec. 8.3.1), 
the QED effective action (expressed as a perturbation series in the fine structure constant, see Sec. 8.3.2), 
and the quantum-mechanical double-well problem (Sec. 8.3.3). Within the context of the double- well 
problem, we perform an analytic evaluation of higher-order corrections to the two-instanton effect and 
demonstrate the consistency of numerically determined energy levels with the instanton expansion. 

We conclude with a summary of the results in Chapter 9, where we also explain the interrelations and 
connections between the different subjects treated in this Thesis. 
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Part I 



Quantum Electrodynamic 
Bound— State Calculations 
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Chapter 1 



Numerical Calculation of the 
One— Loop Self Energy (Excited 
States) 

1.1 Orientation 

A nonperturbative numerical evaluation of the one-photon electron self energy for the K- and L-shell 
states of hydrogenlike ions with nuclear charge numbers Z = 1 to 5 is described. Our calculation for the 
lSi/2 state has a numerical uncertainty of 0.8 Hz in atomic hydrogen, and for the L-shell states (2Si/2! 
2Pi/2i E^nd 2P3/2) the numerical uncertainty is 1.0 Hz. The method of evaluation for the ground state 
and for the excited states is described in detail. The numerical results are compared to results based on 
known terms in the expansion of the self energy in powers of Za. 

1.2 Introduction to the Numerical Calculation of Radiative Cor- 
rections 

The nonperturbative numerical evaluation of radiative corrections to bound-state energy levels is inter- 
esting for two reasons. First, the recent dramatic increase in the accuracy of experiments that measure 
the transition frequencies in hydrogen and deuterium [35, 41, 42] necessitates a numerical evaluation 
(nonperturbative in the binding Coulomb field) of the radiative corrections to the spectrum of atomic 
systems with low nuclear charge Z . Second, the numerical calculation serves as an independent test of 
analytic evaluations which are based on an expansion in the binding field with an expansion parameter 
Za. 

In order to address both issues, a high-precision numerical evaluation of the self energy of an electron 
in the ground state in hydrogenlike ions has been performed [9, 12, 43]. The approach outlined in [9] is 
generalized here to the L shell, and numerical results are obtained for the (n = 2) states 2Si/2, 2Pi/2 and 
2P3/2- Results are provided for atomic hydrogen. He"*", Li^"*", Be^^, and B''+. 

It has been pointed out in [9, 43] that the nonperturbative effects (in Za) can be large even for low 
nuclear charge and exceed the current experimental accuracy for atomic transitions. For example, the 
difference between the sum of the analytically evaluated terms up to the order of a (Za)^ and the final 
numerical result for the ground state is roughly 27 kHz for atomic hydrogen and about 3200 kHz for He'*'. 
For the 2S state the difference is 3.5 kHz for atomic hydrogen and 412 kHz for He"*". The large difference 
between the result obtained by an expansion in Za persists even after the inclusion of a result recently 
obtained in [44] for the logarithmic term of order a {ZaY ln(ZQ;)~^. For the ground state, the difference 
between the all-order numerical result and the sum of the perturbative terms is still 13 kHz for atomic 
hydrogen and 1600 kHz for He"''. For the 2 S state, the difference amounts to 1.6 kHz for atomic hydrogen 
and to 213 kHz for He+. 
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These figures should be compared to the current experimental precision. The most accurately measured 
transition to date is the 1S-2S frequency in hydrogen; it has been measured with to a relative accuracy 
of 1.8 parts in lO^'' or 46 Hz [42]. This experimental progress is due in part to the use of frequency 
chains that bridge the range between optical frequencies and the microwave cesium time standard. The 
accuracy of the measurement is likely to be improved by an order of magnitude in the near future [42, 45] . 
With trapped hydrogen atoms, it should be feasible to observe the 1S-2S frequency with an experimental 
linewidth that approaches the 1.3 Hz natural width of the 2S level [46, 47]. 

The apparent convergence of the perturbation series in Za is slow. Our all-order numerical calculation 
presented here essentially eliminates the uncertainty from unevaluated higher-order analytic terms, and 
we obtain results for the self energy remainder function Gse with a precision of roughly 0.8 x Hz for 
the ground state of atomic hydrogen and 1.0 x Z'^Hz for the 2S state. 

In the evaluation, we take advantage of resummation and convergence acceleration techniques. The re- 
summation techniques provide an efficient method of evaluation of the Dirac-Coulomb Green function to 
a relative accuracy of 10~^* over a wide parameter range [4-3]. The convergence acceleration techniques 
remove the principal numerical difficulties associated with the singularity of the relativistic propagators 
for nearly equal radial arguments [23]. 

The one-photon self energy is about two orders of magnitude larger than the other contributions to 
the Lamb shift in atomic hydrogen. Comprehensive reviews of the various contributions to the Lamb 
shift in hydrogenlike atoms in the full range of nuclear charge numbers Z = 1-110 have been given in 
Refs. [48, 49, 50, 51]. 

1.3 Method of Evaluation 
1.3.1 Status of Analytic Calculations 

The (real part of the) energy shift Ai?sE due to the electron self energy radiative correction is usually 
written as ^ 

AEsE^-^-^F{nl,,Za)m,c^ (1.1) 

where is a dimensionlcss quantity. In the following, the natural unit system with h = c = nic = 1 and 
= 47ra is employed. Note that F{nlj, Za) is a dimensionlcss function which depends for a given atomic 
state with quantum numbers n, I and j on only one argument (the coupling Za). For excited states, 
the (nonvanishing) imaginary part of the self energy is proportional to the (spontaneous) decay width of 
the state. We will denote here the real part of the self energy by Ai^sE, exclusively. The semi-analytic 
expansion of Finlj, Za) about Za = for a general atomic state with quantum numbers n, I and j gives 
rise to the expression, 

F{nlj,Za) = A4i{nlj) ln{Zay^ 

+A4oinlj) + (Za) Ao(n^j) 
+ {Zaf [Ae2{nlj) \n{Za)-^ 

+Aeiinlj) \n{Za)-^ + GsEinlj, Za)] . (1.2) 

For particular states, some of the coefficients may vanish. Notably, this is the case for P states, which are 
less singular than S states at the origin [see Eq. (1.4) below]. For the nSi/2 state {I = 0, j = 1/2), none 
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of the terms in Eq. (1.2) vanishes, and we have, 

F{7jSi/2,Za) = ^41(7181/2) ln(Za)~^ 



+A4o(nSi/2) + (Za) Ar^oinS 1/2) 



+ {ZaY [A62(nSi/2) \n\Za)- 



+A6i(r7,Si/2) ln(ZQ;) + Gs^{nSi/2, Za)] 



(1.3) 



The A coefficients have two indices, the first of which denotes the power of Za [including those powers 
imphcitly contained in Eq. (1.1)], while the second index denotes the power of the logarithm ln(Za)~^. 
For P states, the coefficients A41. A50 and Aq2 vanish, and we have 



F{7iPj, Za) ^ AiQ{nPj) + {Zaf [^6i(nPj) ln(Za)-2 + GsE(nPj, Za)] . 
For S states, the self energy remainder function Gse can be expanded semi-analytically as 
Gse(»^Si/2, Za) = A6o(riSi/2) + {Za) [^71(71,81/2) ln(ZQ;)"^ + ^70(77.81/2) + o{Za)] 



(1.4) 



(1.5) 



For the "order" symbols o and O we follow the usual convention [")2, -53]: the requirement is 0{x)/x 
const, as a; ^ 0, whereas o{x) fulfills the weaker requirement o(a;) ^ as a; ^ 0. For example, the 
expression [{Za) ln(ZQ;)] is o{Za) but not 0{Za). Because logarithmic terms corresponding to a nonva- 
nishing ^83-coefficient must be expected in Eq. (1.5), the symbol o{Za) should be used to characterize 
the remainder, not 0{Za). 

For P states, the semi-analytic expansion of Gse reads 



GsE{nFj,Za) = Aeo{nPj){Za) [^70(77?^) + o{Za)] 



(1.6) 



The fact that ^71(77?^) vanishes has been pointed out in [44]. We list below the analytic coefficients 
and the Bethe logarithms relevant to the atomic states under investigation. For the ground state, the 
coefficients A41 and A4Q were obtained in [33, 54, 55, 56, 57, 58, 59], the correction term ^50 was found 
in [60, 61, 62], and the higher-order binding corrections Ae2 and Aqi were evaluated in [32, 63, 64, 65, 
66, 67, 68, 69, 70, 71]. The results are. 



^41(181/2) 



4 
3 ' 



^40(181/2 



^-^lnfco(lS), 



A50(18i/2 



in 



139 1 , „ 

In 2 

128 2 



^62(lSi/2) = -1 



28 21 

A6i(iSi/2) - yi^2--. 

The Bethe logarithm Infco(lS) has been evaluated in [72] and [73, 74, 75, 76, 77] as 

lnA:o(lS) = 2.984 128 555 8(3). 



(1.7) 
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For the 2S state, we have 



v44i(2Si/2) 



3' 



^40(281/2) 



^-^lnfco(2S), 



^50(2Si/2) - 4^ 



139 1 , „ 

n 2 

128 2 



^2(281/2) - -1 



^6l(2Si/2) 



16 , „ 67 

— In 2 H . 

3 30 



(1.9) 



The Bethc logarithm hrfco(2S) has been evaluated (see [72, 73, 74, 75, 76, 77], the results exhibit varying 
accuracy) as 

In /so (2S) = 2.811 769 893(3). (1.10) 

It might be worth noting that the value for In /co(2S) given in [78] evidently contains a typographical error. 
Our independent re-evaluation confirms the result given in Eq. (1.10), which was originally obtained in [72] 
to the required precision. For the 2Pi/2 state we have 



A40(2Pl/2) 



-i-^ lnfco(2P). 



A6i(2P 



1/2J 



103 
108 



(1.11) 



Note that a general analytic result for the logarithmic correction Aqi as a function of the bound state 
quantum numbers n, I and j can be inferred from Eq. (4.4a) of [68, 69] upon subtraction of the vacuum 
polarization contribution implicitly contained in the quoted equation. The Bethe logarithm for the 2P 
states reads [72, 79] 

lnfco(2P) = -0.030 016 708 9(3) . (1.12) 

Because the Bethe logarithm is an inherently nonrelativistic quantity, it is spin-independent and therefore 
independent of the total angular momentum j for a given orbital angular momentum /. For the 2P3/2 
state the analytic coefRcients are 



^40(2P3/2) 



^-^lnM2P). 



A6l(2P 



3/2) 



29 

90 



(1.13) 



We now consider the limit of the function GsE{Za) as Za 0. The higher-order terms in the potential 
expansion (see Fig. 1.3 below) and relativistic corrections to the wavefunction both generate terms of 
higher order in Za which are manifest in Eq. (1.2) in the form of the nonvanishing function GsEiZa) 
which summarizes the effects of the relativistic corrections to the bound electron wave function and of 
higher-order terms in the potential expansion. For very soft virtual photons, the potential expansion fails 
and generates an infrared divergence which is cut off by the atomic momentum scale, Za. This cut-off 
for the infrared divergence is one of the mechanisms which lead to the logarithmic terms in Eq. (1.2). 
Some of the nonlogarithmic terms in relative order (Za)^ in Eq. (1.2) are generated by the relativistic 
corrections to the wave function. The function Gse does not vanish, but approaches a constant in the 
limit Za 0. This constant can be determined by analytic or semi-analytic calculations; it is referred 
to as the Ago coefficient, i.e. 



(1.14) 
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The evaluation of the coefficient A6o(lSi/2) has drawn a lot of attention for a long time [32, 68, 69, 70, 71]. 
For the 2S state, there is currently only one accurate analytic result available, 

Ao(2Si/2) = -31.84047(1) (sec Ref. [32]). (1.15) 

For the 2Pi/2 state, the analytically obtained result is 

^6o(2Pi/2) = -0.998 91(1) (see Ref. [1]), (1.16) 
and for the 2P3/2 state, we have 

^6o(2P3/2) = -0.503 37(1) (see Ref. [1]), (1.17) 

The analytic evaluations essentially rely on an expansion of the relativistic Dirac-Coulomb propagator 
in powers of the binding field, i.e. in powers of Coulomb interactions of the electron with the nucleus. In 
numerical evaluations, the binding field is treated nonperturbatively, and no expansion is performed. 



1.3.2 Formulation of the Numerical Problem 



Numerical cancellations are severe for small nuclear charges. In order to understand the origin of the 
numerical cancellations it is necessary to consider the renormalization of the self energy. The renormal- 
ization procedure postulates that the self energy is essentially the effect on the bound electron due to the 
self interaction with its own radiation field, minus the same effect on a free electron which is absorbed 
in the mass of the electron and therefore not observable. The self energy of the bound electron is the 
residual effect obtained after the subtraction of two large quantities. Terms associated with renormaliza- 
tion counterterms are of order 1 in the Za-expansion, whereas the residual effect is of order {Za)^ [see 
Eq. (1.1)]. This corresponds to a loss of roughly 9 significant digits at Z ~ 1. Consequently, even the 
precise evaluation of the one-photon self energy in a Coulomb field presented in [80] extends only down 
to Z = 5. Among the self energy corrections in one-loop and higher-loop order, numerical cancellations 
in absolute terms are most severe for the one-loop problem because of the large size of the effect of the 
one-loop self energy correction on the spectrum. 

For our high-precision numerical evaluation, we start from the regularized and renormalized expression 
for the one-loop self energy of a bound electron, 



SE 



lim 



'Re 



duj f d?k 



27r J (27r)3 



X ( V' 



1 



t/j)— Am 



lim 

A— ♦00 



^Re / ^ 



did f d^k 



2n J (27r) 



1 



X (V- 1 e"=-^ G{En - to) a" c-'^"' \ ^j) - Am }> , 
where G denotes the Dirac-Coulomb propagator, 

G{z) 



a- p + (3 + V ~ z ' 

and Am is the A-dependent (cutoff-dependent) one-loop mass-counter term. 



(1.18) 



(1.19) 



(1.20) 
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Im((D) 



-CL 



.E„, 



.2m 



Figure 1.1: Integration contour C for the integration over the energy lu ~ En — z of the 
virtual photon. The contour C consists of the low-energy contour Cl and the high-energy 
contour Ch- Lines shown displaced directly below and above the real axis denote branch 
cuts from the photon and electron propagator. Crosses denote poles originating from the 
discrete spectrum of the electron propagator. The contour used in this work corresponds to 
the one used in [Ml. 



The photon propagator D^^(fc^, A) in Eq. (1.18) in Feynman gauge reads 



fc2 



A2 



(1.21) 



The contour Cp in Eq. (1.18) is the Feynman contour, whereas the contour C is depicted in Fig. 1.1. The 
contour C is employed for the cj-integration in the current evaluation [see the last line of Eq. (1.18)]. The 
energy variable z in Eq. (1.19) therefore assumes the value 



E„ 



(1.22) 



where En is the Dirac energy of the atomic state, and w denotes the complex- valued energy of the virtual 
photon. It is understood that the limit A ^ oo is taken after all integrals in Eq. (1.18) are evaluated. 

The integration contour for the complex-valued energy of the virtual photon lo in this calculation is 
the contour C employed in [80, 81, 82, 83] and depicted in Fig. 1.1. The integrations along the low- 
energy contour Cl and the high-energy contour Ch in Fig. 1.1 give rise to the low- and the high-energy 
contributions AEi^ and AE^ to the self energy, respectively. Here, we employ a further separation of the 
low-energy integration contour Cl into an infrared contour Cir and a middle-energy contour Cm shown 
in Fig. 1.2. This separation gives rise to a separation of the low-energy part ASl into the infrared part 
Aii'iR and the middle-energy part AEm, 



AEl = A-BiR + AEi 



M 



(1.23) 



For the \ow-Z systems discussed here, all complications which arise for excited states due to the decay 
into the ground state are relevant only for the infrared part. Except for the further separation into the 
infrared and the middle-energy part, the same basic formulation of the self energy problem as in [81] is 
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lm{uj) ' 


CiR 


Cm 




XKXX X 















' 1 — 

CiR 


Cm 





Figure 1.2: Separation of the low-energy contour Cl into the infrared part Cir and the middle- 
energy part Cm- As in Fig. 1.1, the lines directly above and below the real axis denote branch cuts 
from the photon and electron propagator. Strictly speaking, the figure is valid only for the ground 
state. For excited states, some of the crosses, which denote poles originating from the discrete 
spectrum of the electron propagator, are positioned to the right of the line Rccj = 0. These poles 
are subtracted in the numerical evaluation. 



used. This leads to the following separation, 

Lo e (0, Vio -En) ± i e : infrared part AE^ir, 

oj G (Vio Em £"„) ± i e : middle-energy part Ai?M, 

oj e En + i (— oo, -l-oo) : high-energy part AE'h. 

Integration along these contours gives rise to the infrared, the middle-energy, and the high-energy con- 
tributions to the energy shift. For all of these contributions, lower-order terms are subtracted in order to 
obtain the contribution to the self energy of order (Za)^. We obtain for the infrared part, 



IR 



a 

TT 



21 



200" ' 600 

where Fi^(nlj, Za) is a dimensionlcss function of order one. The middle-energy part is recovered as 



(1.24) 



A^M = - 



and the high-energy part reads [81, 82] 
AEu = Am + 



-l{f3)-l{V)+^-^F^{nl„Za) 



(1.25) 



(1.26) 



The infrared part is discussed in Sec. 1.4.1. The middle-energy part is divided into a middle-energy 
subtraction term -Fma and a middle-energy remainder Fmb- The subtraction term -Fma is discussed in 
Sec. 1.4.2, the remainder term Fmb is treated in Sec. 1.4.3. We recover the middle-energy term as the 
sum 

Fuinlj^Za) = FuAinlj, Za) + Funinlj, Za) . (1-27) 

A similar separation is employed for the high-energy part. The high-energy part is divided into a subtrac- 
tion term Fha, which is evaluated in Sec. 1.5.1, and the high-energy remainder Fhb, which is discussed 
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in Sec. 1.5.2. The sum of the subtraction term and the remainder is 



Fninlj, Za) = FhaHj, Za) + Fusinlj, Za) . (1.28) 
The total energy shift is given as 

A£;sE = A£;iR, + I^Em + -Eh - Am 

= --^^ [Fininlj,Za) + FM{nlj,Za) + F]i{nlj, Za)] . (1.29) 

The scaled self energy function F defined in Eq. (1.1) is therefore obtained as 

F{nlj, Za) = FiR(nZj, Za) + Fu{nlj, Za) + F}i{nlj, Za) . (1.30) 

In analogy to the approach described in [80, 81, 83], we define the low-energy part as the sum of the 
infrared part and the middle-energy part, 

ASl = ASir + ASm 



a 

TT 



4 



lm + l{V)+^-^F^{nl„Za) 
2 6 n-^ 



(1.31) 



where 

Fi^{nlj,Za) = Fin{nlj,Za) + Fuinlj^Za) . (1.32) 
The limits for the functions F\^{nlj, Za) and FB_{nlj, Za) as Za were obtained in [43, 82, 84]. 

1.3.3 Treatment of the divergent terms 

The free electron propagator, 

F= 1- , (1.33) 

a ■ p + p — z 

and the full electron propagator G defined in Eq. (1.19), fulfill the following identity which is of particular 
importance for the validity of the method used in the numerical evaluation of the all-order binding 
correction to the Lamb shift, 

G = F - FV F + FVGV F . (1.34) 

This identity leads naturally to a separation of the one-photon self energy into a zero-vertex, a single- 
vertex and a many- vertex term. This is represented diagrammatically in Fig. 1.3. 

All ultraviolet divergences which occur in the one-photon problem (mass counter term and vertex diver- 
gence) are generated by the zero- vertex and the single-vertex terms. The many-vertex term is ultraviolet 
safe. Of crucial importance is the observation that one may additionally simplify the problem by replac- 
ing the one-potential term with an approximate expression in which the potential is "commuted to the 
outside" . The approximate expression generates all divergences and all terms of lower order than a {Za)^ 
present in the one-vertex term. Unlike the raw one-potential term, it is amenable to significant further 
simplification and can be reduced to one-dimensional numerical integrals which can be evaluated easily 
(a straightforward formulation of the self energy problem requires a i/iree-dimensional numerical integra- 
tion). Without this significant improvement, an all-order calculation would be much more difficult at low 
nuclear charge, because the lower-order terms would introduce significant further numerical cancellations. 

Furthermore, the special approximate resolvent can be used effectively for an efficient subtraction scheme 
in the middle-energy part of the calculation. In the infrared part, such a subtraction is not used because 
it would introduce infrared divergences. 

We now turn to the construction of the special approximate resolvent, which will be referred to as Ga 
and will be used in this calculation to isolate the ultraviolet divergences in the high-energy part (and 
to provide subtraction terms in the middle-energy part). It is based on an approximation to the first 
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Figure 1.3: The exact expansion of the bound electron propagator in powers of the binding 
field leads to a zero-potential, a one-potential and a many-potential term. The dashed lines 
denote Coulomb photons, the crosses denote the interaction with the (external) binding 
field. 



two terms on the right-hand side of Eq. (1.34). The so-called one-potential term FVF in Eq. (1.34) is 
approximated by an expression in which the potential terms V are commuted to the outside, 

-FVF « {V,F^} . (1.35) 

Furthermore, the following identity is used, 

F^ = ^ ' 



a. ■ p + f3 — z 



2z{f3 + z) 



2 z {ex. ■ p) 



(p2 + 1 _ ^2 



(1.36) 



In 2 X 2 spinor space, this expression may be divided into a diagonal and a non-diagonal part. The 
diagonal part is 



.2^_ 1 , 2z(/3 + z) 



diag(i^") = — — + ^ ^ . (1.37) 

p^ + l-z^ (p2 + l-^2)2 



The off-diagonal part is given by 

F'^ - diag(F^) 



2 J = „„.^.2^_ 2 Z {a ■ p) 



(p2 + 1 _ ^2 

We define the resolvent Ga as 

Ga = F-^ {V,diag{F^)} . (1.38) 
All divergences which occur in the self energy are generated by the simplified propagator Ga- We define 



16 



the propagator Gb as the difference of G and Ga, 

Gb = G — Ga 



= ^ {V,di&g{F^)} ^ FV F + FVG VF . (1.39) 

Gb docs not generate any divergences and leads to the middle-energy remainder discussed in Sec. 1.4.3 
and the high-energy remainder (Sec. 1.5.2). 

1.4 The Low-Energy Part 

1.4.1 The Infrared Part 

The infrared part is given by 

X (t/. I e"=-^ G{En ~ t^) o!" e"""'^ | i>) , (1.40) 

where relevant definitions of the symbols can be fomid in Eqs. (1.18-1.21), the contour Gir, is as shown 
in Fig. 1.2. and the unregularizcd version of the photon propagator 

/?..(fc^)=-T#^ (1.41) 
A; + 1 e 

may be used. The infrared part comprises the following integration region for the virtual photon (contour 
Gir in Fig. 1.2), 

^ e (0,Vio^„) ±ie 



(7l0Sn,£^n) ±ie 



infrared part ASir . (1-42) 



Following Sees. 2 and 3 of [M], we write ASjr as a three-dimensional integral [see, e.g., Eqs. (3.4), (3.11) 
and (3.14) ihid\ 

A-BiR = --^--(P.V.) / dz / dxix? / dxsx^ Xir(.t2,xi,z), (1.43) 

10 J°/ioE„ Jo Jo 

where 

2 

= E E .Ma^2)G:f"(x2,xi,z)/,>i)Ajf"(x2,xi). (1.44) 

K i,j = l 

Here, the quantum number k is the Dirac angular quantum number of the intermediate state, 

K^2{l-j){j + l/2), (1.45) 

where I is the orbital angular momentum quantum number and j is the total angular momentum of 
the bound electron. The functions fi{x2) {i = 1,2) are the radial wave functions defined in Eq. (A. 4) 
in [81] for an arbitrary bound state (and in Eq. (A. 8) in [81] for the IS state). We define i = 3 — i. The 
functions G^^(x2,xi,z) (i,j = 1,2) are the radial Green functions, which result from a decomposition of 
the electron Green function defined in Eq. (1.19) into partial waves. The explicit formulas are given in 
Eq. (A. 16) in [81], and we do not discuss them in any further detail, here. 

The photon angular functions A]^ = 1,2) are defined in Eq. (3.15) of Ref. [81] for an arbitrary 
bound state. In Eq. (3.17) in [81], specific formulas are given for the IS state. In Eqs. (2.2), (2.3) and 
(2.4) of [n.')], the special cases of S1/2, P1/2 and P3/2 states are considered. Further relevant formulas for 
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excited states can be found in [n")]. The photon angular functions depend on the energy argument z, but 
this dependence is usually suppressed. The summation over k in Eq. (1.44) extends over all negative and 
all positive integers, excluding zero. We observe that the integral is symmetric under the interchange of 
the radial coordinates X2 and xi, so that 



dz j dxixf / dx2X2 Mm,{x2,xi,z) . 



Ai^iH = ^ - ^ (P.V.) 
The following variable substitution, 

r = X2/xi, y = axi, 
is made, so that r G (0, 1) and y G (0, oo). The scaling variable a is defined as 



The Jacobian is 



d{x2,xi) 



d{r, y) 



dx2 dxi 

dr dr 

dx2 dxi 

~3y ~dy 



The function Syr is given by. 



2 „,5 



c I \ r . . fry y \ 

^m{r,y,z) = g— XirI — 

a° \ a a y 



(1.46) 
(1.47) 
(1.48) 

(1.49) 



2^2 y5 



OO z 



|k|=1 K=±|re| i,j=l 



V a a / Va/ \ a a 



2^2 y5 



|k| = 1 



where in the last line we define implicitly the terms Tir for |k| = 1, . . . ,cxd as 

2 

7iR,|«|(r,?/,z) = 

k=±|k| i,j=l 



(1.50) 



V a / \ a a / \a/ \ a a) 
Using the definition (1.50), we obtain for Ai?iR, 



Ai;iR - 



TT 10 



(P.V. 



(1.51) 



dr 



dy 5'iR(r,?/, z) . 



(1.52) 



The specification of the principal value (P.V.) is necessary for the excited states of the L shell, because 
of the poles along the integration contour which correspond to the spontaneous decay into the ground 
state. Here we are exclusively concerned with the real part of the energy shift, as specified in Eq. (1.40), 
which is equivalent to the specification of the principal value in (1.52). Evaluation of the integral over 
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z is facilitated by the subtraction of those terms which generate the singularities along the integration 
contour (for higher excited states, there can be numerous bound state poles, as pointed out in [2, 8 j]). 
For the 2S and 2Pi/2 states, only the pole contribution from the ground state must be subtracted. For the 
2P3/2 state, pole contributions originating from the IS, the 2S and the 2Pi/2 states must be taken into 
account. The numerical evaluation of the subtracted integrand proceeds along ideas outlined in [83, 85] 
and is not discussed here in any further detail. 

The scaling parameter a for the integration over y is chosen to simplify the exponential dependence of 
the fmiction S defined in Eq. (1.50). The main exponential dependence is given by the relativistic radial 
wave functions (upper and lower components). Both components \]\(x) and /2(a^)] vary approximately 
as (neglecting relatively slowly varying factors) 

exp (— a2:/2) (for large x) . 

The scaling variable a, expanded in powers of Za, is 



2^1^ 



2n2 



O {{ZoiY 



= 2 — + 0\{Zaf] . (1.53) 
Therefore, a is just twice the inverse of the Bohr radius n/{Za) in the nonrelativistic limit. The product 

A (^) ^ h (9 for arbitrary J € {1, 2} 

[which occurs in Eq. (1.50)] depends on the radial arguments approximately as 

X exp [Y2 (1 — r) y] (for large y) . 

Note that the main dependence as given by the term exp(— y) is exactly the weight factor of the Gauf5- 
Laguerre integration quadrature formula. The deviation from the exact cxp(— y)-type behavior becomes 
smaller as r ^ 1. This is favorable because the region near r = 1 gives a large contribution to the integral 
in (1.52). 



Z i^iR(lSi/2,Za) FiR(2Si/2,Za) Fir(2Pi/2, Za) Fir(2P3/2, Za) 

1 7.236 623 736 8(1) 7.479 764 180(1) 0.085 327 852(1) 0.082 736 497(1) 

2 5.539 002 119 1(1) 5.782 025 637(1) 0.086 073 669(1) 0.083 279 461(1) 

3 4.598 155 821 8(1) 4.840 923 962(1) 0.087 162 510(1) 0.084 091 830(1) 

4 3.963 124 140 6(1) 4.205 501 798(1) 0.088 543 188(1) 0.085 140 788(1) 

5 3.493 253 319 4(1) 3.735 114 958(1) 0.090 180 835(1) 0.086 403 178(1) 

Table 1.1: Infrared part for the K and L shell states, Fir(1Si/2, Za), Fir(2Si/2, Za), Fir(2Pi/2, Za), 
and i^iR(2P3/2, -^a), evaluated for low-Z hydrogcn(-like) ions. The calculations are performed with the 
numerical value of = 137.036 for the fine structure constant. 



The sum over \k\ in Eq. (1.50) is carried out locally, i.e. for each set of arguments r,y,z. The sum over 
\k\ is absolutely convergent. For \k\ — > 00, the convergence of the sum is governed by the asymptotic 
behavior of the Besscl functions which occur in the photon functions A]^ = 1, 2) [see Eqs. (3.15) and 
(3.16) in [81]]. The photon functions contain products of two Bessel functions of the form J7;(p2/i) where 
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Ji stands for either ji or the index I is in the range I € — 1,|k|,|k| + 1}. The argument is either 
P2 = [En — z) X2 or pi = [En — z) xi- The asymptotic behavior of the two relevant Bessel functions for 
large I (and therefore large is 



X {2i + iy. 



1 + 



and 



(1.54) 



(2^ + 1)!! 



1 + 



(1.55) 



This implies that when min{p2, Pi} = P2 < I, the function J^i{p2) vanishes with increasing I approximately 
as (ep2/2/)'. This rapidly converging asymptotic behavior sets in as soon as I k, \k\ > p2 ^ ruoy/a [see 
Eqs. (1-22) and (1.51)]. Due to the rapid convergence for |k| > p2, the maximum angular momentum 
quantum number \k\ in the numerical calculation of the infrared part is less than 3 000. Note that because 
z e {^/ 10 En, En) in the infrared part, lo < Vio En- 

The integration scheme is based on a crude estimate of the dependence of the integrand 5'iR(r, y, z) defined 
in Eq. (1.50) on the integration variables r, y and z. The main contribution to the integral is given by 
the region where the arguments of the Whittaker functions as they occur in the Green function [see 
Eq. (A. 16) in [81]] are much larger than the Dirac angular momentum, 

2c ^» \k\ 
a 

(see also p. 56 of [82]). We assume the asymptotic form of the Green function given in Eq. (A. 3) in [82] 
applies and attribute a factor 

exp[— (1 — r) cy/ a] 

to the radial Green functions GJ_^ as they occur in Eq. (1.50). Note that relatively slowly varying factors 
are replaced by unity. The products of the radial wave functions /j and /j, according to the discussion 
following Eq. (1.53), behave as 

e-y exp[V2(l-r)2/] 

for large y. The photon functions A^i* in Eq. (1.50) give rise to an approximate factor 

sin[(l-r) {En''z)y/a] 

■ ^ ' 

Therefore [sec also Eq. (2.12) in [82]], we base our choice of the integration routine on the approximation 

c _ r 

a 2. 



exp 



(1 -r)y 



^^ sin[(l-r) {{En-z)y/a}] 
(l-r) 

for 5iR . The three-dimensional integral in (1.52) is evaluated by successive Gaussian quadrature. Details 
of the integration procedure can be found in [4.3]. 

In order to check the numerical stability of the results, the calculations are repeated with three different 
values of the fine structure constant a, 

a< = 1/137.036 000 5, 

ao = 1/137.036 000 and (1.58) 
a> = 1/137.035 999 5. 

These values of the fine-structure constant are close to the 1998 CODATA recommended value of = 
137.035 999 76(50) [8(3]. The calculation was parallelized using the Message Passing Interface (MPI) 
and carried out on a cluster of Silicon Graphics workstations and on an IBM 9276 SP/2 multiprocessor 
system. The results for the infrared part. Fir defined in Eq. (1.24), are given in Table 1.1 for a value of 
a^^ = a^^ = 137.036. This value of a will be used exclusively in the numerical evaluations presented 
here. For numerical results obtained by employing the values of a< and Q!> [see Eq. (1.58)] we refer 
to [43]. 
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1.4.2 The Middle-Energy Subtraction Term 



The middle-energy part is given by 



X (V' I a" e"=-" G{En - lo) a" c-"' '' \ ^/.) , (1.59) 

where relevant definitions of the symbols can be found in Eqs. (1.18)-(1.21) and in Eq. (1.41). The middle- 
energy part comprises the following integration region for the virtual photon (contour Cm in Fig. 1.2), 

uje (Vio^„,^n) ±ie 1 

> middle-energy part AEm ■ (1.60) 

(0,7io£;„) ±ie J 

The numerical evaluation of the middle-energy part is simplified considerably by the decomposition of 
the relativistic Dirac-Coulomb Green function G as 



G = Ga + Gb, (1.61) 

where Ga is defined in (1.38) and represents the sum of an approximation to the so-called zero- and 
one-potential terms generated by the expansion of the Dirac-Coulomb Green function G in powers of 
the binding field V. We define the middle-energy subtraction term Fma sls the expression obtained upon 
substitution of the propagator Ga for G in Eq. (1.59). The propagator Gb is simply calculated as the 
difference of G and Ga [see Eq. (1.39)]. A substitution of the propagator Gb for G in Eq. (1.59) leads to 
the middle-energy remainder ^mb which is discussed in Sec. 1.4.3. We provide here the explicit expressions 

. ^ 9 /" /■ d^fc 9, 



X (V- 1 e''^-'" Ga(K - w) a'' e^''"'^ | V) 



(1.62) 



and 



c„2^y (2^ 



MB = -le^ / — l7—^D,,,{k'') 



X (tA I a'' e"=-=" Gb(£;„ - uj) a"" e'^^'^ | V) • (1.63) 

Note that the decomposition of the Dirac-Coulomb Green function as in (1.61) is not applicable in the 
infrared part, because of numerical problems for ultra-soft photons (infrared divergences). Rewriting 
(1.62) appropriately into a three-dimensional integral [43, 81, 82], we have 

ASma = --r7En / dz 

TT 10 TT ./n 



poo fXi 

/ dxix\ I dx2x\ Mma{x2tXi,z) . (1.64) 
Jo Ja 

The function A^MA(a;2, 2^1, ^;) is defined in analogy to the function A4m{x2,Xi,z) defined in Eq. (1.44) 
for the infrared part. Also, we define a function SMAix2, xi, z) in analogy to the function 5ir(x2, xi, z) 
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^iven in Eq. (1.50) for the infrared part, which wiU be used in Eq. (1.67) befow. We have, 

c I \ 2r^y^ fry y 

a° \ a a 



= 1 k=±|k| i,3 = l 



' \ a a y \a/ \ a aJ 



= — 2_/ '^MA,|K|('',y,z) ■ (1-65) 

" |k| = 1 

The expansion of the propagator Ga into partial waves is given in Eqs. (5.4) and (A. 20) in [81] and in 
Eqs. (D.37) and (D.42) in [4:;]. This expansion leads to the component functions G)^ ^. The terms Tmajk] 
in the last line of Eq. (1.65) read 

a / ■ \ a a / \a/ \ a a 



k=±|k.| = l 



With these definitions, the middle-energy subtraction term AEma can be written as 

AEma = - En + - / dz dy dr SMA{r,y, z) . (1-67) 
7^ 10 Jo Jo Jo 

The subtracted lower-order terms yield, 



A£;ma = - 



+ ^(v^) + FMA{nL,Za) 
200^^' 200^ ' n3 ^ ^' ' 



(i.e 



The three-dimensional integral (1.67) is evaluated by successive Gaussian quadrature. Details of the 
integration procedure can be found in [43]. The numerical results are summarized in the Table 1.2. 

1.4.3 The Middle-Energy Remainder 

The remainder term in the middle-energy part involves the propagator Gb defined in Eq. (1.39), Gb = 
G— Ga, where G is defined in (1.19) and Ga is given in (1.38). In analogy to the middle-energy subtraction 
term, the middle-energy remainder can be rewritten as a three-dimensional integral, 

A£^MB = — dz dr dy Susir^y, z) , (1.69) 

Jo Jo Jo 

where 

s„e(.-...) ^ -^E E i:/.(?)'=t(?.!-)/,(!)^n?-!) 

|k| = 1 k=±|k| i,j=l 



2 y5 



X] Tmb, \K\{r,y,z) 
|k| = 1 



(1.70) 
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z 


Fma(1Si/2, ^a) 


Fma(2Si/2, Za) 


FMA(2Pi/2,^a) 


FMA(2P3/2,^a) 


1 


2 fiQQ '^7Q Qn4 Rl'l'l 


2 720 878 SlSfl") 


083 207 '^^i(^) 


701 705 2401' 1 ~1 


2 


9 fi^Q ^fil 3S1 1 n "1 


2 681 890 660('1 


084 208 832 n 


701 850 024^1 "1 


3 


2 f!2S 77Q 45S nn^l 


2 647 262 568 fl 


085 831 658^1 "1 


702 091 147('l"l 


4 


2.591 151 010 1(1) 


2.616 290 432(1) 


0.088 040 763(1) 


0.702 426 850(1) 


5 


2.561 096 522 1(1) 


2.588 297 638(1) 


0.090 803 408(1) 


0.702 854 461(1) 


Z 


-Fmb(1Si/2, Za) 


Fmb(2Si/2, Za) 


FMB(2Pi/2,^a) 


FMB(2P3/2,^a) 


1 


1 685 993 923 2(Vi 


1 784 756 705(2) 


771 787 771f2"l 

l-I- lOI f t A-\£U 1 


—0 094 272 681(2) 


2 


1.626 842 294 5(1) 


1 725 583 798('2") 


770 778 394f2"l 


—0 094 612 071('2l 


3 


1 571 406 090 7fl") 


1 670 086 996('2") 


769 153 314f2") 


—0 095 165 248(2) 


4 


1.519 082 768 6(1) 


1.617 650 004(2) 


0.766 954 435(2) 


-0.095 922 506(2) 


5 


1.469 482 409 0(1) 


1.567 873 140(2) 


0.764 220 149(2) 


-0.096 874 556(2) 




Fm(1Si/2, Za) 


Fm(2Si/2,^q) 


FM(2Pi/2,^a) 


FM(2P3/2,Za) 


1 


4 385 373 827 7fll 


4 505 635 023('2") 


854 995 085 f2") 


607 432 559f2"l 


2 


4 286 403 675 7fll 


4 407 404 458('2") 


854 987 226f2"l 


607 237 953f2l 


3 


4 195 185 543 6(1) 


4 317 349 564('2") 


854 984 972f2"l 


606 925 899f2l 


4 


4.110 233 778 8(1) 


4.233 940 436(2) 


0.854 995 198(2) 


0.606 504 344(2) 


5 


4.030 578 931 1(1) 


4.156 170 778(2) 


0.855 023 557(2) 


0.605 979 905(2) 



Table 1.2: Numerical results for the middle-energy subtraction term Fma, the middle-energy re- 
mainder term i^MB and the middle-energy term Fyi- The middle-energy term Fn is given as the 
sum Fuinlj, Za) = FuAinlj, Za) + FuBinlj, Za) [see also Eqs.(1.25), (1.68) and (1-72)]. 



where we implicitly define the terms Ty[-Q ^^^{r,y,z) in analogy with the infrared part Eq. (1.50). The 
functions Gg ^ are obtained as the difference of the expansion of the full propagator G and the simplified 
propagator Ga into angular momenta, 

Gl^^G^^^Gl^ (1.71) 

where the Gj/ are listed in Eq. (A. 16) in [81] and in Eq. (D.43) in [4.3], and the G^ ^ have already been 
defined in Eqs. (5.4) and (A. 20) in [81] and in Eqs. (D.37) and (D.42) in [4.3]. There are no lower-order 
terms to subtract, and therefore 

a (Za)^ 

AEmb = - FMB{nlj,Za) . (1.72) 

The three-dimensional integral (1.69) is evaluated by successive Gaussian quadrature. Details of the 
integration procedure are provided in [43]. Numerical results for the middle-energy remainder -FmB; are 
summarized in Table 1.2 for the K- and L-shell states. 

For the middle-energy part, the separation into a subtraction and a remainder term has considerable 
computational advantages which become obvious upon inspection of Eqs. (1.68) and (1.72). The sub- 
traction involves a propagator whose angular components can be evaluated by recursion [43, 82], which 
is computationally time-consuming. Because the subtraction term involves lower-order components [see 
Eq. (1.25)], it has to be evaluated to high precision numerically (in a typical case, a relative accuracy 
of 10~^^ is required). This high precision requires in turn a large number of integration points for the 
Gaussian quadratures, which is possible only if the numerical evaluation of the integrand is not compu- 
tationally time-consuming. For the remainder term, no lower-order terms have to be subtracted, and the 
relative accuracy required of the integrals is in the range of 10~^^ . . . 10~^. A numerical evaluation to this 



23 



smaller level of precision is feasible although the calculation of the Green function Gb is computationally 
more expensive than that of Ga [43, 81, 82]. The separation of the high-energy part into a subtraction 
term and a remainder term, which is discussed in Sec. 1.5, is motivated by analogous considerations as 
for the middle-energy part. In the high-energy part, this separation is even more important than in the 
middle-energy part, because of the occurrence of infinite terms which need to be subtracted analytically 
before a numerical evaluation can proceed [see Eq. (1.82) below]. 

We now summarize the results for the middle-energy part. The middle-energy part is the sum of the 
middle-energy subtraction term Fma and the middle-energy remainder Fmb [see also Eq. (1.27)]. Numer- 
ical results are summarized in Table 1.2 for the K- and L-shell states. The low-energy part Fl is defined 
as the sum of the infrared contribution i^iR and the middle-energy contribution Fm [sec Eq. (1.32)]. The 
results for Fl arc provided in the Table 1.3 for the K- and L-shell states. The limits for the low-energy 
part as a function of the bound state quantum numbers can be found in Eq. (7.80) of [4-')], 

FL{nlj,Za) = ^ Si,o \n{Za)-^ - j lnfco(n, /) + (^ln2 - ^) ^ 

o , 16\ 1 , „ 7\ 1 

2 In 2 — h - In 2 - 



15 J 21 + 1 \2 A) k{21 + 1) 



3 , „ 9 
-- ln2 + - 
2 4 



The leadings asymptotics for the states under investigation arc, 

FL(lSi/2,Za) = (4/3) ln(Za)-2 _ 1.554 642 -fO(Za), 



FL(2Si/2,Za) = {A/?,)\n{Zay^ ^l.l'ilAQl + 0{Za), 



FL(2Pi/2,Za) = 0.940 023 + 0(Za), 



FL(2P3/2,Za) = 0.690 023 + 0(Za). (1.74) 

These asymptotics are consistent with the numerical data in Table 1.3. For S states, the low-energy 
contribution Fl diverges logarithmically as Za — > 0, whereas for P states, Fl approaches a constant 
as Za 0. The leading logarithm is a consequence of an infrared divergence cut off by the atomic 



z 


Fl(1Si/2,Z«) 


FL(2Si/2,Za) 


FL(2Pi/2,^a) 


FL(2P3/2,Za) 


1 


11.621 997 564 5(1) 


11.985 399 203(2) 


0.940 322 937(2) 


0.690 169 056(2) 


2 


9.825 405 794 7(1) 


10.189 430 095(2) 


0.941 060 895(2) 


0.690 517 414(2) 


3 


8.793 341 365 4(1) 


9.158 273 526(2) 


0.942 147 482(2) 


0.691 017 729(2) 


4 


8.073 357 919 4(1) 


8.439 442 234(2) 


0.943 538 386(2) 


0.691 645 132(2) 


5 


7.523 832 250 6(1) 


7.891 285 736(2) 


0.945 204 392(2) 


0.692 383 083(2) 



Table 1.3: Low-energy part Fl for the K- and L-sheU states Fl(1Si/2, Za), Fl(2Si/2, Za), Fl(2Pi/2, Za), 
and Fl(2P3/2, Za), evaluated for low-Z hydrogenlike ions. 
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momentum scale. It is a nonrelativistic effect which is generated by the nonvanishing probability density 
of S waves at the origin in the nonrelativistic limit. The presence of the logarithmic behavior for S 
states [nonvanishing A41 -coefficient, see Eqs. (1.2) and (1.3)] and its absence for P states is reproduced 
consistently by the data in Table 1.3. 

1.5 The High-Energy Part 

1.5.1 The High-Energy Subtraction Term 

The high-energy part is given by 

X (V' I e"=-=^ GiEn - lu) a" e-""'^ | V) , (1.75) 

where relevant definitions of the symbols can be foimd in Eqs. (1.18)-(1.21), and the contour Ch is as 
shown in Fig. 1.1. The high-energy part comprises the following integration region for the virtual photon, 

uj G {En — i 00, £"„ + i 00) I 

> high-energy part AEn . (1-76) 
z e (— ioo, ioo) J 

The separation of the high-energy part into a subtraction term and a remainder is accomplished as in 
the middle-energy part [see Eq. (1.61)] by writing the full Dirac-Coulomb Green function G [Eq. (1.19)] 
as G = Ga + Gb- We define the high-energy subtraction term F^a as the expression obtained upon 
substitution of the propagator Ga for G in Eq. (1.75), and a substitution of the propagator Gb for G in 
Eq. (1.75) leads to the high-energy remainder ^hb which is discussed in Sec. 1.5.2. The subtraction term 
(including all divergent contributions) is generated by Ga, the high-energy remainder term corresponds 
to Gb. We have 

9 /" duo f d^k „ 9 . ^ 
Ai?HA = - lim ie^ / — -—D^,{k\K) 



(V' 



a^e''=-^ GA{En-oj)a'' c 



— ifc-; 



(1.77) 



and 



HB 



-le 



du f d'^k 



Ch J (2^) 



X (V I a'' e''^-^ Gb(S„ - uj) a" c-"=-^ | ■ 



(1.78) 



The contribution AEha corresponding to Ga can be separated further into a term AE^J^, which contains 
all divergent contributions, and a term AE^j^, which comprises contributions of lower order than (Za)^, 
but is convergent as A ^ 00. This separation is described in detail in [81, n4]. We have 



A^HA - AeI^I + A^^'l . 



(1.79) 



We obtain for AE^^]^, which contains a logarithmic divergence as A — > cxd. 



AE, 



(1) _ 

HA ~ 



I m A^ - I) iP) + ln2 - 11) (V) + F«(./„ Za) 



(1.80) 
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For the contribution Pj^J^, an explicit analytic result is obtained in Eq. (4.15) in [s]]. This contribution 
is therefore not discussed in any further detail, here. The contribution AE^j^ contains lower-order terms, 



Altogether we have 

AEha 



(2) 
HA 



a 








(-il„2 4 


4) 









= AE, 



HA 



AE, 



(2) 
HA 



(V) 



F^i{nl,,Za) 



(1.81) 



a 

TT 



^ In - ^) (/?) ^1{V) + ^ i^HA(n/„ Za) 



The scaled function Fn^inlj, Za) is given as 



(1.82) 



FnAinlj, Za) = F^'linlj, Za) + F^l{nlj,Za) 



(1.83) 



The term AE^j^ falls naturally into a sum of four contributions [81], 



AE^^l = Ti + + T3 + r4 



(1.84) 



where 



To 



l^{V)+^^h,{nl,,Za), 
7 1 



20 2 



In 2) {V)+''-^^h2{nl,,Za), 



T3 = ^^—^h'i{nlj,Za), 



{Za)'^ , , , 

r — hAnlj^ Za) . 



(1.85) 



The functions hi {i = 1,2,3,4) are defined in Eqs. (4.18), (4.19) and (4.21) in [81] (see also Eq. (3.6) 
in [83]). The evaluation of the high-energy subtraction term proceeds as outlined in [81, 82, 83], albeit 
with an increased accuracy and improved calculational methods in intermediate steps of the calculation 
in order to overcome the severe numerical cancellations in the \ow-Z region. We recover F-^J^ as the sum 



{2) 

F^J^ijilj, Za) = hi{nlj, Za) + h2{nlj, Za) + h^lnlj, Za) + h4{nlj, Za) . 
The scaled function F^Ainlj, Za) [sec also Eqs. (1.26) and (1.28)] is obtained as 

FiiA{nlj,Za) = F^linlj, Za) + F^l{nlj, Za) . 



(1.86) 



(1.87) 



The limits of the contributions F^J^lnlj, Za) and F^p^(nlj, Za) as {Za) — > have been investigated 

in [81, 83, 84]. For the contribution Fj^p^{nlj, 0), the result can be found in Eq. (3.5) in [83]. The limits of 
the functions hiinlj, Za) (i ~ 1,2,3,4) as Za — > are given as a fimction of the atomic state quantum 
numbers in Eq. (3.8) in [83]. For the scaled high-energy subtraction term i^HA, the limits read (see 
Eq. (3.9) in [83]) 



(2), 



FiiA{nlj,Za) = 



li_ln2U 
10 / n 



i^-21n2)-^ 
15 J 21 + 1 



2 4 



1 



k(2; + 1) 
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Therefore, the exphcit forms of the Umits for the states under investigation are, 

i^HA(lSi/2,^a) = -1.219 627+ 0(Za) , 



FHA(2Si/2,^a) 



1.423 054+ 0(Za) , 



Fha{2Pi/2, Za) 



1.081 204+ 0(Za) , 



FHA(2P3/2,^a) 



0.524 351 + 0(Za). 



(1.89) 



Numerical results for ^ha, which arc presented in Table 1.4, exhibit consistency with the limits in 



z 


^ha(1Si/2, Za) 


FHA(2Si/2,Za) 


FHA(2Pi/2,^a) 


FHA{2P3/2,Za) 


1 


-1.216 846 660 6(1) 


-1.420 293 291(1) 


-1.081 265 954(1) 


-0.524 359 802(1) 


2 


-1.214 322 536 9(1) 


-1.417 829 864(1) 


-1.081 451 269(1) 


-0.524 385 053(1) 


3 


-1.212 026 714 1(1) 


-1.415 635 310(1) 


-1.081 760 224(1) 


-0.524 427 051(1) 


4 


-1.209 942 847 4(1) 


-1.413 693 422(1) 


-1.082 192 995(1) 


-0.524 485 727(1) 


5 


-1.208 059 033 6(1) 


-1.411 992 480(1) 


-1.082 749 845(1) 


-0.524 561 017(1) 


Z 


-Fhb(1Si/2, Za) 


Fhb(2Si/2, Za) 


FHB(2Pi/2,^a) 


Fhb(2P3/2,^«) 


1 


-0.088 357 254(1) 


-0.018 280 727(5) 


0.014 546 64(1) 


-0.042 310 69(1) 


2 


-0.082 758 206(1) 


-0.012 729 99(1) 


0.014 574 21(1) 


-0.042 296 81(1) 


3 


-0.076 811 229(1) 


-0.006 861 02(1) 


0.014 620 51(1) 


-0.042 273 58(1) 


4 


-0.070 590 991(1) 


-0.000 746 40(1) 


0.014 685 82(1) 


-0.042 240 92(1) 


5 


-0.064 146 139(1) 


0.005 567 16(1) 


0.014 770 52(1) 


-0.042 198 76(1) 


Z 


FH(lSi/2,Za) 


FH(2Si/2,Za) 


FH(2Pi/2,Za) 


FH(2P3/2,Za) 


1 


-1.305 203 915(1) 


-1.438 574 018(5) 


-1.066 719 31(1) 


-0.566 670 50(1) 


2 


-1.297 080 743(1) 


-1.430 559 85(1) 


-1.066 877 06(1) 


-0.566 681 86(1) 


3 


-1.288 837 943(1) 


-1.422 496 33(1) 


-1.067 139 72(1) 


-0.566 700 63(1) 


4 


-1.280 533 839(1) 


-1.414 439 82(1) 


-1.067 507 18(1) 


-0.566 726 65(1) 


5 


-1.272 205 173(1) 


-1.406 425 32(1) 


-1.067 979 33(1) 


-0.566 759 78(1) 



Table 1.4: Numerical results for the high-energy subtraction term ^ha and the high-energy remainder 
term i^HB- The high-energy term Fn is given as the sum Fii{nlj, Za) = FnAin-lj, Za) + FuBinlj-, Za). 



1.5.2 The High-Energy Remainder 

The remainder term in the high-energy part involves the propagator Gb defined in Eq. (1.39), Gb = 
G — Ga, where G is defined in (1.19) and Ga is given in (1.38). The evaluation proceeds in complete 
analogy to the calculation of the middle-energy remainder term (Sec. 1.4.3). The only difference lies in 
the different integration region for the photon energy, which is given - for the high-energy part - in 



Eq. (1.89). 
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Eq. (1.76). The photon energy mtegration variable z is conveniently expressed as 



where 



(1.90) 



The method of integration is described in [12, 43], and we do not discuss any further details, here. We 
focus instead on the convergence acceleration technique used in the evaluation. 

In analogy with the middle-energy remainder, we may write the integrand S'hb which is defined in 
complete analogy to (1.70) as a sum over angular momenta ("partial waves") 



S'hb oc ^ 7hb,|k1 
|«|=i 



(1.91) 



where the 7hb,|k| are defined in analogy to Eq. (1.70). Here, |k| represents the modulus of the Dirac 
angular momentum quantum number of the virtual intermediate state. The asymptotic behaviour of the 
Thb.IkI for large \k\ is [see Eq. (4.7) in [82]] 



T, 



r.2 l«l 



HB,|k| 







const. - 









(1.92) 



where "const." is independent of The series in Eq. (1.91) is slowly convergent for r close to one, and 
the region near r = 1 is known to be problematic in numerical evaluations. 

It is found that the convergence of the scries (1.92) series near r = 1 can be accelerated very efficiently 
using the combined nonlinear-condensation transformation (see [l':!] and Sec. 7.1.6) applied to the series 
Sfc°=o^fc where ifc = iHB.fe+i [see Eqs. (1.91) and (1.92)]. The combined transformation (combination of 
the condensation transformation and the Weniger transformation) was found to be applicable to a wide 
range of slowly convergent monotone series (series whose terms have the same sign) , and many examples 
for its application were given in Ref. [23, 31]. For the numerical treatment of radiative corrections in low-Z 
systems, the transformation has the advantage of removing the principal numerical difficulties associated 
with the slow convergence of angular momentum decompositions of the propagators near their singularity 
for equal radial arguments. 

All that remains to be discussed in the current section is the low-Z limit of the energy shift 



Aii;HB = -^^i^HB(r^?„^«). 

For the high-energy remainder ^hb, the limits as Zol read [see Eq. (4.15) in [83]] 

1 



(1.93) 



Fw,{nlj, Za) 



21 + 1 



17 4 



^-21n2U 
2 K 





n~2l-l' 








n 



0{Za). 



For the atomic states under investigation, this leads to 

^HB(lSi/2,Za) = -OmS A57 +0{Za) , 



(1.94) 



FHB(2Si/2,^a) 



-0.023 364 + 0(Za) 



^HB(2Pi/2,^a) = 0.014 538 + 0(Za) 



-FHB(2P3/2,2'a) = -0.042 315 + 0(Za). 



(1.95) 
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1.5.3 Results for the High-Energy Part 

The limit of the function Fa as Za can be derived easily from the Eqs. (1.88), (1.94) as a function 
of the bound state quantum numbers. It reads 

Fii(nl,,Za) = I — -ln2 I - + I — - 2 ln2 ) — — 
^ \10 J n \15 J 21 + 1 

/ 3 , „ 5\ 1 /3 , „ 9\ 1 

+ (S - 5 ^) *•« + - ^) ^TT^ + «(^"' ■ 

For the atomic states investigated here, this expression yields the numerical values, 

Fii{lSi/2,Za) = -1.313 085 + 0(Zq;), 

FH(2Si/2,Za) = -1.446 418 + 0(Za), 

FH(2Pi/2,Za) = -1.066 667 + 0(Za), 

FH(2P3/2,Za) = -0.566 667 +0(Za). 

Numerical results for the high-energy part 

Fii{nlj, Za) = FiiA{nlj, Za) + Fhb("-'j, Za) 

are summarized in Table 1.4. Note the apparent consistency of the numerical 
their analytically obtained low-Z limits in Eq. (1.97). 

1.6 Comparison to Analytic Calculations 

The numerical results for the scaled function F{nlj,Za) describing the self energy defined in Eq. (1.1) 
are given in Table 1.5, together with the results for the nonperturbative function GsEinlj, Za), which is 
implicitly defined in Eq. (1.2). Results are provided for K and L shell states. The numerical results for 
the remainder Gse are obtained by subtracting the analytic lower-order terms listed in Eq. (1.2) from 
the complete numerical result for the scaled function F{nlj, Za). No additional fitting is performed. 

Analytic and numerical results at low Z can be compared by considering the remainder function Gse • Note 
that an inconsistency in any of the analytically obtained lower-order terms would be likely to manifest 
itself in a grossly inconsistent dependence of GsEinlj, Za) on its argument Za; this is not observed. 
For S states, the following analytic model for Gse is commonly assumed, which is motivated in part 
by a renormalization-group analysis [S7] and is constructed in analogy with the pattern of the analytic 
coefficients Aij in Eq. (1.2) and (1.3), 

GsEinSi/2, Za) = AGo(nSi/2) + (Za) [Aji{nSi/2) lri{Za)^^ + Aro{nSi/2)] 
+ {Za)^ [As3{nSi/2) \n\Za)-^ + As2{nSi/2) In^Za)'^ 
+A8i(nSi/2) ln(Za)-2 + A8o(nSi/2)] . (1.99) 



(1.97) 
(1.98) 

results in Table 1.4 with 
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The (probably nonvanishing) Ag^ coefficient, which introduces a triple logarithmic singularity at Za ~ 0, 
hinders an accurate comparison of numerical and analytic data for Gse- A somewhat less singular behavior 
is expected of the difference 

AGsEiZa) = Gse(2Si/2, Za) - Gse(1Si/2, Za) , (1.100) 

because the leading logarithmic coefficients in any given order of Za are generally assumed to be equal 
for all S states, which would mean in particular 

An(lSi/2) = A7i(2Si/2) (1.101) 

and 

A83(lSi/2) = ^83(2Si/2). (1.102) 

Now we define AAki as the difference of the values of the analytic coefficients for the two lowest S states, 

AAu - Afc,(2Si/2) - v4fci(lSi/2) . (1.103) 

The function AGse defined in Eq. (1.100) can be assumed to have the following semi-analytic expansion 
around Za = 0, 



AGsE(^a) = AAeo + {Za) AA^ + [Za)^ [AAg2 \^{Z^ 



a) 2 



-fAAgi In(Za)"^ + AAgo + o{Za)\ . (1.104) 

In order to detect possible inconsistencies in the numerical and analytic data for Gse, we difference the 
data for AGse, i-e-, we consider the following finite difference approximation to the derivative of the 
function AGse, 

giZ) = AGse((^ + 1) a) - AGse(^«) ■ (1-105) 

It is perhaps not a priori obvious why this combination leads to a sensible comparison of numerical and 
analytic data. This will be explained in the sequel. 

We denote the analytic and numerical limits of AGse(^Q;) as Za ^ as Aylgg** and A^g™, respectively. 
Of course, we hope that these will turn out to be equal. However, we temporarily leave open the possibility 
of an inconsistency between numerical and analytic data by keeping A^g™ and A^Iq'* as distinct 
variables. In order to illustrate how a discrepancy could be detected by investigating the function g{Z), 
we consider special cases of the function AGse('^Q;) and g{Z). For Z = 0, the exact result can be inferred 
exclusively and uniquely using the analytic approach of Ref. [-32], and we have 

AGse(O) = A^-- , (1.106) 
whereas for Z = 1, which is determined by numerical data, 

AGsE(a) « AA^r + " ^^70 , (1.107) 
with a possibly different limit A^gQ™ as Z — ^ 0, and for Z — 2, 

AGsE(2a)« AAgJ"^ + 2aA^7o, (1-108) 
where we ignore higher-order analytic terms. Hence, for Z = 0, we have 

g{0) = AGse(«) - AGse(O) = AA'™" - AA^- + a [AAjo + o{Za)] . (1.109) 

A possible inconsistency (nonvanishing AAgp™ — AA||}^) would influence the value of ^(0). For Z = 1, 
the value of g is determined solely by numerical data, 

5(1) = AGsE(2a) - AGsE(a) = a [AA^ + o^Za)] , (1.110) 

and for Z = 2, we have 

g(2) = AGsE(3a) - AGsE(2a) = a [AA^o + o(Za)] . (1.111) 
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Of course, analogous relations hold for Z > 2. This means that a meaningful comparison of the analytic 
and numerical approaches can be made by comparing the value 5(0), which is influenced by the term 
AAqq'^-AAqq'', to the other values 5(1), . . . , g(4), which are independent of the difference AAgo™-A^|o''. 
We recall that the numerical data from Table 1.5 lead to the evaluation of the five values (?(0), .g(l), .9(2), 
5(3) and (?(4). A plot of the function g{Z) serves two purposes: (i) the values g{l), ■ ■ ■ ,5(4) should exhibit 
apparent convergence to some limiting value aAA^Q as Z 0, and this can be verified by inspection 
of the plot, and (ii) a discrepancy between the analytic and numerical approaches - as explained above 
- would result in a nonvanishing value for AAg™ — AAIJ}'', and it would introduce an inconsistency 
between the trend in the values of 5(1), . . . , <7(4), and the value of g{0) [see Eq. (1.109)]. 



z 


F(lSi/2,Za) 


F{2S,/2,Za) 


F(2Pi/2,Za) 


F(2P3/2,Z«) 


1 


10.316 793 650(1) 


10.546 825 185(5) 


-0.126 396 37(1) 


0.123 498 56(1) 


2 


8.528 325 052(1) 


8.758 870 25(1) 


-0.125 816 16(1) 


0.123 835 55(1) 


3 


7.504 503 422(1) 


7.735 777 20(1) 


-0.124 992 24(1) 


0.124 317 10(1) 


4 


6.792 824 081(1) 


7.025 002 41(1) 


-0.123 968 79(1) 


0.124 918 48(1) 


5 


6.251 627 078(1) 


6.484 860 42(1) 


-0.122 774 94(1) 


0.125 623 30(1) 


Z 


Gse(1Si/2, Za) 


G'se(2Si/2, Za) 


GsE(2Pi/2,Za) 


G'se(2P3/2, Za) 


1 


-30.290 24(2) 


-31.185 15(9) 


-0.973 5(2) 


-0.486 5(2) 


2 


-29.770 967(5) 


-30.644 66(5) 


-0.949 40(5) 


-0.470 94(5) 


3 


-29.299 170(2) 


-30.151 93(2) 


-0.926 37(2) 


-0.456 65(2) 


4 


-28.859 222(1) 


-29.691 27(1) 


-0.904 12(1) 


-0.443 13(1) 


5 


-28.443 372 3(8) 


-29.255 033(8) 


-0.882 478(8) 


-0.430 244(8) 



Table 1.5: Numerical results for the scaled function F [defined in Eq. (1.1)] and the remainder function 
GsE implicitly defined in Eq. (1.2). 



0.0220 



0.0215 



0.0210 



0.0205 



0.0200 

12 3 4 

Z 

Figure 1.4: Plot of the function g{Z) defined in Eq. (1.105) in the region of 
low nuclear charge. For the evaluation of the data point at Z = 0, a value of 
^6o(lSi/2) = -30.92415(1) is employed [9, 32, 88]. 



I 



I 



I 



J L 



31 



Among the different evaluations of Aqq for the ground state, the result in [■■)'.!] has the smallest quoted 
uncertainty. In Fig. 1.4 we display a plot of g{Z) for low nuclear charge Z. A value of ^60(181/2) = 
A|[J^(lSi/2) = -30.92415(1) [9, 32, 88] is used in Fig. 1.4. The results indicate very good agreement 
between the numerical and analytic approaches to the Lamb shift in the \ow-Z region up to the level of 
a few Hz in frequency units for the low- lying atomic states (where n is the principal quantum number). 
The error bars represent the numerical uncertainty of the values in Table 1.5, which correspond to an 
uncertainty on the level of 1.0 x Z*Hz in frequency units. 

Analytic work on the correction A^o has extended over three decades [32, 68, 69, 70, 71]. The complication 
arises that although the calculations are in general analytic, some remaining one-dimensional integrations 
could not be evaluated analytically because of the nature of the integrals [sec e.g. Eq. (6.96) in [32]]. 
Therefore a step- by-step comparison of the analytic calculations is difficult. An additional difficulty is the 
isolation of those analytic terms which contribute in a given order in Za, i.e. the isolation of only those 
terms which contribute to Aqq. The apparent consistency of the numerical and analytic data in Fig. 1.4 
represents an independent consistency check on the rather involved analytic calculations. 
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Figure 1.5: Comparison of numerical data and analytically evaluated higher- 
order binding corrections for the 2Pi/2 state. We plot the function gi/2{Z) 
defined in Eq. (1.114) in the region of low Z. The numerical data appear to 
be consistent with the analytic result of ^6o(2Pi/2) = —0.99891(1) obtained 
in [1]. 



Our numerical results are not inconsistent with the analytic result [44] for a higher-order logarithm, 

A71 -TT (^^-ln2^ =4.65, (1.112) 

although they do not necessarily confirm it. As in ['.)], we obtain as an estimate A71 = 5.5(1.0) (from the fit 
to the numerical data for both S states). Logarithmic terms corresponding to the probably nonvanishing 
As3 coefficient should be taken into account for a consistent fit of the corrections to Gse- These highly 
singular terms are difficult to handle with a numerical fitting procedure. The terms Aga, Ag2 and Agi 
furnish three more free parameters for the numerical fit, where only five data points are available (in 
addition to the quantities Aqq, Ah and A-^q, which may also be regarded as free parameters for the 
fitting procedure). The determination of Aqq by a fit from the numerical data is much more stable 
than the determination of the logarithmic correction A71 . We briefly note that our all-order evaluation 
essentially eliminates the uncertainty due to the unknown higher-order analytic terms. It is interesting 
to note that the same numerical methods are employed for both the S and P states in our all-order (in 
Za) calculation, whereas the analytic treatment of S and P states differs [1, 32]. 

The comparison of numerical and analytic result is much less problematic for P states, because the 
function Gse is less singular [see Eqs. (1.4) and (1.6)]. For the 2P states, we observe that the function 
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GsE{'2Pj,Za) has the same semi- analytic expansion about Za = as the function AGse(-Z^Q^) defined 
for S states in Eq. (1.100). We have 

Gse(2Pj, Za) = A6o(2P,) + (Za) ^7o(2Pj) + (Zaf [^82(2?^) ln\Za)-^ 

+A8i(2Pj) ln(Za)-2 + A8o(2P,) + o{Za)] . (1.113) 

Hence, we plot the function 

gj{Z) - Gse(2P„ (Z + 1) a) - Gse(2Pj, Za) (1.114) 

for j = 1/2 in the region of low Z, with the notion that an inconsistent analytic result for yl6o(2Pj) would 
lead to irregularity at Z = 0, in analogy with the S states. The numerical data shown in Fig. 1.5 appears 
to be consistent with the analytic results of 

Ao(2Pi/2) = -0.998 91(1). (1.115) 

obtained in [I]. For the 2P3/2 state, see Fig. 6 of Ref. [\'2]. In this context it may be interesting to note 
that analytic results obtained in [1, 2] for the higher-order binding corrections to 2P, 3P and 4P states 
have recently been confirmed indirectly [89]. 
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Chapter 2 



Analytic Self— Energy Calculations 
for Excited States 

2.1 Orientation 

Analytic calculations of the one-loop self energy in bound systems have a long history, starting from 
Bethe's seminal paper [■!>!], and have meanwhile extended over more than five decades. The calculations 
arc based on an expansion of the effect into powers of the parameter Za, where Z is the nuclear charge 
and a is the fine-structure constant (see e.g. [f , 2, 32, 43]). The expansion is semi-analytic, i.e. it involves 
powers of Za. and ln[(ZQ;)~^]. Terms which are of higher than leading order in Za. arc commonly referred 
to as the binding corrections (see also Ch. 3). A certain nonlogarithmic binding correction, known as the 
Ago-coefficient, has been of particular interest. For relevant formulas and definitions regarding ^eo: we 
refer to Ch. 1, especially to Eqs. (1.1) - (1.6). 

In this chapter, we present results of an analytic evaluation of higher-order binding corrections to the 
one-loop self energy for excited bound states of hydrogenlike systems. We focus on states with principal 
quantum numbers n = 2 ... 5 and Dirac angular quantum numbers k = 1, —2, 2, . . . , 4, —5. The current 
calculation represents a continuation of previous investigations on the subject [1, 2]. 

The improved results for Ago coefficients rely essentially on a more compact code for the analytic cal- 
culations, written in the computer-algebra package Mathematica [90], which enables the corrections to 
be evaluated semi-automatically. Intermediate expressions with some 10,000 terms are encountered, and 
the complexity of the calculations sharply increases with the principal quantum number n, and, as far 
as the complexity of the angular momentum algebra is concerned, with the orbital angular momen- 
tum I = \k + 1/2| — 1/2 of the bound electron. Calculations were performed, in part, on IBM RISC/6000 
workstations, both at the National Institute of Standards and Technology and at the Laboratoire Kastler- 
Brossel, and profited from the availability of 2 GigaBytes of Random Access Memory. 

Of crucial importance was the development of convergence acceleration methods which were used ex- 
tensively for the evaluation of remaining one-dimensional integrals which could not be done analytically. 
These integrals are analogous to expressions encountered in previous work [1, 2]. The numerically eval- 
uated contributions involve slowly convergent hypergeometric series, and - in more extreme cases - 
infinite series over partial derivatives of hypergeometric functions, and generalizations of Lerch's $ tran- 
scendent [91, 92]. We have found the Combined Nonlinear-Condensation Transformation [23] to be a 
versatile method for the efficient numerical evaluation of the slowly convergent series encountered in the 
evaluation of the AgQ-coefficients. 

Results are given below for the AgQ-coefficicnts with an absolute precision of 10^^. As explained in 
detail in [1, 2, 32], the calculation of the one- loop self energy falls naturally into a high- and a low- 
energy part [Fh and F^, respectively). In turn, the low-energy part can be separated naturally into 
the nonrelativistic dipole and the nonrelativistic quadrupole part, and into relativistic corrections to 
the current, to the hamiltonian, to the wavefunction, and to the energy of the bound state. For some 
states, severe numerical cancellations are observed between these (low-energy) contributions to the low- 
energy part of the self energy. In the final step of the calculation, we observe additional numerical 
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cancellations between the high- and the low-energy part for the finite contributions to Aqq, beyond the 
necessary cancellation of the divergent contributions which depend on the scale-separation parameter e; 
this parameter and the mathematical method which is employed in the expansion of the self energy are 
described and explained in Sec. 2.2. For some of the investigated atomic states, the absolute magnitude 
of the ^60-coefhcients is as small as 10~^, although the largest single contribution is always on the order 
of 10~^. This illustrates the cancellations between the finite contributions to Aeo- 



2.2 The epsilon-Method 

We discuss here, by way of example, the e method employed in the analytic calculation of the self energy 
in bound systems. This method is very suitable [1, 2, 32] for the separation of the two different energy 
scales for virtual photons: the nonrclativistic domain, in which the virtual photon assumes values of the 
order of the atomic binding energy, and the relativistic domain, in which the virtual photon assumes 
values of the order of the electron rest mass. Different approximation schemes and different asymptotic 
expansions are adequate for the two different domains. Without these approximations and expansions, 
the analytic evaluation of either the high- or the low-energy part would not be feasible. Therefore, the 
e method is also used for the treatment of the highly excited atomic states which are the subject of 
(this part) of the thesis. At the same time, the model example discussed here is meant to illustrate the 
usefulness of the "e method" in a more general context. 

We will consider here a model problem with one "virtual photon" . The separation into high-and low- 
energy photons necessitates the temporary introduction of a parameter e; the dependence on e cancels 
when the high- and the low-energy parts are added together. We have, 

nonrclativistic domain <C e <C electron rest mass (2.1) 



(Za)^mo <e< , (2.2) 

where a is the fine structure constant, and Z is the nuclear charge. The high-energy part is associated 
with photon energies > e, and the low-energy part is associated with photon energies uj < e. 

In order to illustrate the procedure, we discuss a simple, one-dimensional example: the evaluation of 



where the integration variable lo might be interpreted as the "energy" of a "virtual photon" . The integral 
/ can be expressed in terms of special functions. 



where E is the complete elliptic integral of the second kind, and 2F1 denotes a hypergeometric function. 
An alternative integral representation reads /(/3) = J^^'^ \J P'^ ^ sin^(ijj) dw. 

The purpose of the calculation is to derive a semi-analytic expansion of /(/3) in powers of /3 and In (3. 
The fine structure constant a takes the role of the expansion parameter /3 in actual calculations of the 
self energy. We discuss first the "high-energy part" of the calculation. It is given by the expression 



For w > e, we may expand 



v/^^^^+T?^ = ^ + 1^ + + 0(/3') , (2.6) 

but this expansion is not applicable in higher orders to the domain < < e because of the appearance 
of inverse powers of cj (analogous to an "infrared divergence" in QED). 

The separation parameter e acts an infrared regulator. After expanding in (3 [see Eq. (2.6)], the resulting 
integrals in each order of (3 can be evaluated analytically. Subsequently, we expand every term in the 
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/9-expansion in powers of e up to the order e*^, i.e. we keep only the divergent and constant terms in e. 
The result is 

(2.7) 

Here, the "C'-symbol identifies a contribution for which 0{x)/x — *■ const, as a; ^ 0, whereas the "o"- 
symbol identifies the weaker requirement o(a;) ^ as a; ^ 0; this is consistent with the standard notation 
(see e.g. [53]). 

The contribution /h (/?) corresponds to the "high-energy part" in an analytic calculation of the self energy, 
where the propagator of the bound electron is explicitly expanded in powers of the fine structure constant 
a. Now we turn to the "low-energy part" . The expression for the low-energy part (0 < a; < e) reads 



The expansion (2.6) is not applicable in this energy domain; we therefore have to keep the numerator of 
the integrand ^ uP- ^ f3'^ in uncxpanded form. However, we can expand the denominator \/l — uj"^ of the 
integrand in powers of w; because < a; < e (with e small), this expansion in lu is in fact an expansion 
in /3 - although the situation is somewhat problematic in the sense that every term in the w-expansion 
gives rise to terms of arbitrarily high order in the /3-expansion [see also Eq. (2.10) below]. 



The term y^uj^ + f]^ is analogous to the Schrodinger-Coulomb propagator in the evaluation of the self 
energy - it has to be kept in an uncxpanded form ~, whereas the expansion 

^ =1 + ^ + + (2.9) 



VT^7J2 2 

corresponds to the expansion into the (ZQ;)-expansion in the low-energy part. 

Every term in the expansion (2.9) gives rise to arbitrarily high-order corrections in /?, but it starts with 
the power For example, we have for the leading term of order a;° = 1 from Eq. (2.9), 

^'V^J^T^dc.^/?^ |ilnge^+i + 0(,)| 

+'^' { 1^ + ^(^)} + {"61? + ^(^)} + '^^^'^ ■ ^^-^^^ 

Note that the terms generated in the orders (3^ and are needed to cancel divergent contributions in 
respective orders of /3 from the high-energy part given in Eq. (2.6). The term of order uj'^ from (2.9) 
results in 

i^^^V^d. ^ H In (|.) +! + (.(.)} 



+/3' i~-^ + Oie)}+OiP^). (2.11) 
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Altogether, we obtain for the low-energy part, 

1 . /2 \ 1 



16 £2 16 \f] 



64 



1 1 



64 £4 64 £2 128 \P J 512 



Oie] 



In/?). 



(2.12) 



When the high-energy part (2.7) and the low-energy part (2.12) are added, the dependence on e cancels, 
and we have 

/(/?) = /H(/3,£)+/L(Ae) 



= 1 /3M - In 



/3 



+ 



16 \(3 J 64 



128 \PJ 128 J 



©(Z?*^ ln/3). 



(2.13) 



In order to illustrate the analogy with the self energy calculation presented here, we would like to point 
out that the dependence on e cancels out in the final result which is the sum of the high-energy part Eu 
given in Eq. (2.14) [for the two-loop case see Eq. (3.42)] and the low-energy part i?L in Eq. (2.15) [for 
the two-loop case sec Eq. (3.58)]. 



2.3 Results for High— and Low— Energy Contributions 



2.3.1 P_l/2 states (kappa = 1) 

Here and in the following subsections, we present results which have been obtained for the high- and 
low-energy energy parts of the self energy or all states with principal quantum number n ~ 2 ... 5 
and orbital angular quantum number I > 0. These results, combined with numerical calculations for 
low-, intermediate and high nuclear charge, allow to reliably estimate quantum electrodynamic (QED) 
radiative corrections for any atomic state with principal quantum number n < 5. For the definition of 
the scaled self energy function F and of the analytic A-coefficicnts, see Eqs. (1.1) - (1.6). The results are 
given in terms of the expansion parameter Za and the scale-separation parameter e, whose role in the 
calculation is explained in Sec. 2.2. 

For Pi/2 states, we obtain for the following results for the high- and low-energy parts. The 2Pi/2 high- 
energy part is: 



FHi2P,/2)^-l + {Zaf 



4177 2 103 , , 

ln(2€) 

21600 9£ 180 ^ ' 



2Pi/2 low-energy part: 

i^L(2Pi/2) - -| hifco(2P) + (Za)' 



2 103 

-0.795 649 812(1) + - + - m^^^^^. 



(2.14) 



(2.15) 
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3Pi/2 high-energy part: 



i^H(3Pi/2) = -g + (Za)' 



6191 20 268 , , 
— ln(2e) 



24300 81 e 405 



(2.16) 



3Pi/2 low-energy part: 



i^L(3Pi/2) = --lnfco(3P) + (^a)' 



^0.944 288 447(1) + ^ + ^lnf ' 



81 e 405 \iZa)^ 



(2.17) 



4Pi/2 high-energy part: 



^^h(4Pi/2) = -g + (^«)' 



24409 23 499 



86400 90 e 720 



hi(2e) 



(2.18) 



4Pi/2 low-energy part: 



i^L(4Pi/2) = --hrA:o(4P) + (Zaf 



23 499 / e 
-0.997 810 211(1) + - +-hr(^^ 



5Pi/2 high-energy part: 



(2.19) 



20129 292 796 



67500 1125 e 1125 



ln(2e) 



(2.20) 



5Pi/2 low-energy part: 



i^L (5Pi/2 ) = - 3 In fco (5P) + (Za) 



2Q2 796 / f 
_1^023 991781,l) + — +— ln(^ 



(2,21) 



Here and in the following, In fco is the Betlie logarithm for which very accurate values are available [72, 
73, 74, 75, 76, 77, 79]. 

2.3.2 P_3/2 states (kappa = -2) 

The obtained results read as follows. The 2P3/2 high-energy part reads: 



6577 2 29 , / 



(2.22) 



2P3/2 low-energy part: 



Fi(2P3/2) = --lnfco(2P) + (Z«)2 



2 29 / e 
-0.584 516 780(1) + - + - ln(^^ 



3P3/2 high-energy part: 



(2.23) 



FH{3Py2) = j^ + {Zaf 



67903 20 148 



194400 81 e 405 



ln(2e) 



(2.24) 



3P3/2 low-energy part: 



^l(3P3/2) = --lnfco(3P) + (Za)'' 
4P3/2 high-energy part: 



-0.693 566 427(1) + — + —ln( y-^ 
^ ' 81 e 405 \iZa) 



(2.25) 



i^ff(4P3/2) - + (Za)' 



31399 23 137 
86400 ~ 90e " 360 



(2.26) 



38 



4P3/2 low-energy part: 

^^L(4P3/2) = -tlnfco(4P) + (^c 



5P3/2 high-energy part: 



-0.730579 137(1) + ^ + HI l„(^ 



(2.27) 



199387 292 436 



540000 1125 e 1125 



ln(2e) 



(2.28) 



5P3/2 low-energy part: 



ft(5P3/2) = --lnto(5P) + (Zaf 



_„.™M5 653(l) + ^ + ii|,„(^ 



(2.29) 



2.3.3 D_3/2 states (kappa = 2) 

The results obtained for states with n = 2 read as follows. 3D3/2 high-energy part: 



Fh(3D3/2) = + (^«)2 



1721 4 4 



194400 405 e 405 



ln(2e) 



(2.30) 



3D3/2 low-energy part: 



i^L(3D3/2) = --lnA:o(3D) + {Zaf 



0.021 250 354(1) + + 7^ In (7^ 

405 e 405 \{Za) 



4D3/2 high-energy part: 



(2.31) 



i^i/(4D3/2) = - ^ + (^«)' 



829 1 1 



86400 90 e 90 



ln(2e) 



(2.32) 



4D3/2 low-energy part: 



^l(4D3/2) = --lnfco(4D) + [Zaf 



0.022 882 528(1) + ^ + ! ln(^ 



5D3/2 high-energy part: 



^^ff(5D3/2) = + {Z^^f 



35947 92 92 



3780000 7875 e 7875 



ln(2e) 



5D3/2 low-energy part: 



(2.33) 



(2.34) 



^^l(5D3/2) = --lnfco(5D) + [Zaf 



92 9*5 / e 
0.023 759 683(1) + ^ + ^ ln(^ 



(2.35) 



2.3.4 D_5/2 states (kappa = -3) 

For D5/2 states with Dirac quantum number k = —3, the results read as follows. For the 3D5/2 high-energy 
part, we have: 



371 4 4 , , 

In 2e 

24300 405 e 405 ^ ' 



(2.36) 



3D5/2 low-energy part: 



^l(3D5/2) = --lnA:o(3D) + {Zaf 



0.019 188 397(1) + + — In f — ^ 
^ ' 405 e 405 \{Za 



(2.37) 
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4D5/2 high-energy part: 



53 1 1 



2880 90 e 90 



hi(2e) 



4D5/2 low-energy part: 



^^l(4D5/2) = --lnA:o(4D) + {Zaf 



0.020 710 720(1) + ^ + ! ln(^ 



5D5/2 high-energy part: 



FH{^D^/2) = ^ + [Zaf 



3097 92 92 



157500 7875 e 7875 



ln(2e) 



5D5/2 low-energy part: 

Fl{^D^/2) = -^lnA:o(5D) + {Zaf 



92 92 / f 

0.021 511 798(1) + ^ + ^ ln(— 



(2.38) 



(2.39) 



(2.40) 



(2.41) 



2.3.5 F_5/2 states (kappa = 3) 

Now we turn our attention to F5/2 states. We have for the 4F5/2 high-energy part: 



J^H(4F5/2) = -- + (^a)' 



493 1 1 , ; 

ln(2e) 



235200 630 e 630 



4F5/2 low-energy part: 



J^l(4F5/2) = --lnA:o(4F) + [Zaf 



0.005 523 310(1) 



5F5/2 high-energy part: 

i^H(5F5/2) 

5F5/2 low-energy part: 
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{Za) 



2657 



1102500 1125 e 1125 



630 e 630 \{Zaf 



ln(2e) 



i^L(5F5/2) - --lnfco(5F) + {Zaf 



0.006 045 398(1) 



1125 e 1125 \{Za) 



(2.42) 



(2.43) 



(2.44) 



(2.45) 



2.3.6 F_7/2 states (kappa = -4) 

States with k = —4 yield the following results. We obtain for the 4F7/2 state (high-energy part): 



^i/(4F7/2) = ^ + (^«)' 



25357 1 1 
8467200 " 630 e " 630 ^' 



4F7/2 low-energy part: 



J^l(4F7/2) = --lnA:o(4F) + [Zaf 
5F7/2 high-energy part: 



0.005 180 461(1) + ^ + ^ ln(^ 



J^h(5F7/2) - ^ + {Zaf 



774121 



211680000 1125 e 1125 



ln(2e) 



5F7/2 low-energy part: 

i^L(5F7/2) = -^lnfco(5F) + + [Zaf 



0.005 662 248(1) + ^ln(^ 



(2.46) 



(2.47) 



(2.48) 



(2.49) 
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2.3.7 G_7/2 states (kappa = 4) 



G7/2 and Gg/2 states involve the most problematic angular momentum algebra of all atomic states 
considered here. We obtain, 



1 



4397 



6048000 4725 e 4725 



and for the low-energy part, 



J^l(5G7/2) = --lnfco(5G) + {Zaf 



0.001 834 827(1) 



4725 e 4725 



ln(2e) 



In 



{Zaf 



(2.50) 



(2.51) 



2.3.8 G_9/2 states (kappa = -5) 

We conclude by investigating states, where the high-energy part reads (n — 5) 



i^ff (5G9/2 

and the low-energy part reads, 



90 



+ [Zaf 



269 



283500 4725 e 4725 



= --lnfco(5G) + {Zaf 



0.001 757 471(1) 



ln(2e) 



In 



4725 e 4725 



{Zaf 



(2.52) 



(2.53) 



2.4 Summary of Results 

By adding the high- and low-energy parts given in Sees. 2.3.1 ~ 2.3.8, wc obtain the following results 
for the self energy. For all atomic states, the dependence on the scale-separating parameter e cancels 
when the low- and the high-energy parts are added, in agreement with the principal idea of the e-method 
described in Sec. 2.2. The coefficients A40, Agi and Ago can then be read off according to the definitions 
contained in Eqs. (1.1) - (1-6). 









Pi/2 state (k ~ 


1) 




n 




A40 


A40 (num) 




Aqi (num) 


Aqq (num) 


2 


1 

6 


-| lnfco(2P) 


-0.126644 


103 
108 


0.572 222 


-0.998 904 402(1) 


3 


1 

6 


- 1 lnfco(3P) 


-0.115746 


268 

405 


0.661728 


-1.148 189 956(1) 


4 


1 

6 


~| lnfco(4P) 


-0.110727 


499 
720 


0.693 056 


-1.195 688 142(1) 


5 


1 

6 


- 1 lnfco(5P) 


-0.107954 


796 
1125 


0.707 556 


-1.216 224 512(1) 



Table 2.1: Results for P1/2 states. The coefficients A40 and Aqi are rational numbers; in the 
columns entitled "(num)" , wc give the numerical values to six decimals. The Ago-coefficicnts 
in the last column constitute the main result of the investigation presented in the current 
chapter. As explained in Sec. 2.1, these coefficients involve numerical integrals which cannot 
be done analytically. 



2.5 Typical Cancellations 

As already emphasized in Sec. 2.1, severe numerical cancellations are encountered in the calculation of 
the Ago-coefficicnts. These cancellations will be illustrated here in Tables 2.8 and 2.9 for the 5G7/2 state. 
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P3/2 state (k = 


-2) 




n 




A40 (num) Aei 


Aei (num) 


^60 (num) 



2 - I lnA:o(2P) 0.123355 §§ 0.322222 -0.503 373 465(1) 

3 - I lnA:o(3P) 0.134254 i§ 0.365432 -0.597 569 388(1) 

4 - I lnA:o(4P) 0.139273 i§ 0.380556 -0.630 945 796(1) 

5 - I lnA:o(5P) 0.142046 ^ 0.387556 -0.647 013 509(1) 



Table 2.2: Analytic results for P3/2 states. 



D3/2 state (k = 2) 

n A40 A4Q (num) Aqi Aqi (num) Ago (num) 

3 - ^\nko{3D) -0.043024 0.009877 0.005 551 575(1) 

4 - I lnfco(4i:') -0.041012 ^ 0.011 111 0.005 587 985(1) 

5 - I lnfco(5L') -0.039866 ^ 0.011682 0.006 152 175(1) 

Table 2.3: Analytic results for D3/2 states. 



2.6 Observations 

It is hoped that the analytic calculations described in the previous sections will be supplemented in the 
near future by a numerical treatment of the problem, which follows the lines indicated in Ch. 1. Work 
performed recently on relevant angular functions [93] is expected to become useful in that context. In 
order to facilitate the comparisons and consistency checks which will be of relevance for the numerical 
calculations, we report here on a number of observations which should facilitate the comparison of nu- 
merical and analytic approaches to the self energy problem. They are relevant for all atomic states under 
investigation here. 

(i) Our analytic calculations for the high-energy part establish the structure of nonvanishing A4Q, Aqi 
and Ago-coefficients for the states under investigation. 



n 






^40 


D5/2 state (k = 
A40 (num) Aei 


-3) 
Aqi (num) 


Ago (num) 


3 
4 
5 


1 

30 
1 

30 
1 

30 


4 
3 

4 
3 

4 
3 


lnA:o(3i:') 
In koiiD) 


0.040310 
0.042 321 
0.043 468 


4 

405 
1 

90 

92 
7875 


0.009 877 
0.011111 
0.011682 


0.027 609 989(1) 
0.031 411 862(1) 
0.033 077 571(1) 



Table 2.4: Analytic results for D5/2 states. 
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F5/2 state (k = 


3) 




n 




Am 


A4Q (num) 


Aei 


Aei (num) 


(num) 


4 


1 

42 


-llnko{4F) 


-0.021498 


1 

630 


0.001587 


0.002 326 988(1) 


5 


1 

42 


- 1 lnA:o(5F) 


-0.020 873 


2 

1125 


0.001778 


0.002 403 159(1) 






Table 2.5: Analytic 


results for F5/2 states. 










Fy/2 state (k = 


-4) 




n 




^40 


A40 (num) 


^61 


Aei (num) 


^60 (num) 


4 


1 

56 


1 lnfco(4F) 


0.020169 


1 

630 


0.001587 


0.007 074 960(1) 


5 


1 

56 


1 lnfco(5F) 


0.020 793 


2 

1125 


0.001778 


0.008 087 021(1) 



Table 2.6: Analytic results for F7/2 states. 



(ii) For a given k quantum number, the coefficients increase in absolute magnitude with increasing 
principal quantum number n. This is the case for Bethe logarithms, for Aqi coefficients and also for 
Aqq coefficients. This behavior is consistent with the phenomenon that relativistic corrections to the self 
energy become more important as n increases. 

(iii) For a given principal quantum number n, the coefficients decrease in absolute magnitude with in- 
creasing K. As a function of the nuclear charge number Z, the plots of F{Za) against Z should become 
"more flat" as k increases for a given n [94]. 

(iv) The Aqi coefficients, which represent the dominant relativistic correction to F{Za), all have the 
same (positive) sign for states with k quantum numbers k = 2, —3, 3, —4, 4, —5. This "same-sign" pattern 
is also true of the Aqq coefficients. Because Agi is the dominant correction in the order (Za)^ due to the 
enhancement by the large logarithm, this observation explains why, in general, the qualitative dependence 
(on the nuclear charge number Z) of the self energy correction for the highly excited states is the same 
for the states with k = 2, —3, 3, —4, 4, —5 (see [94]). 

(v) Agi-coefficients are different for P1/2 and P3/2 states, but the spin dependence vanishes for higher 







G7/2 state (k = 


4) 




n 


^40 


A40 (num) Aei 


Aei (num) 


Aqq (num) 


5 


- 1 ln/co(5G) 


-0.012860 4^ 


0.000423 


0.000 814 414(1) 






G9/2 state {k = 


-5) 




n 


^40 


A40 (num) Aei 


Aei (num) 


Aqo (num) 


5 


lnfco(5G) 


0.012141 ^ 


0.000423 


0.002 412 929(1) 



Table 2.7: Analytic results for G7/2 and G9/2 states. 
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Contributions to the Low-Energy Part (5G7/2) 



160 



Aeo-contribution due to F„q 0.002 875 830 9(5). 

Aeo-contribution due to Fsy -0.001 083 109 4(5). 

Ago-contribution due to Fsn -0.008 917 782 1(5). 

Aeo-contribution due to Fg^. 0.004 920 556 0(5). 

Aeo-contribution due to Fs^ 0.004 039 332 1(5). 
[from low-energy, cf. Eq. (2.51)] 0.001 834 827(1) 



Table 2.8: As explained in Refs. [1, 2], the low-energy part of the self energy can be separated 
naturally into the nonrelativistic dipole term (which does not contribute to ^eo)- and into the 
following contributions: the nonrelativistic quadrupole part -Fnq, the relativistic corrections 
to the current F^y, and relativistic corrections to the hamiltonian Fsyi, to the bound-state 
energy F^e and to the wavefunction Fg^. Observe that the total contribution to Ago to the 
low-energy part which read 0.001 834 827(1) is almost an order-of-magnitude smaller than 
the biggest individual contribution (from -FIsh), due to mutual cancellations. 



A6o{Fh) -0.001 020 413 
Ao(-Fl) 0.001 834 827(1) 
Ago 0.000 814 415(1) 

Table 2.9: An additional numerical cancellation occurs when the finite contributions to 
Ago originating from the low-energy part [Ago(Fi), see Eq. (2.51)] and the high-energy 
part [Aqq{Fh), see Eq. (2.50)] are added. The result for the total Ago is also contained in 
Table 2.7. 



orbital angular momentum L. Consequently, Agi -coefficients are the same e.g. for -D3/2 and -D5/2 states, 
or for i^5/2 and -F7/2 states etc. This can be easily checked against numerical results. 

Preliminary results indicate satisfactory consistency between the analytic and numerical approaches to 
the self energy problem for highly excited states [95]. 
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Chapter 3 

The Two— Loop Self— Energy 



3.1 Orientation 

In the current Chapter, we investigate two-loop higher-order binding corrections to the fine structure, 
which contribute to the spin-dependent part of the Lamb shift. Our calculation focuses on the so-called 
"two-loop self energy" involving two virtual closed photon loops. For bound states, this correction has 
proven to be notoriously difficult to evaluate. The calculation of the binding corrections to the bound- 
state two-loop self energy is simplified by a separate treatment of hard and soft virtual photons. The two 
photon-energy scales are matched at the end of the calculation. We present results for the fine-structure 
difference of the two- loop self energy through the order of a* . 

3.2 Introduction to the Two— Loop Self— Energy 

As already discussed in Sees. 1.3.1 and 2.1 - 2.3, radiative corrections can be described - for atomic 
systems with low nuclear charge number - by a nonanalytic expansion in powers of the three parameters 
(i) a (the fine-structure constant), (ii) the product Za (Z is the nuclear charge number), and (iii) the 
logarithm ln[(Za)~^]. The expansion in powers of a, which is the perturbation theory parameter in 
quantum electrodynamics (QED), corresponds to the number of loops in the diagrams. The bound-state 
effects are taken into account by the expansions in the two latter parameters (see also Ch. 2). Higher-order 
terms in the expansions in powers of Za and ln[(ZQ!)~^] are referred to as the "binding corrections". One 
of the historically most problematic sets of Feynman diagrams in the treatment of the Lamb shift for 
atomic systems has been the radiative correction due to two closed virtual-photon loops shown in Fig. 3.1. 

Let us recall at this point that even the evaluation of higher-order binding corrections to the one-loop 
self energy - see also Ch. 2 -, which a 'priori should represent a less involved calculational challenge, 
has represented a problem for analytic evaluations for over three decades [1, 2, 9, 32, 70, 71, 96]. The 
energy shifts of the bound states due to the radiative corrections are conveniently expressed by expansion 
coefficients corresponding to the powers of Za and ln[(ZQ;)~^]; the naming convention is that the power of 
Za and the power of the logarithm are indicated as indices to the analytic coefficients [see also Eq. (3.1)] 
below. Because the expansion in both the one-loop and the two-loop case starts with the fourth power of 
Za, the non-vanishing coefficients carry indices A}a and Bfej for the one- and two-loop cases, respectively 
(with fc > 4 - see [49, 97] for comprehensive reviews). 

Logarithmic corrections with I > 1 can sometimes be inferred separately in a much simplified approach, 
e.g. by considering infrared divergent contributions to electron form factors. By contrast, the higher- 
order non-logarithmic coefficients represent a considerable calculational challenge. Realistically, i.e. with 
the help of current computer algebra systems [9(i], one can hope to evaluate non-logarithmic coefficients 
of sixth order in Za. Complete results for the one- loop higher-order correction Aqq for S and P states 
have only been available recently [1, 2, 32]. Calculational difficulties have by now precluded a successful 
evaluation of the corresponding coefficient B^q for the two-loop effect. Ground-work for the evaluation of 
Bqo was laid in Here, we are concerned with the evaluation of the fine-structure differences of the 
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Figure 3.1: Feynman diagrams representing the two-photon elec- 
tron self energy. The double line denotes the bound electron prop- 
agator. The arrow of time is from left to right. 



logarithmic and non-logarithmic coefficients BeL (where L = 0, 1, 2), i.e. with the ftP3/2~'^Pi/2 difference 
of these coefhcients. 

Using natural Gaussian units (ft = c = cq = l)j a-s it is customary for the current type of calculation, we 
write the two-photon self energy in the Za-expansion for P states in terms of ^-coefficients as 



AEsE = 



{Zaf 



B40 + [Zaf [Bq2 \n[Za)-^ + Bqi \n{Zc 



B, 



60 



n 



(3.1) 



where the remainder TZ is of order 0{Za)^ . Relevant Feynman diagrams arc shown in Fig. 3.1. 

Here, m denotes the electron mass (we write Eq. (3.1) in the non-recoil limit, i.e. for an infinite nuclear 
mass). The double logarithmic i?62-cocfficient is spin-independent, so that we have Afsi?62 = 0. Here, we 
evaluate the fine-structure differences 



fs^ei 



S6l(?lP3/2) - ^6l('^Pl/2) ■ 



AfeSeo = B6o(nP3/2) -S6o(nPi/2). (3.2) 

We follow the convention that AfgX = X(77,P3/2) — -''^(^Pi/2) denotes the "fine-structure part" of a 
given quantity X. For AtsiJgi and AfsBgo- we provide complete results. It is perhaps worth noting 
that the two-loop self energy for bound states has represented a considerable challenge for theoretical 
evaluations. Our investigation represents a continuation of previous work on the two-loop problem (see 
e.g. [98, 99, 100, 101]). It is probably a triviality to express that technical difficulties in the calculation 
and its description in the following Sections cannot be avoided. 

For the description of the self energy radiative effects - mediated by hard virtual photons -, we use the 
modified Dirac hamiltonian 

H^:"^ =a-[p-cFi(A)A]+pm + eFi{A)(t) + F2{A)^ {ij-E-f3(7-B), (3.3) 

2 m 

which approximately describes an electron subject to an external scalar potential = (t>(r) and an 
external vector potential A = A{r). This modified hamiltonian is still local in coordinate space. The 
Dirac matrices in (3.3) are to be understood in the standard (Dirac) representation [.S-l] (in the sequel, 
we will also use the non-covariant notation /3 = 7" and a* = 7^7'). 
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The argument A of the electron form factors Fi and F2 in Eq. (3.3) is to be interpreted as a Laplacian 
operator acting on all quantities to the right (but not on the wave function of the bound electron in 
evaluating H^^^lip)). In momentum space, the action of the hamiltonian iJ^"' is described by the convo- 



lution 



H^'^p {p') = J d'^p/(27r)'^ H^'{p' — p) ip{p). The form factors - in momentum space - assume 

arguments according to the replacement A = — (p' — p)^. In Eq. (3.3), radiative corrections are 

taken into account in the sense of an effective theory via the inclusion of the on-shell form factors Fi 
and F2. Although the bound electron is not an on-shcU particle, the modified hamiltonian (3.3) can still 
approximately account for significant radiative systems with low nuclear charge number Z. Of course, 
the hamiltonian (3.3) cannot offer a complete description of the bound electron. Recoil effects cannot be 
described by a one-particle equation in principle, and vacuum-polarization effects are not contained in 
Eq. (3.3). However, the effective description of self energy radiative corrections mediated by hard virtual 
photons given by Eq. (3.3) will turn out to be useful in the context of the current investigation of the 
two-loop self energy. 

Both of the form factors Fi and F2 entering in Eq. (3.3) are infrared divergent, but this divergence is 
cut off in a natural way at the atomic binding energy scale {Za)'^ m. The fact that on-shell form factors 
can describe radiative corrections to the fine structure - mediated by high-energy virtual photons - has 
been demonstrated explicitly in [89]. The modified Dirac hamiltonian (3.3) and the associated modified 
Dirac equation have been introduced - in the one-loop approximation - in Ch. 7 of [34] [see for example 
Eqs. (7-77) and (7-103) ibid]. The low-energy part of the calculation is carried out using nonrelativistic 
approximations in the spirit of the simplified treatment introduced in the previous one- and two-loop 
calculations [1, 2, 32, 98, 102]. This approach was inspired, in part, by various attempts to formulate 
simplified low-energy ( "nonrelativistic" ) approximations to quantum electrodynamics ( "NRQED" ) , see 
e.g. [103, 104]. Both the high-energy and the low-energy contributions are matched at the separation scale 
e whose role in the calculation is illustrated by the mathematical model example discussed in Sec. 2.2. 

In a two-loop calculation, cither of the two virtual photons may have a high or low energy as compared to 
the separation scale e. A priori, this necessitates [!)(S] a separation of the calculation into three different 
contributions: (i) both photon energies large, (ii) one photon with a large and one with a small energy, and 
(iii) both photons with small energies. For the particular problem at hand (the fine-structure differences 
of Bqi and Bqq), we are in the fortunate position that effects caused by hard virtual photons (i) are 
described by the modified Dirac hamiltonian (3.3), whereas the low-energy part discussed in Sec. 3.5 
below comprises both remaining contributions (ii) and (iii). 

We will continue by analyzing the two- loop form factors which enter into Eq. (3.3); these will be treated 
in Sec. 3.3. The calculation will be split into two parts: the high-energy part discussed in Sec. 3.4 and 
the low-energy part, which is treated along ideas introduced in [104] in Sec. 3.5. 



3.3 Two-loop Form Factors 

In order to analyze the modified Dirac hamiltonian (3.3) through two-loop order, we first have to in- 
vestigate certain expansion coefficients of the electronic Fi and F2 form factors which are thoroughly 
discussed in the seminal papers [105, lOG]. For the momentum transfer which is the argument of the 
two functions Fi = Fi{q^) and F2 = F2{q^), wc use the convention — q^q^ — {q^Y — ■ The variable t 
in [105, lOG] is given ast = q^ . When we evaluate radiative corrections to the binding Coulomb field which 
is mediated by space-like virtual photons, we have q^ = —q^ because for = 0. Wc use the conventions 
(see Eq. (1.2) in [10-5]): 

F^{t)^i+Y.{^yFr{t), ^^.(o=EO"^i'"^w- (3-4) 

n— 1 n— 1 

One and two-loop effects are denoted by upper indices 2 and 4, respectively. This notation is motivated 
by the observation that two-loop effects are of forth order in the quantum electrodynamic interaction 
Lagrangian —e^^^A^ip (in the Furry picture, which is used for the description of bound states, the 
Coulomb interaction is taken out of the interaction Lagrangian). 

There are two different points of view regarding the choice of diagrams to be included in the two-loop 
form factors, depending on whether the self energy vacuum polarization diagram Fig. 3.2 is included in 
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the calculation or not. Wc will discuss both cases and give results with and without the diagram shown 
in Fig. 3.2 taken into account. 




{V) 



Figure 3.2: Combined self energy vacuum-polarization diagram (denoted "y 
in the text). 



First, we discuss results obtained for Fi including the combined self energy vacuum polarization diagram. 
In this case, the known results for the slope i^{(0) and for ^2(0), through two-loop order, read as follows. 
From Eq. (1.11) of [105], we have: 

m^F[{0) = - 







1 












^ 8 



4819 49 3 

_ ^(2) + 3 C(2) In 2 - - C(3) 



where the forth-order coefficient has the numerical value 

rn^Fi'^'^^ (0) = 0.469 941 487 460 . 
According to Eq. (1.7) in [105], the value of ^2(0), through two-loop order, reads 



^^2(0) 



1 a 

2 n 



197 1 



+ -C(2)-3C(2) In2 + -C(3) 



144 2 

where the two-loop coefficient has the numerical value 

F^'^^ (0) = -0.328 478 965 579 . 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



We now turn to the discussion of the slope F2'^\o). In view of Eq. (1.20) of [1U5] (see also [107]), we have 
(up to two-loop order) 



A 



ln-B(t).Ff (t)+.Fr(t) 



(4)^ 



where the coefficients J-^ are by definition infrared safe and 



^f(0) = -, B{t) = -- ^ 
2 3 



1 



t 



20 m2 



(3.9) 



(3.10) 



Equations (3.9) and (3.10) uniquely determine the infrared divergent contribution to i^2*'''''(0). An analytic 
expressions for JF^^-'(t), t spacelike, has recently been obtained [108, 109, 110, 111] in terms of harmonic 
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polylogarithms [112, 113]. As a byproduct, an analytic expression for the slope J^2^''^(0) was found. The 
result reads 



™^^.'«(0) = i|l + |«2)-|c(2)l„2 + gc(3). (Ml) 

A numerical result for J-'2^\o), complementing the above analytic expression, can easily be derived in 
combining Eq. (1.20), Eq. (1.30), and Eq. (3.2) in [1U5], as will be explained in the sequel. The dispersion 
relation (1.30) in [105] reads, 

oo 

4m2 , , 1 t f dt' t' - All? , ,^ , ^ 

where P denotes the Cauchy principal value. Equation (3.12) applies also if we single out the two-loop 

- 2 



effect and differentiate at zero momentum transfer, and we obtain for the slope -F'2''^^(0) the relation 



F',''' (0) = \ iO) + ^P J dt' Im^^^(^) (0 = J (Q) + T , (3.13) 



where F2^^(0) is given in Eq. (3.8). The second term on the right-hand side, denoted by T, can be 
evaluated using the result for ImF2 {^) presented in Eq. (3.2) in [ ]; it reads 

1 ^ 

- / da: ~ Tm ^ (x) = --ln(—\ + 0.030 740 507 833(1) . (3.14) 
7 a; (1-1- x) 6 \rn J 



Here, the last error is due to numerical integration, and use is made of the natural variable [105] 



(3.15) 

1 + Vl-4m2/i 

In combining the resuh of Eq. (3.8) with Eqs. (3.13) and (3.14), the result tti^ J^!}'^\0) = 
—0.051 379 233 561(1) is obtained which is in agreement with (3.11). 

Now we will provide results for the form factors obtained excluding the self energy vacuum-polarization 
graph V shown in Fig. 3.2. These results refer to the pure two-photon self energy diagrams shown in 
Fig. 3.1. The two- loop self energy diagrams independently form a gauge-invariant set. They represent 
a historically problematic correction, and are the main subject of our investigation. The combined self 
energy vacuum-polarization diagram, according to Eqs. (1.9) and (1.10) in [106] - taking into account 
the subtracted dispersion relation (1.30) of [105] - leads to the following corrections: 

^^{^'^'^(0) = -i^ + ^C(2)- 0.031588972474, 



F2^'*)'^(0) = ^ - 2C(2) = 0.015687421859, 

oD 



F^^'^^'^iO) = —- -C(2) = 0.000 497 506 323. (3.16) 
216 8 
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For the pure self energy graphs, which we would like to denote by the symbol S, we therefore obtain the 
following results, 



^ - ^ C(2) + 3 C(2) In 2 - J C(3) = 0.438 352 514 986 , (3.17) 



(4),S 



(0) = + ^C(2)- 3 C(2) In2 + ^C(3) = -0.344166 387438, 



(3.18) 



m2F2'(^''^(0) 



1, /A\ 151 61 23 , „ 23 

-- n — h — C 2 C 2) n2H C(3) 

6 \m 240 40 ^^ 10^^^ 40^^'^ 



V TO 



0.051876 739 885: 



b V TO 



/(4).S 



(0) 



(3.19) 



where the latter equality defines ^^^^'^■'''^(0) in analogy with Eqs. (3.9) and (3.11). 



3.4 High— Energy Part 

Based on the modified Dirac hamiltonian (3.3), corrections to the energy of the bound Dirac particle can 
be inferred. We will refer to the energy corrections attributable to the Fi and F2 form factors as Ei and 
E2, respectively. For Ei, we have 

i?i = ([^^i(-g')-l] e0)f^ , (3.20) 

where the index fs refers to the fine-structure terms, i.e. to the result obtained by subtracting the value 
of the matrix element for a nP3/2 state from the value of the same matrix element evaluated on a nPi/2 
state. A matrix element {A)f^ of a given operator A is evaluated as 



A 



Note that in evaluating these matrix elements, (the hermitian conjugate of the Dirac wave function 
ip) should be used (not the Dirac adjoint ip = ijj'^j'^). The Dirac wave functions ip are expanded in powers 
of (Za) up to the order relevant for the current investigation of two-loop effects. This expansion avoids 
potential problems associated with the logarithmic divergence of the Dirac wave function at the origin. 

For El, up to the order of (Za)^, we have 



Ei^AirZaF^^^^Q) {6^^\r) 



(3.21) 



For P states, the nonrelativistic (Schrodinger) wave function - the leading term in the Za-expansion of 
the Dirac wave function - vanishes at r = 0, but the first relativistic correction gives a finite contribution, 
resulting in 



{Za)^m^ 



4 77,5 



This leads - again up to the order of (Za)^ - to the following result for Ei, 



El 



a^^2 {Zaf 



(3.22) 



(3.23) 



Observe that the derivative of the Fi form factor has a physical dimension of 1/to^ in natural units, 
giving the correct physical dimension for Ei. The correction due to F2 in (3.3) reads. 



E2^{F2i-q^) — n-E)f,. 



(3.24) 
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A particle in an external binding Coulomb field feels an electric field E = i (Ze) q/q^ - in momentum space 
- 01 E = — (Ze) r/(47rr^) in coordinate space. Vacuum polarization corrections to E = — (Ze) r/(47rr'^) 
lead to higher-order effects. The correction E2 splits up in a natural way into two contributions E2a and 
£^26 which are associated with ^2(0) and the slope ^2(0)1 respectively. £^2a reads 



(3.25) 



The evaluation of the matrix element leads to 



. 7 ■ f" 



{Za) 



487 
360 



5_ 

An 



23 



10n2 



[Zaf 



(3.26) 



For the purpose of the current investigation, the (Za)^-component of i?2a is selected only: 



E2a — 



{Zaf 



^2 \j2Q^^ 



23 \ 
20n2 / 



m . 



(3.27) 



The matrix element £^26 can be expressed as 



E2b — 



AirZa 
2m 



Ff'(O) (7-g)fs. 



(3.28) 



A transformation into coordinate space leads to 



(7 • q)is = i 



- 1 



{Za)^m^ 



x=OJs 



(3.29) 



As a function of the principal quantum number n, the result for £26 reads: 



E2b — 



a\2 [Zoif 



1 



(3.30) 



This result involves the infrared divergent slope of the F2 form factor [see Eqs. (3.11) and (3.19)]. We 
are thus faced with the problem of matching the infrared divergence of the slope of the F2 form factor, 
expressed in terms of the fictitious photon mass A, with the usual (energy matching parameter) e intro- 
duced originally in [-Si]. This can be done in two ways: (i) by matching the infrared divergence of the rate 
of soft bremsstrahlung, calculated with a fictitious photon mass A, to a result of the same calculation, 
carried out with an explicit infrared cut-off e for the photon energy. This way of calculation is described 
on pp. 361-362 of [34]. It leads to the result 



2e 



(3.31) 



The matching procedure (ii) consists in a comparison of the result of the application of the formalism 
considered above, and its application to the high-energy part of the ground state Lamb shift, which is in 
leading order given by the infrared divergence of the F\ form factor, and the result obtained by direct 
calculation of this high-energy part in a non-covariant formalism with an explicit energy cut-off e, as it 
has been carried out in [32] . This second matching procedure leads to the following result - in agreement 
with (3.31) 

"1 , m 11 , , 

In — + — . (3.32) 



In- 



A 



So, we are led to the replacement 



InA 

m 



2e 



m 

In 

2e 



11 
24 

5 
6 



(3.33) 



A comparison with the results in Eqs. (3.11), (3.19), and (3.36) reveals that the logarithmic divergence 
for the fine-structure difference is given by a term 



I — 1 m 
"3^ Ye' 



(3.34) 
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so that we may anticipate at this stage the resuh for AfsBgi, 

n2 - 1 
Afs-Bei — „ o ■ 

Based on (3.30) and (3.33), we can express i?2h in terms of e and J^2'^''(0), 



E2b = 



a\2 {Zaf 



1 



1 7Tl 

— In 



2e 



18 



+ 2.F2*''^(0)to 



71^ - 1 



(3.35) 



(3.36) 



There is a third correction due to the effect of two one-loop corrections on the electron vertices. Because 
we are only interested in the fine structure, we isolate the terms which are proportional to the spin-orbit 
coupling, and we obtain 



where 



Es = ( [2F2(0)] i/fc 



1 



E- H 



Za 



cr ■ L . 



(3.37) 



(3.38) 



and \/{E — H)' is the nonrelativistic, spin-independent reduced Schrodinger-Coulomb Green func- 
tion [114, 115]. The only spin-dependence in (3.37) occurs in the coupling cr ■ L, and it can be taken 
into account by an overall factor, 



(3.39) 



We are therefore led to consider the "spin-independent version" of the matrix element which occurs in 
Eq. (3.37) and obtain the following result, 



Za 



1 



4^2j,3 \ E — H J Avii^r^ 



Za 



nP 



227 



1 



1 



8640 96n SOn^ 



[ZaY m 



(3.40) 



The spin-dependence can be easily restored by considering Eq. (3.39). The index "nP" in Eq. (3.40) 
means that the matrix element is evaluated with the nonrelativistic, spin-independent (Schrodinger) 
wave function. Alternatively, on may evaluate with cither the nPi/2 or the nP3/2 Dirac wave function 
and expand up to the leading order in (Za). 

The evaluation of (3.40) can proceed e.g. by solving the differential equation which defines the correction 
to the wave function induced by Hfs , and subsequent direct evaluation of the resulting standard integrals 
using computer algebra [90]. The final result for E3 reads. 



E. = 



[Zaf 



227 



1 



(3.41) 



2880 32?i 80n2_ 

This concludes the discussion of the high-energy part. The final result for the high-energy part is 

Eh = Ei+ E2a + E2b + E3 , (3.42) 
where Si, E2a, £^26, E3 are given in Eqs. (3.23), (3.27), (3.36), (3.41), respectively. 



3.5 Low— Energy Part 

The low-energy part consists essentially of two contributions. Both effects, denoted here E4 and E5, can 
be obtained by a suitable variation of the low-energy part of the one-loop self energy, by considering the 
spin-dependent effects introduced by a further one- loop electron anomalous magnetic moment interaction. 
The first of the two terms, E4, is caused by spin-dependent higher-order effects in the one- loop self energy, 
which receive additional corrections due to the anomalous magnetic moment of the electron. The second 
term, E^, is due to an anomalous magnetic moment correction to the electron transition current, which 
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can also be seen as a correction to the radiation field of the electron due to its anomalous magnetic 
moment. 

The leading-order low-energy part (see [32]) reads 

2a 



= -- 



duj L 



1 



P 



H-{E~uj) 



P 



(3.43) 



In order to isolate the fine-structure effects, we should now consider corrections to the wave function, to 
the current, to the hamiltonian and to the energy of the bound state due to the spin-dependent relativistic 
(spin-orbit) hamiltonian 

e . _ ol(Zol) — i7 • r 



U 



F2{Q) — i-i-E 
2m 



(3.44) 



47rm 

The above hamiltonian 7i is the last term in the modified Dirac hamiltonian [right-hand side of Eq. (3.3)], 
approximated for a particle bound in a Coulomb field with the F2 form factor evaluated at zero momen- 
tum. The electric field E in (3.44) corresponds to the binding Coulomb interaction. The hamiltonian 
(3.44) describes the modification of the spin-orbit interaction due to the anomalous magnetic moment of 
the electron. 

The nonrelativistic limit of Ti is the spin-orbit coupling His given in Eq. (3.38), multiplied by a factor 
2F2 (0) = q/tt (the additional factor 2 finds an explanation in [116]). The resulting hamiltonian 



a a Za 
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(3.45) 



takes into account magnetic vertex corrections in the framework of an effective theory. Denoting the 
variation of the expression (3.43) mediated by iJoff with the symbol ^off ~ in the spirit of the notation 
introduced in [98] -, we obtain the contribution 
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(3.46) 



Following the notation introduced in [1, 2], the contribution £"40 is the sum of the fine-structure effects 
created by the wave- function-correction -F^^, the first relativistic correction to the energy Fse, and the 
correction due to the relativistic hamiltonian Fsh-, each multiplied by a factor a/'K. The final result for 
Ei is 

- {Zaf — {AisFs^ + AfsFsE + Af^FsH) ■ (3.47) 

There is a further correction to the nonrelativistic effective coupling to the radiation field due to the 
"anomalous spin-orbit hamiltonian" (3.44). The correction, in the nonrelativistic limit, can be derived by 
considering a Foldy Wouthuysen transformation which by definition diagonalizcs the hamiltonian (3.44) 
in spinor space and also leads to corrections to the current according to 



a 



Ua'U-\ C/==cxpf-i^ 
\ 2m, 



(3.48) 



Here, f3 and a' are standard Dirac matrices [34], i is a spatial index, and H is given in (3.44). The 
calculation is carried out along ideas introduced in [1] and leads to the result 
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(3.49) 



as a relativistic correction to the electron current which is simply a* in the relativistic formalism and 
p^/m in the leading nonrelativistic approximation. Again, following the notation introduced in [1, 2], the 
resulting additional contribution is 



E. 



4b 



- (Zaf^AisFsy 



(3.50) 



The sum of (3.47) and (3.50) is just the (Za)^ -component of the fine-structure difference of the one-loop 
self energy from [1, 2], multiplied by an additional factor a/ir. It can also be written as 



E. 



E. 
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(3.51) 
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where Afs^4(n) could be interpreted as a relativistic generalization of a Bethe logarithm, which is n- 
dependent. However, a significant numerical fraction of the n-dependence can be eliminated if the factor 
{n'^ — l)/n? is taken out of the final result. The evaluation of Afs^4(n) has recently been performed in [117] 
with improved numerical methods (see e.g. [23]), and the following results have been obtained: 



Afe£4(2) 


= 0.512 559 768(1), 


Afs4(3) 


= 0.511 978 815(1), 


Afs£4(4) 


= 0.516 095 539(1), 


Afe£4(5) 


= 0.519 976 941(1), 



(3.52) 

where the uncertainty is due to numerical integration. 

There is, as stated above, a further correction due to the explicit modification of the transition current 
due to the anomalous magnetic moment; it can be obtained through the replacement 

a'^a' + F2iO)^ (3.53) 

Zm 

and must be considered in addition to the correction (3.48). A careful consideration of the nonrelativistic 
limit of this correction to the current, including retardation effects, leads to the result 

27r 2mr-' 

Consequently, we find that the correction is effectively ^2(0) times the retardation corrections to the 
transition current Fgy found in [1, 2]. We obtain 

E.= n\zarII^^. (3.55) 
In analogy with E^, this correction can favorably be rewritten as 



^5 = 

On the basis of [1, 2, 117], we obtain 



- 1 



Afe4(n) 



(3.56) 



Afe£5(2) = 


-0.173 344 868(1) , 


Afe£5(3) = 


-0.164 776 514(1), 


Afe4(4) = 


-0.162 263 216(1), 


Afe4(5) = 


-0.161 165 602(1). 



(3.57) 

The final result for the low-energy part is 

Ei^^Ei + E^, (3.58) 
with Ei and E^, being given in Eqs. (3.51) and (3.56), respectively. 

We can now understand why it was possible to join the two contributions with "mixed" and "low-and- 
low" energy virtual photons (ii) and (iii), which were discussed in Sec. 3.2, into a joint "low-energy part". 
The reason is simple: The effective hamiltonian (3.45) has no infrared divergence, because it involves the 
low-energy limit of the magnetic form factor _F2, which is infrared safe in one-loop order according to 
Eq. (3.9). Because the main contribution to the quantity i^2(0) is caused by hard virtual photons, it is 
also justified to say that the contribution of "low-and-low" energy virtual photons vanishes at the order 
of interest for the current calculation (fine-structure difference). In higher- loop order, the further infrared 
divergence acquired by F2 would lead to an infrared divergence in the effective hamiltonian constructed in 
analogy with Eq. (3.45); this infrared divergence would have to be attributed to a "mixed" contribution. 
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3.6 Results for the Two— Loop Corrections 



Wc have obtained analytic results for higher-order correction to the two-loop self energy of P states in 
hydrogen-like systems. In our calculation, we have analyzed the electron form factors through two-loop 
order in Sec. 3.3, and we have split the calculation into a high-energy part with two hard virtual photons 
discussed in Sec. 3.4, and a low-energy part with at least one soft virtual photon analyzed in Sec. 3.5. 
The final result for the contribution to the fine-structure energy difference is obtained by adding the 
high-energy contributions Ei - E3 given in Eqs. (3.23), (3.27), (3.36), (3.41), and the low-energy effects 
E4 and E5 from Eqs. (3.51) and (3.56). The dependence on e cancels out in the final result which is the 
sum of the high-energy part En given in Eq. (3.42) and the low-energy part i?L defined in Eq. (3.58). This 
is also evident when considering explicitly the Eqs. (3.36) and (3.51). The final results for the analytic 
coefficients of order a^{Za)^ read 

AfeSei = . (3.59) 

[see also Eq. (3.35)] and 



f: 



'('^'^(0) + 2.Ff ^'^(0)) + AtMn) + Af,4(n) 



(3.60) 



where explicit numerical results for F^^^'^'^ {0), F2^^'^{0) and ^2^^'^ (0) can be found in Eqs. (3.17), (3.18) 
and (3.19), respectively. This result refers to the pure self energy diagrams in Fig. 3.1. The result reads 
numerically for the principal quantum numbers n = 2-5, 



Afsi?6o(2) — 


-0.361 196 470(1), 


(3.61) 


Afsi?6o(3) = 


-0.411 156 068(1), 


(3.62) 


Af3B6o(4) = 


-0.419 926 624(1), 


(3.63) 


AfsB6o(5) = 


-0.419 832 876(1). 


(3.64) 



If it is desired to add in the combined self energy vacuum-polarization diagram from Fig. 3.2, then 
the form-factor results from Eqs. (3.6), (3.8) and (3.11) instead of the pure self energy results given in 
Eqs. (3.17), (3.18) and (3.19) have to be used in evaluating (3.60). When including the combined self 
energy vacuum-polarization diagram from Fig. 3.2, there is no further low-energy contribution, so that 
the alternative set of numerical values for the form factors from Eqs. (3.6), (3.8) and (3.11) fully takes 
into account the additional effect of the diagram in Fig. 3.1 on the fine-structure in the order of {Za)^. 
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Chapter 4 



Spinless Particles in Bound— State 
Quantum Electrodynamics 

4.1 Orientation 

We describe in this Chapter a simpHfied derivation of the relativistic corrections of order a'^ for a bound 
system consisting of two spinless particles. We devote special attention to pionium, the bound system of 
two oppositely charged pions. The leading quantum electrodynamic (QED) correction to the energy levels 
is of the order of and due to electronic vacuum polarization. Wc analyze further corrections due to the 
self energy of the pions, and due to recoil effects, and we give a complete result for the scalar-QED leading 
logarithmic corrections which are due to virtual loops involving only the scalar constituent particles (the 
pions); these corrections arc of order In a for S states. 

4.2 Introduction to Spinless QED and Pionium 

Exotic bound systems like pionium [118, 119] (the bound system of two oppositely charged pions) offer 
interesting possibilities for studies of fundamental properties of quantum mechanical bound states: the 
interplay between strong-interaction corrections and quantum electrodynamic corrections is of prime 
interest, and the small length scales characteristic of the heavy particles make it possible to explore effects 
of the virtual excitations of the quantum fields in previously unexplored kinematical regimes [3, 4, G]. 
We do not wish to hide the fact that any potential high-precision experiments in this area are faced with 
various experimental difficulties. Our calculations address QED corrections to the spectrum of bound 
systems whose constituent particles are spinless; relativistic corrections to the decay lifetime of pionium 
have recently been discussed in [120] in the context of the DIRAC experiment at CERN. 

Here, we report on results regarding the spectrum of a bound system consisting of two spinless particles. 
We apply the simplified calculational scheme employed in [4i i] for the relativistic and recoil corrections to 
a bound system of two "non-Dirac" particles to the case of two interacting spinless particles (see Sec. 4.3). 
We then recall known results on leading-order vacuum polarization corrections in Sec. 4.4 and clarify the 
relative order-of-magnitude of the one- and two-loop electronic vacuum polarization, the relativistic and 
recoil corrections and the self energies in pionium (also in Sec. 4.4). We then provide an estimate for the 
self energy in Sec. 4.5, and we analyze the leading recoil correction of order (the Salpeter correction) 
which leads us to complete results for the scalar-QED logarithmic corrections of order In a. 

4.3 Breit hamiltonian for Spinless Particles 

We start from the Lagrangian for a charged spinless field coupled to the electromagnetic field [see equa- 
tions (6-50) - (6-51b) of [34]], 

£(x) = [{d^ - ieA^) rix)] (a^ + ie^^) 0(x) - m^ri^) <P{.x) - ^ F^,{x) Pf^^x) , (4.1) 
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where the field strength tensor Ff^^, reads F^^{x) = d^A^{x) — d,yAfj_{x). We use natural Gaussian units 
with h = c ^ = 1. The transition current for a free spinless particle {A^^ ~ 0) can be inferred from 
(4.1); it reads in momentum space 

fip\p)^r{p')ip'^+p^)cl,{p). (4.2) 

This current now has to be expressed in terms of nonrelativistic wave functions. Specifically, the j"- 
component has to reproduce the normalization of the nonrelativistic (Schrodinger) wave function. By 
contrast, according to Eq. (4.2) the zero-component of the current reads 2m(j)*4) in the nonrelativistic 
limit p"^ m ,p'^ m. The nonrelativistic wave functions arc normalized according to 



d^x4>s{x)(l)s{x) = 1. (4.3) 

It is therefore evident that we cannot simply associate the relativistic wave function 4> with (jjg; rather, 
we should define according to Eqs. (13) - (14) of [40] 

^iP)-^^^ (4.4) 



where p^ = \/p'^ + trfi w m is the energy of the free nonrelativistic particle (in deriving low-energy 
effective interactions, one always expands about /ree-particle amplitudes; all interactions are treated as 
perturbations; note the analogy to nonrelativistic QED - NRQED - for spinor particles [104]). The 
Klein-Gordon current, in the presence of external fields, reads in contrast to (4.2) 

f{p',p) = 4>*{p') ip'^+p'' -2eAn m . (4.5) 

The zero-component of this current can be interpreted as a charge density, which is not necessarily 
positive definite. Questions related to the normalization of the Klein-Gordon wave functions in this case 
are discussed in detail in [121, 122, 123, 124]. 

In terms of the Schrodinger wave function, the current is given as 

fip'.p) - 0s(p')0s(p), (4.6) 



/(P',P) = 0g(p')^^0s(p), (4.7) 

where m is the mass of the particle. The atomic momenta p^ and p" in Eq. (4.7) are of order Za. 
As shown below, interactions involving the spatial components -p of the transition current give rise to 
relativistic contributions of order {ZaY to the spectrum. This is exactly the order of magnitude that 
is the subject of the current investigation. Therefore, although Eq. (4.7) is only valid up to corrections 
of relative order {Zaf^, these can be neglected because the further corrections contribute to the energy 
levels at the order of {Zoif . Specifically, we can expect corrections proportional to {p^ p^) to the current 
when a systematic expansion of the nonrelativistic current is performed; these terms are analogous to 
those obtained for relativistic corrections to the current of spinor particles which can be obtained via a 
Foldy-Wouthuysen transformation [1, 125]. 

In the following, the index S on the wave function will be dropped, and the nonrelativistic amplitudes 
describing the two interacting particles (with electric charges ei and 62) will be denoted as 0i and 02, 
respectively. Following [40], the Brcit hamiltonian U{pi,P2,q) in momentum space is related to the 
invariant scattering amplitude M and to the photon propagator Df^i,{q) in the following way [see also 
equation (83,8) in [126]]: 

M = eie2jf(p'i,Pi)i^^.(q)j2"(P2,P2) 



^l{p'i)4>*2{p'2: 



ei 62 



uiPi,P2,q) 



hiPl) (l)2iP2) 



(4.8) 
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where q = P2 P2 = ^{Pi ^ P'l)- We employ a Coulomb-gauge photon propagator, 



Dooiq) = , Dij{q) = 





1 




q . 




q . 



(4.9) 



where we can neglect the energy of the virtual photon in Dij for the derivation of ncxt-to-leading order 
relativistic corrections. The invariant scattering amplitude M then reads 



M 
ei 62 



1 



1 (Pl ) 02 (P2 ) -J (Pl ) 02 (P2 ) 



+0Kp'i)0;(P2) 



P\ + P'l P2 + P2 J_ 

2toi 2m2 



MPi)MP2)- (4.10) 



We therefore identify 



u{Pi,P2,q) = - 



6162 (2pl-90(V2+9^) 



4mim2 



9' 



6162 / (2pi - q) ■ {2p2 + q) 
4mim2 \ 



(2pi-q-q2)(2p2-q + g2) 



61 62 

mi 7712 



Pi • P2 (Pi • g) (P2 • q) 



We now transform to the center-of-mass frame in which Pi 
U{pi,P2,q) becomes even simpler, 



-P2 



(4.11) 

p, so that the expression for 



U{p, p, q) 



61 62 

7771 fn2 



P 



ip-q? 



The formula (83,13) of [12(i] can now be employed in evaluating the Fourier transform. 



d^g 
(27r)3 



exp {iq- r] 



47r(a ■q){h- q) 



1 

2^ 



a ■ b 



[a ■ r){b ■ r) 



(4.12) 



(4.13) 



The Breit hamiltonian, which we would like to denote by i?Bj is obtained by adding to the Fourier 
transform of (4.12) the relativistic correction to the kinetic energy. Denoting with p = —id/dx the 
momentum operator in the coordinate-space representation, we obtain 



H-air.p) = - 



P 

8 ml 



P 

87773 



61 62 p 



61 62 r ■ {r ■ p) p 



Snr 7771 TTT'2 Sttv^ 



777 1 7)7,2 



(4.14) 



In the order of (Za)^, there is no contribution due to virtual annihilation for spinless particles; corrections 
of this type would enter only for positronium and dimuonium [G] because they are caused by the spin- 
dependent part of the transition current [see Eqs. (83,20) and (82,22) of [12(i]], which is absent for spinless 
particles. For S states, virtual annihilation is altogether prohibited by angular momentum conservation. 

The matrix elements of the Breit hamiltonian (4.14) for spinless particles can be evaluated on nonrela- 
tivistic bound states via computational techniques outlined in Sec. A3 of Ch. 1 of [11 (>]. For 7771 = 7772 = 777 
and 61 62 = —AnZa, we obtain 



[Zaf 777 {ZaY 



2773 



2/ + 1 



1 e 



11 

32^ 



(4.15) 
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as the Brcit energy for the energy levels of the bound system of two spinless particles, including relativistic 
corrections of order (Za)^. Here, we keep Z as a parameter which denotes the nuclear charge number in 
a bound system. Of course, for two particles each carrying an elementary charge, Z has to be set to unity. 
The fine-structure constant is denoted by a. The result (4.15) agrees with previous calculations [40, 127, 
128, 129, 130, 131], notably with Eq. (38) of [129]. 

It is instructive to compare the result (4.15) with the known result for a single-particle system of mass 
to/2 satisfying the Klein-Gordon equation, bound to a nucleus with charge Ze. According to Eq. (2-86) 
of [34], we obtain the "Klein-Gordon energy" (KG) 



Ekg 



(Za)^ m {Za)'^m 



4n2 



2n3 



1 



21 + 1 



8^ 



(4.16) 



The two results (4.15) and (4.16) are manifestly different in the order of (Za)^. 

From (4.14) we conclude that the zittcrbcwegung term is absent for spinless particles. However, this 
statement is in need of further explanation because a considerable variety of physical interpretations 
exists in the literature with regard to the zitterbewegung term. We briefly expand: The Dirac a-matrices 
fulfill a = i[H-[:i,x] (i/o is the Dirac hamiltonian) as the relativistic generalization of the velocity operator. 
By contrast, in the nonrelativistic formalism, we have the analogous relation p/m = i[Hs-, x] where Hg is 
the Schrodinger hamiltonian. Since the a-matrices have eigenvalues ±1, the magnitude of the velocity 
of the electron - at face value - is equal to the velocity of light at any given instant. On p. 106 of [132], 
it is argued that "the explanation for this fact is that the electron carries out a fast irregular motion 
( "zitterbewegung" ) - which is responsible for the spin - whereas the mean velocity is given by the 
momentum p/m" . Note that the introduction of the Dirac matrix formalism is necessitated by the need 
to describe the internal degrees of freedom of the particle - the spin. On p. 71 of [34], it is shown that the 
zitterbewegung term can be traced back to the positional fluctuations {Sr"^) ~ l/m^ of the electron, and a 
connection is drawn to the Darwin term which results naturally in the context of the Foldy-Wouthuysen 
transformed Dirac hamiltonian. On pp. 117-118 of [ l ■ '>'■>] and p. 62 of [.34], it is argued that the momentum 
p of a Dirac wave packet can be associated in a natural way with the group velocity, but that in addition to 
the group velocity term, there exist highly oscillatory terms which represent the zitterbewegung. Similarly, 
on pp. 139-140 of [133], it is shown that the zitterbewcgimg term can also be interpreted as arising from 
the interaction of the atomic electron with virtual electron-positron pairs created in the Coulomb field of 
the nucleus. This virtual electron-positron pair-creation is subject to the uncertainty principle and can 
occur only for time intervals of the order of At ^ h/ {2m(?) (where we temporarily restore the factor h). 
At the time the original atomic electron fills up the vacated negative-energy state (the bound-electron 
wave- function has negative-energy components), the escalated electron (which forms part of the virtual 
pair) is at most a distance cAt ~ 1/m away from the original electron. This distance is precisely of the 
order of magnitude of the fluctuations of the electron coordinate and consistent with the discussion on 
p. 71 of [34]. All these interpretations elucidate different aspects of the same problem. 

In the context of the Breit hamiltonian, we would like to adhere to the definition that the zitterbewegung 
term is the term of order [Za)'^ in the Breit hamiltonian generated by a contribution which is manifestly 
proportional of 5{r) in coordinate space (or a constant in momentum space). Such a term is absent in 
the result (4.14). For spin- 1/2 particles, such a term is generated by the multiplication of the photon 
propagator (proportional to with the zero-component of the transition current which is given for a 

spin- 1/2 particle as [see equation (4) of [40]] 



-/ 0- * 

u J u ~ w 



1 



8m2 



icr ■ p X p 
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(4.17) 



Here, u is the bispinor amplitude for the bound particle, and w is the bound-state Schrodinger wave 
function related by 
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(4.18) 



according to equation (3) of [40]. One might wonder why a term proportional to (5;o, apparently generated 
by a (5-function in coordinate space, prevails in the Breit energy (4.15). This term arises naturally when 
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evaluating a matrix element of the structure {4>s\{f ■ {f ■ p) p)/'>'^\4's) (last term of equation (4.14)) on 
the nonrelativistic wave function 0g and should not be associated with the zitterbewegung. 

Vacuum polarization corrections and the self energy, as well as corrections due to the strong interaction, 
are not included in (4.15). These corrections will be discussed in the two following sections. 



4.4 Vacuum Polarization Effects 

As pointed out by various authors (e.g. [5, 7, 134, 135, 136, 137, 138]), the electronic vacuum polarization 
enters already at the order of [more precisely, a (Zot)^\ in bound systems with spinless particles, 
because the spinless particles are much heavier than the electron, which means that the Bohr radius of 
the bound system is roughly of the same order of magnitude as the Compton wavelength of the electron. 
The Compton wavelength of the electron, however, is the fundamental length scale at which the charge 
of any bound particle is screened by the electronic vacuum polarization. 

The vacuum polarization (VP) correction to energy levels has been evaluated [5, 6, 134, 136, 137, 138] 
with nonrelativistic wave functions. We recall that the leading-order VP correction (due to the Uehling 
potential) can be expressed as 

(X 

t.E^(ii,\Vv\i,)^-CEE^,, (4.19) 



^.--^^ (4.20) 



where 

471^ 

is the Schrodinger binding energy for a two-body system with two particles each of mass m [first term on 
the right-hand side of (4.15)]. For the Ce coefficients, we recall the following known results [6, 136, 137], 



C£;(1S) = 0.22, C£;(2S) = 0.10. (4.21) 

Two-loop vacuum polarization effects enter at a relative order in pionium and are therefore of the 
same order of magnitude as the relativistic corrections mediated by the Breit interaction (discussed in 
Sec. 4.3). The self energy correction which is discussed in the following section is even smaller, but of 
considerable theoretical interest. 



4.5 Effects due to Scalar QED 

As shown in [lil2], the leading logarithmic correction to the self energy can be obtained, in nonrelativistic 
approximation, from second-order perturbation theory based on nonrelativistic quantum electrodynam- 
ics [104] (see also [103]). We will investigate here, in a systematic way, the leading logarithms generated 
for S states by self energy and rclativistic-recoil effects (the so-called Salpeter correction), and show that 
these are spin- independent. 

The quantized electromagnetic field is [see Eq. (5) of [102]], 

/d'^fc r 1 

/(27r)3 2fc ^^^^^ [a^^cxp(-ifc-r) + afc^exp(ifc-r)J , (4.22) 



A=l,2 ' 



and the nonrelativistic interaction hamiltonian for an atomic system with two spinless particles (charges 
ei and ei and masses m\ and m-i) reads 

i?i = - — Pi • A(ri) + -^A(ri)2 - • ^(ra) + -^A{v2f . (4.23) 

mi ZTTil 1712 27712 

For two spin-1/2 particles, the terms 

-^a,-Biri)-^cT2-Bir2) (4.24) 

77T.1 777,2 

have to be added to Hj [see Eq. (7) of [102]]. We will carry out the calculations for the general case of 
one particle of charge ei = e and the other having a charge 62 = —Ze (we follow the convention of [34] 
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that for hydrogen, e is the physical charge of the electron, i.e. e = — |e|). The unperturbed hamiltonian 
of the system of the two particles and the electromagnetic field reads [see e.g. Eq. (6) of [1(12]], 

where r = ri — r2. The eigcnstatcs of the "atomic part" Hq of this hamiltonian in the center-of- 
mass system Pi + P2 = are the nonrelativistic Schrodinger-Coulomb wave functions for a reduced 
mass ra-c = mim^/ {tni + 771,2) [here, the "atomic part" Hq excludes the photon field, i.e. the last term 
of (4.25)]. We denote p = p^= P2- 

Given that the first-order perturbation ((/is|-ffi|0s) vanishes, the second-order perturbation yields the 
dominant nonvanishing perturbation. When evaluated on an atomic state, it is given by 

5Es^ = ((f>s\Hi—^—Hi\^s) ■ (4.26) 

The interaction hamiltonian (4.23) gives rise to QED corrections that involve both particles (in the 
current context, these are recoil corrections involving the product ei 62), and also to terms which involve 
only a single particle and are proportional to ef or e^. The latter effects correspond to the self-energies 
of the two particles. 

The low-energy part of the self energy in leading order [32] can be inferred directly from (4.26), and it can 
be seen that the spin-dependent parts from (4.24) vanish in leading order in the (ZQ;)-expansion [102]: 

P 1 P 



6 7r^ 7o 



where 



mi - {Es - k) mi 



{Za)^'mi- 



</> ) + (ei ^ e2,mi <-> m2) . (4.27) 



Es--^^ (4.28) 

is the Schrodinger energy (rrir in the reduced mass of the atomic system under investigation). Starting 
from the spin-independent expression (4.27), it is now relatively straightforward to show that the leading 
"self energy logarithm" for S states is given by 



^4hi(z«r! 

3 TT n-^ 



(4.29) 



This result is by consequence spin- independent. The derivation is simplified when using the e- method 
developed and used in various bound-state calculations [1, 2, 32]. The two terms in square brackets 
in (4.29) correspond to the two self-energies of the two constituent particles with charges ei = e and 
62 = —Ze and masses toi and m2, respectively. It has been pointed out [97] that in contrast to the self 
energy corrections, the vacuum polarization corrections given in Eq. (4.21) must not be double-counted. 
The "double-counting" of self energy corrections (and lack of it in the vacuum-polarization case) finds 
a natural explanation in our formalism: whereas the vacuum-polarization correction mainly leads to a 
modification of the l/r-type Coulomb attraction in (4.25) within a nonrelativistic effective theory, the 
structure of the interaction hamiltonian (4.23) implies the existence of the two self-energies of the two 
constituent particles of the atomic system. 

It might be instructive to point out that the formula (4.29) is consistent with Welton's argument for 
estimating the self energy of a bound particle which is based on analyzing the influence of the fluctuating 
electromagnetic field [a detailed discussion is given on pp. 80-82 of [34]]. For a system with two particles 
of equal mass mi ~ 1112 = m, we have m^ = in/2. 

The leading-order recoil correction (Salpeter correction) can also be inferred from the interaction hamil- 
tonian (4.23) via second-order perturbation theory, by "picking up" terms that involve products ei 62- 
It has been shown in [102] that the leading logarithm (for S states) of the Salpeter correction is spin- 
independent (just like the leading logarithm of the self energy correction). The Salpeter correction is 
usually referred to as a relativistic recoil (RR) correction. By following [102], we obtain for the leading 
logarithm of this effect 
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This correction involves only the products ei 62 = —AirZa and can therefore be written as a function of 
Za alone. 

For pionium, we have Z = 1, mi = m2 = m = m^, = The leading logarithmic correction from 

scalar QED for pionium in the order of In a is obtained by adding the corrections (4.29) and (4.30), 



The non-logarithmic term of order is spin-dependent, and its evaluation requires a relativistic treatment 
of the self energy of a bound spinless particle; such a calculation would be of considerable theoretical 
interest, but the size of the effect for pionium, which is roughly two orders of a smaller than the leading 
vacuum polarization correction, precludes experimental verification in the near future. However, we would 
like to point out here that a fully relativistic treatment of this problem, including a detailed discussion 
of the rcnormalization of the self energy of the spinless particle, has not yet been accomplished. Scalar 
QED is a renormalizable theory [34] . 

The dominance of vacuum polarization over the self energy in pionium is expressed, in particular, by 
the fact that even two-loop vacuum polarization of order has a stronger effect on the spectrum of 
pionium than the leading logarithm from Eq. (4.31) according to Ref. [17], and that the strong-interaction 
correction of order [139] has to be well understood before any experimental verification of (4.31) appears 
feasible. Finally, we remark that for a manifestly non-elementary particle like the pion which has a finite 
charge radius, form- factor corrections have to be taken into account. 




(4.31) 
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Chapter 5 

QED Calculations: A Summary 



Calculations in the area of bound-state quantum electrodynamics have a long history, starting with Hans 
Bethe's first evaluation of the hydrogen Lamb shift [■>!]. Since then, the field has developed considerably, 
and the immense progress in high-resolution spectroscopy (for a recent culmination of the activities, 
see [i.-]) has necessitated the inclusion of higher- and higher-order corrections into the theoretical frame- 
work, leading to increased complexity of the theoretical calculations. Why should such an effort be carried 
out? A more accurate understanding of the QED corrections is of crucial importance for the determi- 
nation of the fundamental constants and the test of quantum electrodynamics. The theory of the QED 
bound states can, at little risk to over-statement, be described as one the most developed theories in all 
of physics. It combines the intricacies of quantum field theory (divergences, renormalization) with the 
experimental possibilities implied by the accuracy of high-resolution spectroscopy. 

For decades, there has been a tremendous amount of work on analytic self-energy calculations, notably 
on the coefficient Aqq for S states. However, an accurate comparison of the analytic results to numerical 
data was impossible, simply because no numerical data were available in the region of low nuclear charge. 
The semi-analytic expansion in Za is in part problematic because there is evidence that it may represent 
a divergent series. This argument may be based in part on the analogy between Figs. 5.1 and 5.2, and 
on the observed rapid growth of the analytic self energy coefficients {A40 is of order one for S states, 
A50 is about 10, and the magnitude of Aqq is as large as 30). The QED effective action gives rise to a 
divergent asymptotic expansion in appropriately defined effective coupling parameters. It is therefore a 
priori difficult to associate a finite remainder to a truncated series that involves all known (but disregards 
the unknown) analytic coefficients. 

The results reported here in Ch. 1 answer this conceptual question in the following way: the semi-analytic 
expansion in powers of Za is in full agreement with the obtained numerical results. We may therefore 
conclude that the semi-analytic expansion probably gives rise to a generalized asymptotic series whose 
remainder is of the order of the first neglected semi-analytic term as Za — > 0. Therefore, we may assume 
that a truncated semi-analytic expansion in powers of Za leads to reliable theoretical predictions in the 
realm of low nuclear charge numbers Z. However, although the numerical results confirm the validity 
of the semi-analytic expansion in this Z-region, they also show the limitations of this approach: much 
more accurate results arc obtainable by numerical calculations. The accuracy on the level of 1 Hz for 
atomic hydrogen reported in Ch. 1 would have involved analytic terms of an excessively high order in Za 
[approximately a (Za)^)], and the multitude and complexity of terms would have been a severe obstacle 
for analytic calculations in this order of Za. In general, the related activities have recently sparked a 
rather broad interest in the numerical calculation of rclativistic, QED self energy and two-body corrections 
at low Z and the comparison of analytic and numerical results [87, 140, 141, 142, 143, 144, 145, 14G, 147, 
148, 149, 150, 151]. 

We have described in Ch. 1 a nonperturbative evaluation of the one-photon self energy in hydrogenlike ions 
with low nuclear charge numbers Z = 1 to 5. The general outline of our approach is discussed in Sec. 1.3. 
In Sec. 1.4, the numerical evaluation of the low-energy part (generated by virtual photons of low energy) 
is described. In Sec. 1.5, we discuss the numerical evaluation of the high-energy part, which is generated 
by high-energy virtual photons and contains the formally infinite contributions, which are removed by 
the renormalization. Sec. 1.5 also contains a brief discussion of the convergence acceleration methods 
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Figure 5.1: Expansion of the bound electron self energy in powers of the binding Coulomb field. 
The dashed lines denote Coulomb photons, the crosses denote the interaction with the (external) 
binding field. Each interaction with the binding Coulomb field gives rise to a higher-order term in 
the Za-expansion. 



as employed in the current evaluation. We discuss in Sec. 1.6 the comparison of analytic and numerical 
data for K- and L-shell states in the region of low Z. The main results of the current investigation are 
contained in Table 1.5: numerical data, nonperturbative in Za, for the scaled function F and the self 
energy remainder function Gse for K- and L-shell states at low nuclear charge. The numerical accuracy 
of our data is 1 Hz or better in frequency units for IS, 2S and both 2P states in atomic hydrogen. 

Still, calculations within the semi-analytic expansion in addition to the numerical calculations (a kind of 
"dual strategy") remain important. They serve as a confirmation and check for currently available and 
future numerical approaches. In general, in view of the complexity of the calculations, an independent 
check is always highly desirable. For the two-loop effect, the (ZQ;)-expansion converges more rapidly than 
for the one- loop effect in absolute frequency units because of the additional radiative factor a/n which 
decreases the overall size of the effect (see Ch. 3) . 

It is hoped that the analytic calculations reported in Ch. 3 for low nuclear charge number Z will be sup- 
plemented in the future by an accurate numerical treatment of the two-loop self energy problem (see also 
related recent work in the high-Z region, Refs. [13, 151, 152]). This presupposes that the considerable nu- 
merical problems in the domain of small nuclear charge could be solved by adequate numerical methods, 
and that the further problem of the increased computational demand of the two-loop effect in comparison 
to the one-loop problem [9, 12] can be tackled - possibly by massively parallel computer architectures. 
Note, however, that the most accurate theoretical predictions could only be reached in combining nu- 
merical and analytic results. The reason is the following: All numerical calculations are performed in the 
non-recoil limit which is the limit of infinite nuclear mass. This is not quite sufficient for an accurate 
theoretical treatment because the self energy of a bound-state depends genuinely on the ratio of the or- 
biting particle to the nuclear mass - an effect beyond the recoil correction. For example, the argument of 
the logarithms in Eqs. (1.2) and (3.1) should be replaced according to ln[(Za)~^] — *■ ln[cr (Za)^^], where 
a = m/n\ and n\ is the reduced mass [!)?]. The possibility to include these tiny, but important effects 
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Figure 5.2: Expansion of the QED effective action in powers of the external field. As discussed 
in Sec. 8.3.2, the resulting asymptotic series in powers of the coupling is divergent. 



depends crucially on a reliable knowledge of the analytic coefficients in combination with an accurate 
numerical treatment of the problem (see also the discussion in Sec. 2.6). 

The analytic results can be used to obtain improved theoretical predictions for the hydrogenic fine struc- 
ture as compared to the previous order-a'^-calculations [1, 2], because they remove the principal theo- 
retical uncertainty in the order of due to the problematic two-loop self energy which is represented 
diagrammatically in Fig. 3.1. Our calculation in Ch. 3 illustrates the usefulness of the simplified effective 
treatments of two- loop effects in the analytic approach based on the modified Dirac hamiltonian (3.3) 
and the "e method" . This aspect highlights, as we believe, the need for systematic, simplified treatments 
of higher-order radiative corrections in bound systems. 

In Ch. 3, we primarily address spin-dependent effects in one-electron (hydrogenlike) systems. However, 
the same effects also contribute to the fine-structure splitting in two-electron (heliumlike) systems. There 
is currently remarkable interest in improved measurements of the fine-structure splitting in helium and 
heliumlike atomic systems with low nuclear charge [3(), 37, 38, 39]. The effects addressed here contribute 
to the fine-structure splitting in helium on the level of 100 Hz, which is not much smaller than the 
current experimental accuracy of about 1 kHz, and allows for an estimate of uncalculated higher-order 
contributions. 

All analytic calculations rely on the separation of the virtual photon energy into high- and low-energy 
energy regions via a parameter e. This separation is illustrated in Sec. 2.2. For two-photon problems, 
four energy regions result from such a separation, but for the one-photon problem discussed in Ch. 2, 
one separation parameter e is sufficient. We present in Ch. 2 analytic calculations for excited atomic 
states based on the semi-analytic Za-expansion. These may be complemented in the near future by 
numerical results [95]. We discuss some intricacies of the analytic treatment of highly excited atomic 
states in Sec. 2.1. Results for the individual contributions are presented in Subsections 2.3.1 - 2.3.8 
and summarized in Sec. 2.4. Final values of the analytic coefficients can be found in Tables 2.1 - 2.7. 
The characteristic numerical cancellations between individual contributions are described in Sec. 2.5. 
Finally, in Sec. 2.6, we describe a number of observations regarding the dependence of the Aqi and the 
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AgQ-coefficients on the quantum numbers of the atomic state. 

All of the above calculations are related to the bound electron, which is a spin-1/2 particle. In nature, 
particles with different spin quantum numbers exist, and these, too, can form bound systems. The spin- 
dependence of the corrections is a conceptually interesting issue that forms the subject of Ch. 4. This 
interest is sparked by the projected DIRAC experiment at CERN where the bound system of two oppo- 
sitely charged pious ("pionium") will be studied. Notably, even the rclativistic corrections for a bound 
two-body system consisting of spinless particles are different for spinless particles than for particles with 
nonvanishing spin. They are given by a generalized Breit hamiltonian that we discuss in Sec. 4.3. There, 
we present a simplified derivation for the rclativistic and recoil corrections of order to a bound state of 
two spinless particles. The results agree with previous calculations [f2!)]. As evident from equation (4.14), 
the zitterbewegung term is absent in a bound system of two spinless particles. 

The self energy is suppressed in systems with spinless particles in comparison to the vacuum polarization 
effect as discussed in Sees. 4.4 and 4.5, because the lightest known spinless particle is much heavier than 
the electron, which implies that the vacuum polarization effect is larger by two orders of Za than the self 
energy in bound systems of spinless particles. We provide a complete result for the leading scalar-QED 
correction of order In a in Sec. 4.5. A list of further QED corrections to the IS level of pionium, 
including two- loop vacuum polarization, finite-size corrections and other effects, can be found in [17]. 

To conclude this Chapter, we would like to mention the complementarity of the high-precision calculations 
and experiments at low nuclear charge and the strong-field calculations and related experiments at high 
nuclear charge numbers. While an accuracy [42] of one part in 10^^ currently appears to be out of reach for 
hydrogenlike Uranium [15.3], the experiments on highly charged ions have meanwhile attained a precision 
where a detailed understanding of the nuclear structure and shape is necessary in order to describe QED 
effects rigorously [154]. 

Bound systems described by quantum electrodynamics belong to the most accurately understood physical 
systems today. They convey a certain esthetic appeal because the apparent simplicity of the two- and 
three-particle systems gives rise to a large number of conceptual as well as calculational difficulties when 
a thorough understanding is sought in higher orders of perturbation theory. 
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Part II 



Convergence Acceleration and 
Divergent Series 
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Chapter 6 



Introduction to Convergence 
Acceleration and Divergent Series in 
Physics 



The second part of this thesis will be concerned with convergence acceleration methods and with the 
summation of divergent series. In theoretical physics, we often encounter sequences of numbers, either as 
a result of perturbation theory or as partial sums of an infinite series defining a physical quantity, whose 
numerical properties are disfavourable: in many cases, the sequences converge so slowly that they are 
numerically useless. In other cases, the sequences do not converge at all, and yet we know ~ or assume - 
that the sequences are associated, nevertheless, to the physical quantity in question, for example in the 
form of a divergent asymptotic expansion. 

Let us consider the case where the sequence {s„}^g under investigation is given by partial sums of an 
infinite series, 

71 

Sn^^ttk- (6.1) 

fc=0 

Wc assume, furthermore, that the are real (this condition will be relaxed in the sequel). 

(i) Let us assume that the sequence {s„}^q is convergent. The purpose of convergence acceleration is to 
convert {s„}5^q into a new sequence {s'„}^^q with hopefully better numerical properties. If we denote 
the limit of the sequence by s, that is s = lim„_,oo Sn, then the rate of convergence is obviously increased 
if 

lim = . (6.2) 

(ii) Alternatively, let us assume that the sequence {s„}^q is divergent. The purpose of the resummation 
is to determine the generalized limit (or antilimit) of the sequence. The following question immediately 
arises: how is it possible to associate a finite generalized limit with a divergent scries. L. Eulcr wrote in a 
letter to Goldbach (1745): "Summa cuiusque seriei est valor expressionis illius finitae, ex cuius evolutione 
ilia series oritur." , which means: the generalized limit of any divergent series is the value of the particular 
characteristic finite expression, whose expansion gives rise to the divergent series. This is still the main 
notion of resummation today: Given a divergent series, to "reconstruct" a posteriori the function whose 
expansion gave rise to the divergent series, with a suitable mathematical resummation method. 

We illustrate this idea with a trivial example: A divergent series can be generated, e.g., by expanding 
the function f{z) = 1/(1 + 2:) ~ J2'jLoi^'^y at z = 1 in powers of z. The partial sums of the resulting 
divergent series read 

" I 1 for n even 

^" = E(-l)' = i ■ (6-3) 

j=o [ for n odd 

The resulting sequence, which is 1, 0, 1, 0, ... , can be identified as originating from the discussed fimction 
by expansion at unit argument, and its antilimit is therefore Since this argument may not sound 
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convincing, let us briefly indicate that the result may also be obtained by imposing the postulate that 
the antilimit of the sum of two divergent series should be the sum of the two antilimits each. Then, it 
is easy to derive (by inspection) that the antilimit a of the series X]^o(~l)'' should fulfill 1 — a = a, 
and thus a — ^/^ (again). A more formal investigation of this series and its resummation can be found 
in [155] [see Eq. (1.1.1) and the instructive discussion on the p. 6 ibid]. It is probably a triviality to 
express that the most familiar resummation methods like the Euler transformation, Pade approximants 
and recently developed nonlinear sequence transformations of the Levin- Wcniger type all convert the 
divergent sequence 1,0,1,0,... to the convergent sequence Va^ V27 V21 Vs: ■ • • > in mutual agreement. 
Some of these resummation methods are discussed in the sequel, and we also refer the reader to the 
seminal review [15()]. It has been stressed in J. E. Littlewood's introduction to [155] that the subject of 
divergent series is in no way "mystical or unrigorous" . 

Fundamental differences exist between nonalternating series (without loss of generality we can assume 
flfe > Vfc in that case), and strictly alternating series Ofe — ("I)''' bk where bk > yk (mixed forms will 
also be discussed in the sequel). The acceleration of the convergence and the summation of divergent 
nonalternating series has been historically problematic. 

For example, it has been known for a long time that the acceleration of nonalternating slowly convergent 
series can be highly unstable numerically (see [15G] and references therein). Likewise, the summation of 
nonalternating factorially divergent series has represented a considerable problem, because the series are 
typically non-Borel summable (see [157, 158, 159, 160] for suitable generalizations of the concept of Borel 
summability which provide solutions for large classes of problems of practical interest). 

Recently, considerable progress has been achieved for the nonalternating case, and we will describe some 
of the relevant applications. For the acceleration of convergence of nonalternating series, new power- 
ful algorithms have been devised [23], and generalizations of the Borel summation process have led to 
promising results in the case of nonalternating factorially divergent series [25, 27, 28]. 

One would intuitively assume that the perturbation series, even if it is divergent, should be associated to 
the physical quantity of interest in a more or less unique way. Unfortunately, many divergent series and 
the associated physical observables (e.g., energy levels as a function of the coupling strength) fail to fulfill 
the Carleman criterion [ I (iJ] which guarantees that an asymptotic series can be uniquely associated to a 
given mathematical function. It appears that some ambiguity has to accepted as far as the connection 
between perturbation theory and physical quantities is concerned, although some attempts have been 
made to resolve ambiguities in particular cases (see e.g. [1()2]). 

For completeness, we give here the Carleman criterion which guarantees that there is a one-to-one cor- 
respondence between a function and its associated asymptotic series (see for example [Ki.':)], Theorems 
XII. 17 - XII. 19 and the definition on p. 43 in [1G4], p. 410 in [105], or the comprehensive and elucidating 
review [16G]): 

Carleman Criterion. Let f{z) be a function which is analytic in the interior 
and continuous on a sectorial region 5 = {z|| arg(2;)| < /c7r/2 -|- e,Q < |z| < i?} 
of the complex plane for some e > 0. Let the function / have an asymptotic 
expansion f{z) ~ X^J^Lo Cn 2;" (for z 0). The function / obeys a strong 
asymptotic condition (of order k) if there are suitable positive constants C and 

a such that \f{z) - i2n=0 -^"1 < C'cr'^+l (^-^ ^ |^|m+l j^^^jg f^j. j^jj ^ 

and for all z £ S. The validity of such a condition implies that the function 
f{z) is uniquely determined by its asymptotic series. 

Typically, nonalternating factorially divergent series which entail nonperturbative (imaginary) contribu- 
tions do not fulfill the Carleman condition. This failure gives rise to the inevitable ambiguities. 

Even the generalized Borel resummation methods introduced in [157, 158, 159, IGO] cannot provide 
conclusive answers for all possible physical scenarios. One highly problematic application is given by the 
quantum-mechanical double-well oscillator. The perturbation series for the ground state (as well as for 
any other unperturbed state) is divergent. Yet it is known that the ground state splits into two states of 
opposite parity when the second minimum of the potential approaches the first from infinity [167, 168]. 
The interesting observation is that two distinct quantum states are described by one and the same 
perturbation series. The energy difference is nonperturbative in the coupling and cannot be inferred in 
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principle from perturbation theory. The problem has been the subject of intensive discussion [169, 170, 
171, 172, 173, 174]. The resulting multi-instanton expansion has recently found a more mathematically 
motivated explanation in the theory of resurgent functions [175, 176, 177]. 

We conclude this introduction by expressing that convergence acceleration and resummation methods 
cannot be separated. The connection can be seen as follows. We consider the remainders r„ = .s„ — s, 
i.e. the difference of the element s„ of the sequence and the limit s. The acceleration of convergence is 
achieved by eliminating, as far as possible, the remainders r„ from a necessarily finite set of input data 
{s„}™^g. in order to obtain a better estimate s^^ to s than the obvious estimate which would be provided 
by Sm- Clearly, the "elimination of remainders" is a concept which can be transported to the case of a 
divergent series. The only difference is that for a convergent series, r„ — > as n — > oo, whereas for a 
divergent series, lim„^oo fn does not exist. Otherwise, the concept of the elimination of remainders is the 
same for both cases. 

Returning to our (simple) example, it is 
remainders 

inside the circle of convergence \z\ < 1, could be suitable to accomplish an analogous elimination of 
the r„ if z assumes values outside of the domain of convergence of the series '^JLoi~^y because the 
mathematical structure of the expression which defines r„ docs not change. These entirely heuristic 
considerations will be illustrated and explained in the sequel. 



intuitively obvious that a method which helps to eliminate the 

- E (-^y - ^iT^ (6-4) 
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Chapter 7 



Convergence Acceleration 



7.1 The Concept of Convergence Acceleration 
7.1.1 A Brief Survey 

Wc intend to give a brief overview of the convergence acceleration methods used in this Thesis, without 
any claim for completeness or comprehensiveness with regard to currently known convergence acceleration 
techniques. Indeed, excellent textbooks (e.g. [178, 179, 180]) and review articles [156, 181, 182, 183, 184] 
on the subject are available. Here, we focus on aspects of the methods which will become useful for the 
applications discussed in later Chapters of the Thesis. 

This brief introduction is structured as follows: After a discussion of basic concepts (Sees. 7.1.2 and 7.1.3), 
we introduce Pade approximants (one of the standard mathematical tools for extrapolation and conver- 
gence acceleration) in Sec. 7.1.4. Next, wc discuss nonlinear sequence transformations in Sec. 7.1.5. Both 
of these convergence acceleration methods arc much more powerful and much more numerically stable 
in the case of alternating input series than in the case of input series whose terms have the same sign 
(nonaltcrnating series). We conclude in Sec. 7.1.6 with a discussion of the recently proposed [2'.]] combined 
nonlinear-condensation transformation (CNCT) which provides an efficient algorithm in the historically 
problematic case of nonalternating input series. 

Before we start the discussion of Fade approximants and nonlinear sequence transformations, we would 
like to mention that a large number of other convergence acceleration methods exist which have all been 
documented in the literature: examples include the Euler transformation, Aitken's process, the iterated 
process, and Richardson extrapolation which can be implemented using Wynn's rho algorithm, and 
its iterations. Last, but not least, we mention the Euler-Maclaurin sum formula. The interested reader 
will find related information e.g. in Chs. 5 and 6 of [15(3], App. E.2 of [4.3], and in Eq. (3.6.28) ff. of [185]. 

Concerning the mutual differences between the Fade method (Sec. 7.1.4 and nonlinear sequence trans- 
formations (Sec. 7.1.5), we would like to expand: Given a formal series 

oo 

as input, the Fade method produces rational approximants which depend exclusively on a finite subset of 
the coefficients {cj}°°^Q, which in turn determine the partial sums of (7.1). Nonlinear sequence transfor- 
mations require, in addition to the partial sums of (7.1), also explicit remainder estimates, i.e. estimates 
for the remainder X^j^n+i '^j which remains after the summation of the first n terms of (7.1). 

Thus, the remainder estimates introduce additional degrees of freedom in the construction of sequence 
transformations as compared to Fade approximants. As pointed out in [2()], one may draw an analogy 
between sequence transformations and Fade approximants on the one hand and the Gaussian integration 
and the Simpson rule on the other hand; the variable integration nodes and weight factors of the Gaussian 
integration yield additional degrees of freedom which may be used in order to construct a potentially 
much more powerful algorithm for numerical integration. 
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Except for the additional degrees of freedom, nonlinear sequence transformations have additional advan- 
tages. For example, the delta transformation to be discussed in Sec. 7.1.5 gave very good results in the 
case of the sextic anharmonic oscillator, and there is strong numerical evidence that it is able to sum the 
extremely violently divergent perturbation series for the octic anharmonic oscillator [186]. By contrast, as 
discussed in [180] , Fade approximants are apparently not powerful enough to sum divergent series whose 
coefficients diverge like (2n)! or even (3n)!. These observations stimulated research on sequence transfor- 
mations considerably. The rapid progress in this field is convincingly demonstrated by the large number 
of monographs [180, 187, 188, 189, 190, 191, 192, 193, 194, 19-5] and review articles [156, 196, 197, 198] 
which appeared in recent years. 



7.1.2 The Forward Difference Operator 

An essential ingredient in the construction of sequence transformations, both for convergence acceleration 
and for the resummation of divergent series, is the (forward) difference operator. We define the difference 
operator A acting on the sequence {s„}^o- 

A : {s4- {{Mn}^=o > (As)„ = Sn+l - Sn , (7.2) 

The nth element of the sequence {(As)„}^q is given by (As)„ = Sn+i — s„. This way of writing the nth 
element of the transformed sequence stresses the fact that (As)„ is to be regarded as an element of some 
new sequence. However, the brackets are often left out, and we write 



As„ = s„+i - Sn = (As)„ . (7.3) 

(A'^5)„ = (-1)'^ ^(-I)^YJ) . (7.4) 

A'=s„ = (A'=s)„ . (7.5) 
The fcth power of the difference operator plays a crucial rule in constructing sequence transformations. 



An important relation is. 



Again, we identify 



7.1.3 Linear and Logarithmic Convergence 

As already stressed in Ch. 6, the main notion of convergence acceleration is to extract information "hidden 
in trailing digits" from a necessarily finite number of sequence elements, in order to convert a sequence 
{s„}5^g into a new sequence {s^^o with hopefully better numerical properties. 

We briefly recall the notion of logarithmic convergence and the difficulties associated with the acceleration 
of the convergence of nonalternating series. Let a sequence {sn}^o fulfill the asymptotic condition 

lim ^i+i =: (7.6) 

n^oo Sn ^ S 

where s = Soc is the limit of the sequence as n ^ oo. If p > 1, then the sequence {snli^o divergent. 
For p = 1, the sequence may either be convergent or divergent. A convergent sequence with p = 1 is 
called logarithmically convergent (if p < 1, the series is called linearly convergent). Let us further assume 
that the elements of the sequence {sn}^=o in Eq. (7.6) represent partial sums 

n 

s„ = ^a(fc) (7.7) 

fc=0 

of an infinite series. Here, we will almost exclusively investigate slowly convergent nonalternating se- 
quences {Sn}^=o whose elements are all real and positive (for these sequences, < p < 1). In the case of 
slow convergence, p is either very close or equal to unity. 

As observed by many authors (e.g. [1">()]), the acceleration of the convergence of nonalternating sequences 
is a potentially unstable numerical process. The reason is the following: A sequence transformation can 
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only accelerate convergence if it succeeds in extracting additional information about the index-dependence 
of the truncation errors 

Tn = Sn- S (7.8) 

from a necessarily finite set of partial sums {sn}n=o input series. Normally, this is done by forming 

arithmetic expressions involving higher weighted differences of the s„. The calculation of higher weighted 
differences is a potentially unstable process which can easily lead to a serious loss of numerical significance 
due to cancellation if the input data all have the same sign. 

We had mentioned in Ch. 6 that the rate of convergence can be said to be increased if the transformed 
sequence elements {s^}J5^o ^^"^ ^^le input data {s„}J^g fulfill [see also Eq. (6.2)] 

s' -s 

lim T„ = lim , (7.9) 

n — *oo n — >oo Sji — 5 

where we implicitly define the r„. Let us assume that the input sequence {s„}J^q is logarithmically 
convergent and that the transformed sequence {s^}^o hnearly convergent with lim„_^oo(sJi+i— s)/(s^ — 
s) = p' < 1 [see Eq. (7.6) and the definition of logarithmic convergence]. Then 

lim = hm fki^ Jn_^ =p' <1, (7.10) 

n^oo r„ n^oo s'^ - S S„+i - S 

SO that lim„_,oo Tn ~ 0, and the rate of convergence is increased [according to the definition of convergence 
acceleration in Eq. (7.9)]. 

7.1.4 The Standard Tool: Pade Approximants 

The Fade approximation produces a rational approximation to the original input series; this approxima- 
tion is given by the ratio of two polynomials. The [l/m] Pade approximant to a series 'P{z), 

oo 

7'(z) = 5]c,^^ (7.11) 

is given by the ratio of two polynomials Pi{z) and Qm{z) of degree / and m, respectively: 

ll/m]v{z) ^ = — . (7.12) 

Qm{z) go + 91 2 + . . . + (?m z™ 

The polynomials Pi{z) and Qm{z) are constructed so that the Taylor expansion of the Pade approximation 
Eq. (7.12) agrees with the original input series Eq. (7.11) up to the order of Z + to 

r{z) - [l/m]v{z) = 0(z'+™+i). (7.13) 

In other words, the Pade approximant [l/m]-p{z) reproduces upon re-cxpansion the first / + m + 1 
(co, . . . , ci+m) coefficients of the input series. 

Usually the normalization condition is choosen such as qo is one. Hence there are / -f m -I- 1 equations to 
determine the coefficients of Pi{z) and Qm{z)- In most cases it is quite unattractive to solve this system 
of equations. A more attractive recursive prescription for the evaluation of Pade approximants is the 
e-algorithm given by the formulas [199] 

,(") _ A ,(") _ „ 



(n) _ (n+1) , ,r (ti+1) (n) 



.-.1 - ^.-1 + i/[er^^-4"i (7.14) 

for k and n integer. The e-algorithm is discussed in [199], and a detailed discussion of its implementation 
can be found in [156]. If the input data Cq"^ = s„ are the partial sums 

n 
j=0 
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of the power series Eq. (7.11), then the elements e^^k of the two-dhxicnsional table of transformations 
produce Pade approximants according to 

e^",) = [n + k/k]p{z). (7.16) 

In) 

The elements with odd lower subscripts e-2k+i tabic of transformations arc only auxiliary quantities. 

In practical applications (see e.g. Tables 1 - 3 of [24]), the epsilon algorithm can be used for evaluation 
of the following staircase sequence in the Pade table: 

[0/0], [1/0], [1/1], [2/1], . . . , [v/v], [y + l/y], [y + + . (7.17) 

Because of Eq. (7.16) this corresponds to the sequence 

^(0) ^(1) JO) (1) (0) (1) (0) (0) (1) (0) /„-,ox 

^0 '^0 ' 1*^2 1*^4 1^4 1^6 I ■ ■ ■ ^21^ ; ^2iy ; ^2iy+2' ■ ■ ■ i \'-^°) 

which can be written as {e2[i//2V^^"''-'''^0' ^'^^'^ ll'^/^l denotes the integral part of i^/2, which is the largest 
integer ^ satisfying fi < iy/2. With this notation the Pade sequence in Eq. (7.17) can also be written as 
{[i.-H2l/H2l]p(z)}-o. 

Modern computer algebra packages [90] contain built-in routines for the evaluation of Pade approximants. 
While these are computationally less effective than the epsilon algorithm, they are more flexible with re- 
gard to the degree of the numerator and denominator polynomials. In particular, the sequence allow for the 
evaluation of the sequence of lower-diagonal Pade approximants of the form {[[Jz//2]]/z/ — [[i^/2]]p(2;)}J^g. 

As already mentioned in Sec. 7.1.1, it is a typical feature of all algorithms for Pade approximants that the 
coefficients cq, ci, . . . , q+,„ of the formal power series (7.11) suffice to determine [l/in]-p{g) completely. 
No additional information is needed. At first sight, this seems to be very advantageous. However, there are 
situations in which this apparent advantage becomes a disadvantage. For instance, the truncation error 
of an asymptotic power series which has the Stieltjes property is for real argument bounded in magnitude 
by the first term not included in the partial sum (see Example 9.3 on p. 92 of [91]). Accordingly, the 
truncation error of a strictly alternating divergent hypergeometric series 2Fo{a, f3; —x) with a,(3,x > 0, 
which provides a good mathematical model for many divergent perturbation expansions, is bounded in 
magnitude by the first term not included in the partial sum. In addition, the truncation error and the 
first term not included in the partial sum have the same sign (see Theorem 5.12-5 of [200]). In principle, 
this information should be very helpful in summation processes. Unfortunately, Pade approximants are 
not able to profit from it. This greatly reduces their efficiency if a divergent hypergeometric series 2-F0 of 
that kind or a divergent perturbation series of a similar type is to be summed. 



7.1.5 Nonlinear Sequence Transformations 

A sequence transformation T is a rule which converts a slowly convergent or divergent sequence {s„}J^q 
whose elements may be the partial sums .s„ — X]fe=o ^'^ I^*^*^ (^-1)] '^^ infinite scries, into a new 
sequence {s^jj^g with hopefully better numerical properties. Pade approximants, which were discussed 
in Section 7.1.4, are special sequence transformations since they transform the partial sums of a formal 
power series (7.11) into a doubly indexed sequence of rational approximants (7.17). In the current Section, 
alternative transformations will be discussed. 

Let us now assume that the sequence {s„}^q of partial sums of an infinite series diverges but can be 
summed to its generalized limit s. Then the sequence elements s„ can for all integers n > be partitioned 
into the generalized limit s and remainders r„: 

Sn = s + r„ . (7.19) 

If {sn}^=o diverges, then the remainders r„ do not vanish as n —^ 00. A divergent sequence would be 
summed to its generalized limit s if the remainders r„ could somehow be eliminated from the sequence 
elements s„. For a slowly convergent sequence, one can hope to accelerate convergence by eliminating 
approximations to the r„ obtained from a finite set of sequence elements which are used as input data 
for the transformation. 
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The sequence {snjJ^Lo transformed into a new sequence {s[^}^^q, whose elements can for all n > be 
partitioned into the generalized limit s of the original sequence and a transformed remainder r^: 

< - s + r; . (7.20) 

This approximative elimination process can be termed successful if the transformed remainders either 
vanish as n ^ oo (in the case of a divergent input sequence where the r„ do not vanish at all as n — > oo), or 
if the vanish faster than the original remainders r„ as n — oo [in the case of a slowly convergent input 
sequence, see also Eq. (7.9)]. A sequence transformation is a rule which transforms a slowly convergent 
or divergent sequence {sn}'^=o into a new sequence {s^jj^o with hopefully better numerical properties. 
However, the truncation errors of a divergent series can depend on n in a very complicated way. Con- 
sequently, the direct elimination of approximations to r„ from .s„ can be very difficult. A considerable 
simplification can often be achieved by means of a suitable reformulation. Let us consider the following 
model sequence (see Section 3.2 of [loG]): 

S„ = S + LOn Zn ■ (7.21) 

Here, LJn is a remainder estimate, which has to be chosen according to some rule and which may depend 
on n in a very complicated way, and z,i is a correction term, which should be chosen in such a way that 
it depends on n in a. relcitively smootli Wciy. Moreover, the products cUjiZji should be capable of producing 
sufficiently accurate approximations to the actual remainders r„ of the sequence to be transformed. 

The principal advantage of this approach is that now only the correction terms {^nl^o have to be 
determined. The subsequent elimination of oJnZn from s„ is often much easier than the direct elimination 
of approximations to r„ from s„. The use of remainder estimates {w„}^q is also an efficient way of 
incorporating additional information about the remainders into the summation scheme. 

The model sequence (7.21) has another undisputable advantage: There is a systematic way of constructing 
a sequence transformation which is exact for this model sequence. Let us assume that a linear operator T 
can be found which annihilates the correction term z„. Then, a sequence transformation, which is exact 
for the model sequence (7.21), can be constructed quite easily. Just apply T to [s„ — s]/w„ = z„. Since 
T annihilates Zn according to T(z„) = and is by assumption linear, we find that the following sequence 
transformation T is exact for the model sequence (7.21) (see Eq. (3.2-11) of [15(i]): 

T{Sn,UJn) = - — - = s. (7.22) 
T[l/UJn) 

Simple and yet very powerful sequence transformations are obtained if the annihilation operators are 
based upon the finite difference operator A defined in Sec. 7.1.2. A polynomial Pk-i{n) of degree k — 1 
in n is annihilated by the k-th power of A. Thus, we now assume that the correction terms {zn}^^Q can 
be chosen in such a way that multiplication of z„ by some suitable quantity Wk{n) yields a polynomial 
Pfe_i(n) of degree fc — 1 in n: 

Wk{n)zn = Pk-i{n) . (7.23) 

Since A'=Pfc_i(n) = 0, the weighted difference operator T = A'^ii;fe(n) annihilates Zm and the correspond- 
ing sequence transformation (7.22) is given by the ratio 

rr(n)i I M N A'={u;fc(n)s„/tJ„} 

^ ' [wk{n)\sn,'jJn) = —J- — ^. (7.24) 

A number of sequence transformations can be obtained by specializing Wk{n). For instance, Wkin) = 
(n + Z?)*"'"^ with f3 > yields Levin's sequence transformation [201]: 

A^{(n + /3)^-i.„/^4 
'= ^^'^"''^"^ - AM(" + /3)^-V'^n} 



k 

E 



i-iy 



k\ (P + n + j) 



k-l 



jj i(3 + n + k)''~^ uj„+j 



EM) 

3=0 



jJ (P + n + k)''-^ 



(7.25) 
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We follow the notation of [ISG]. The subscript k denotes the order of the transformation, and n denotes the 
starting index (the initial element) of the input data. Observe that the fc + 1 sequence elements s„, ... Sn+k 

(n) 

and the + 1 remainder estimates w„ . . . ujn+k are needed for the computation of the transformation . 
The shift parameter /3 has to be positive in order to admit n = in (7.25). In most practical applications, 
it is observed that any choice other than /3 = 1 does not lead to a significant improvement of the rate of 
convergence, and therefore the choice /? = 1 is standard [156]. 

The notation £^"'(/3, s„, Wn) is in need of a certain further explanation. As explained above, n represents 
the initial clement of the sequence {snjJ^Lo used in the evaluation of the transform. The specification of 
s„ as the second argument of £[,"''(/3, s„, w„) is therefore redundant as far as the index n is concerned, 

because n already appears as an upper index in >C^"''. The arguments s„ (and a;„) are to be interpreted 
as follows: they rather specify the initial element s„ and the initial remainder estimate w„ which have 
to be used in the calculation of the right-hand side of (7.25). Of course, the evaluation of £^"^(/3, s„,a;„) 
requires the k + 1 sequence elements {s„, . . . , Sn+k, Sn+fc} as input. The knowledge of the "post-initial" 
elements {sn+i, . . . , Sn+k} is implicitly assumed in writing the definition (7.25). 

According to Smith and Ford, who had investigated the performance of sequence transformations for 
a wide range of test problems [202, 203], Levin's transformation is among the most powerful and most 
versatile sequence transformations that are currently known. 

We can also assume that the unspecified weights Wk{n) in (7.24) are Pochhammer symbols according to 
Wk{n) = [n + f})k-i ~ r(n + (3 + k — 1)/T{n + (3) with /3 > 0. This yields the delta transformation (see 
Eq. (8.2-7) of [15(i]): 

, _ A'^{(?l + /3)fc,i Sn/oJn} 

A'={(ri + /3)fc_i/a;„} 



k\ + n + j)k^i Sn+j 



j=0 ^-^ ^ 



E (-1) 



jJ iP + n + k)k-l LOn+j 

. (7.26) 



(/3 + n + j)fc- 



jj iP + n + k)k-l LOn+j 



As in the case of Levin's transformation, the most obvious choice for the shift parameter is /3 = 1. 

The numerator and denominator sums in (7.25) and (7.26) can be computed more effectively with the 
help of the three-term recursions [see Eq. (7.2-8) of [1 "><>]]: 

r{n) ir.. _ An+1)(0) (/? + n + fc) (/3 + n + fc)^'! („) 

Lk+iW-Lk + ^ + + (7.27) 

and (sec Eq. (8.3-7) of [15(i]) 

c(") - {P + n + k){p + n + k-i) („) 

^.+1 - S, (/3) - (^ + „ + 2fc)(/? + n + 2/fc-l) " ^^■^^> 

The initial values Lq"^(/3) = S'q"^(/?) = s„/cj„ and Lq" ■*(/?) = 5'g"^(/3) = l/u}„ produce the numerator and 
denominator sums, respectively, of £^"'' (/3, s„, Wn) and iS^"'^(/3, s„, a;„). 

Let us assume that the terms in (6.1) are strictly alternating, i.e. ak = (— l)'^6fc with bk > 0. In this 
case, the first term (— l)"+^&„+i not included in the partial sum s„ = X]fe=o(~-'^)'^^'= ^^^'^ hest simple 
estimate for the truncation error. This corresponds to the remainder estimate 

= (-l)"+i5„+i = As„, (7.29) 

which was proposed by Smith and Ford [2Uli]. The use of this remainder estimate in (7.25) and (7.26) 
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yields the following sequence transformations [Eqs. (7.3-9) and (8.4-4) of [156]]: 

4")(As„) = 4"'(/3,s„,As„) 



k 



J=0 



jj (/3 + n + fc)''-i As„+j 



fc 



(7.30) 



3=0 

4"H/3,s„) = 5(")(/3,s„,As„) 



k 

E 

J=0 



k\ {f3 + n + j)k-i s„+j 



jJ {f3 + n + k)k-i As„+j 



fc 

E (-1)^ 

J=0 



k\ {(i + n + j)k-i 1 



(7.31) 



j7 {P + n + k)k-i As„+j 



Alternative remainder estimates for the sequence transformations (7.25) and (7.26) were discussed in 
Sections 7 and 8 of [156] or in [204, 205]. 

A final word on sequence transformations: It is observed that sequence transformations constructed 
according to (7.31) often lead to much better numerical results than better known convergence ac- 
celerators like the Euler transformation. Wynn's epsilon or rho algorithms, or Aitken's A^ pro- 
cess [1")0]. In particular, there is overwhelming empirical evidence that d^"'' and S]^^ work very 
well in the case of convergent or divergent alternating series for instance as they occur in spe- 
cial function theory [23, 156, 180, 201, 202, 203, 20G, 207, 208, 209, 210] or in perturbation theory 
[186, 207, 208, 211, 212, 213, 214, 215, 216, 217]. Only for some special model problems could rigorous 
convergence proofs be obtained (see [218] or Sections 13 and 14 of [156]). No general convergence proof 
is known for delta transformation (5^"-* as well as their parent transformations and In view of 
the partially unsatisfactory situation concerning the availability of mathematical proof, we have carried 
out high-precision studies of the convergence of sequence transformations of the type (7.31), in part using 
multi-precision arithmetic (see Sec. 7.2.2). In the absence of a rigorous proof, one of the concerns which 
may be raised against the convergence of the transforms (7.31) is a plausible asymptotic nature of the 
sequence of transforms. In our numerical experiments, we found no indication of asymptotic behaviour. 
The rate of convergence of the transforms remained constant, and we observed apparent convergence 
up to 430 decimal figures (see Sec. 7.2.2 below). In the absence of rigorous proof and the presence of 
considerable numerical evidence for the computational usefulness of sequence transformations of the type 
(7.31), we believe that experimental mathematics may give some hints at the theoretical soundness of 
the concepts involved. 



7.1.6 The Combined Nonlinear— Condensation Transformation (CNCT) 

We assume the scries X^fcLo slowly convergent and nonalternating, i.e. Ofc > Vfc. In order to 

make the notation of involved subscripts more clear, we temporarily identify = a{k). Following Van 
Wijngaarden [219], we transform the nonalternating input series 

oo 

^a(fc), a(fc)>0, (7.32) 

fc=0 

whose partial sums are given by (7.7), into an alternating series X]^o(~l)"'-^j • After the first step of the 
transformation, the limit of the input series is recovered according to 

oo oc 

^a(fc)==^(-l)^A,. (7.33) 

fc=0 i=0 
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The quantities Aj are defined according to 



A,=Ebi'\ (7.34) 

where 

b[^^ = 2^= a(2'^ (j + 1) - 1) . (7.35) 

Obviously, the terms Aj defined in Eq. (7.34) arc all positive if the terms a{k) of the original series are 
all positive. The Aj arc referred to as the condensed series [23], and the series YlJ^Qi—^y Aj is referred 
to as the transformed alternating series, or alternatively as the Van Wijngaarden transformed series. 

The construction of the condensed series reminds one of Cauchy's condensation theorem (sec e.g. p. 28 
of Ref. [220] or p. 121 of Ref. [221]). Given a nonaltcrnating series X)fc°=o'^(^) '^ith terms that satisfy 
\a{k + 1)\ < \a{k)\, Cauchy's condensation theorem states that X)fcLo'^(^) converges if and only if the 
first condensed sum Aq defined according to Eq. (7.34) converges. 

The summation over k in Eq. (7.34) dose not pose numerical problems. Specifically, it can be easily 
shown in many cases of practical importance that the convergence of X^fc^o ^k"^ (^'^ linear even 

if the convergence of J2kLo '^(^) logarithmic. We will illustrate this statement by way of two 

examples. 

Example 1: a logarithmically convergent input series whose terms behave asymptotically as a{k) ^ k~^~^ 
with e > 0. In this case, the partial sums 



a(") = Y: bjf") (7.36) 



converge linearly with 



a("+i)_A- 1 
lim -^-j- = ——. -TTir<'^, a{k)r^k-^-\ k ^ oo. (7.37) 

Example 2: a series with a{k) ~ k^r'' where < ?' < 1 and /? real. Here, we have p = r < 1, and the 
series is (formally) linearly convergent. However, slow convergence may result if p is close to one. In this 
case, the condensed series are very rapidly convergent, 

a("+i) - A- 

lim -^-^ = 0, a(fc)~A:^r^ k-*oo. (7.38) 

Aj 

Therefore, when summing over k in evaluating the condensed series according to Eq. (7.34), it is in 
many cases suflScient to evaluate the condensed series by adding the terms successively, and no further 
acceleration of the convergence is required. 

As shown in [23, 222], the condensation transformation defined according to Eqs (7.33)- (7.35) is essen- 
tially a reordering of the terms of the input series X^fcLo o(fc). Furthermore, Daniel was able to show (see 
the Corollary on p. 92 of Ref. [222]) that for nonaltcrnating convergent series whose terms decrease in 
magnitude (|a(fc)] > \a{k + 1)]), the equality (7.33) holds. This formally justifies the correctness of the 
condensation transformation defined according to Eqs. (7.33) - (7.35). 

Note that the property 

A2,_i-i [A,_i-aj_i] , (j = l,2,...), (7.39) 

derived in [222], facilitates the numerical evaluation of a set of condensed series, by reducing the evaluation 
of condensed series of odd index to a trivial computation. 

In the second step of the CNCT, the convergence of the Van Wijngaarden transformed series 
YlY=o ("l)"* on the right-hand side of Eq. (7.33) is accelerated by a suitable nonlinear sequence trans- 
formation. We start from the partial sums 

n 

S« = E(-l)'Aj (7.40) 
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of the Van Wijngaarden transformed series. 

We use the partial sums S„ defined in Eq. (7.40) as input data for the delta transformation (7.31). This 
leads to the CNC transforms 

TcNc(«)=^i"'(l,So) . (7.41) 
which require as input the elements {So, . . . , S„, S„+i} of the Van Wijngaarden transformed series. 



7.2 Applications of Convergence Acceleration Methods 
7.2.1 Applications in Statistics and Applied Biophysics 

Several slowly convergent series of the type considered in this Section of the Thesis define special functions 
that have important uses in statistics, namely the Riemann zeta, generalized zeta, Jonquiere functions and 
Lerch's transcendent (sec [f)2] for definitions and some properties of these functions). Discrete distributions 
are related to these functions by associating probability mass functions (p.m.f.s) with the terms of the 
particular infinite series defining the special functions. In turn, the normalization constants of the p.m.f.s 
are associated with the sums of the series. 

Consider a discrete distribution with the p.m.f. 

Pr[X = n] = cp{n;e) (7.42) 

and with support of all nonnegative integers n = 0, 1, . . . , where is a (vector of) parameter(s) and c is 
a constant. The probability over the set of all outcomes is unity: 



E P^^t^ = ^] = 1 ' (7-43) 



n=0 

from which we obtain the normalization constant 



-^ = J2p{n;e) (7.44) 



n=0 

(in statistics, the random variable is usually denoted by the symbol n, whether continuous or discrete, 
but we will denote it here by the symbol n in order to follow the usual notation used in the treatment of 
special functions). Thus we expect the sum in Eq. (7.44) to exist and to converge to a finite value. In the 
sequel, we consider several related discrete distributions. The Zipf distribution has the p.m.f. 

Pcc=c^ (7.45) 

with support of all positive integers n — 1,2,.... The normalization constant (7.44) is then 

oo 
n=l 

where s > 1, which is immediately recognized as the Dirichlet series for the Riemann zeta function [Eq. (1) 
on p. 32 of Ref. ['■>-]]■ Basic properties of the Zipf distribution can be found in [22.3]. 

The Zipf-Mandelbrot distribution [224] is the generalization of the Zipf distribution that adds a constant 
to the ranks n and has the p.m.f. 

P. = c— ^ (7.47) 
with the support of all nonnegative integers n = 0, 1, . . . . The normalization constant (7.44) is then 



n=0 
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where w 7^ 0, —1, .... The sum (7.48) is recognized as the generahzed zeta function [Eq. (1) on p. 24 of 
Ref. [')2]]. 

Fmther, the p.ni.f. 

Px = c^ (7.49) 

with the support of all positive integers n = 1,2,... defines the Good distribution ['^?-^-]. The con- 
stant (7.44) is then 

00 ^ 



where \z\ < 1, which is recognized as the Jonquiere's function [Eq. (14) on p. 30 of Ref. [!)2]; see also [22G]]. 
Important properties of the Good distribution were derived in [227, 228]. 

In biophysics, one is interested in establishing the statistical structure of DNA and RNA sequences. 
Indeed, an interest in statistical compositions of DNA and protein sequences originated shortly after 
the discovery that the nucleotide sequences of RNA uniquely determines the amino-acid sequence of the 
proteins. An early analysis was performed in [22!)], where it was observed that the distributions of relative 
abundances of amino-acids and nucleotides were nonrandom, i.e. deviated from models assuming uniform 
distribution of different words in any given text. Zipf distributions and related generalizations proved to 
be very useful in further studies [230, 231, 232]. 

A formal generalization of the Zipf, the Zipf-Mandelbrot and the Good distributions is possible upon 
realizing that the Riemann zeta, the generalized zeta and the Jonquiere functions constitute special cases 
of the Lerch transcendent which is defined by the following series 

00 „ 

^z,s,v) = J2-^-—, (7.51) 

where |z| < 1 and v ^ 0, —1, . . . [Eq. (1) on p. 27 of Ref. [92]]. We use here the the nonstandard notation 
of a dummy index n in order to highlight the connection with statistical theory. The relations between 
<i>(z, s, v) and C(s)j C(sj ''^) ^^nd F{z, s) can be obtained by considering appropriate limits if $(z, s, v) (e.g., 
by letting z and/or u to and/or 1), and by making use of the following functional relation [Eq. (2) on 
p. 27 of Ref. [92]] 

™-l n 

$(z,s,«) = z"$(z,s,m + i;)+ V -— — . (7.52) 

,t^o (n + vy 

In particular, Eq. (7.52) with m ~ I can be used to evaluate sums that start from n = 1 rather than 
from n — 0: 

00 „ 

J2^^-—^zHz,s,l + v) (7.53) 

n— 1 ^ ^ 

In Tab. 7.1, we have summarized the relations between the different statistical distributions and their 
defining special functions. The table illustrates the fact that the three distributions, Zipf (Zeta), Zipf- 
Mandelbrot (Generalized zeta) and Good (Jonquiere), are special cases of the Lerch distribution, and 
that their properties can be expressed in terms of Lcrch's transcendent with special values of parameters. 

We have found that the application of the combined nonlinear-condensation transformation (CNCT) 
described in detail in Sec. 7.1.6 is very suitable for the accelerated evaluation of the Lerch transcendent 
(7.51) in the nonalternating case (z > 0). Details of the implementation, both in C and Mathematica [90], 
are discussed in [30, 31]. For z < 0, the series (7.51) is alternating. In this case, the direct application 
of the delta transformation (7.31) leads to an efficient calculational scheme [31]. A C code for the Lerch 
transcendent is available [233, 234]. 

7.2.2 An Application in Experimental Mathematics 

We would like to begin this section by quoting [235]: "In April 1993, Enrico Au-Yeung, an undergraduate 
at the University of Waterloo, brought to the attention of the author's [David Bailey's] colleague Jonathan 
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Distribution Support p.m.f. c 

(related function) 



Zipf 

(zeta function) 

Zipf-Mandelbrot 
(generalized zeta) 

Good 

(Jonquiere function) 
Lerch 

(Lercli's transcendent) 



1,2,... cn-'* 

0,1,... c{n + vy 

1,2,... cz^n"" 

0,1,... cz"(n + w)-" 



C(s)= $(1,5,1) 

F{z,s) = z$(z,s,l) 
<^{z,s,v) 



Table 7.1: Interrelationships between distributions defined by Riemann zeta, generalized 
zeta, Jonquiere and Lerch functions. 



Borwein the curious fact that 

|(l4 + - + i)^--«<.987-...f«4) = lg!, 

based on a computation of 500,000 terms. Borwein's reaction was to compute the value of this constant 
to a higher level of precision in order to dispel this conjecture. Surprisingly, his computation to 30 digits 
affirmed it. [David Bailey] then computed this constant to 100 decimal digits, and the above equality was 
still affirmed." Many formulas similar to (7.54) have subsequently been established by rigorous proof [236]. 

With the help of a multiprecision system [235, 237, 238] and the CNCT, we have verified (7.54) "experi- 
mentally" to a couple of hundred decimals. The calculation will be sketched in the following. Using the 
definition 

m 

we rewrite (7.54) as follows, 

oc 
fc=0 



-5 



I' 



fe=o \j=o / 



(7.55) 



fc=0 



where C is the Euler-Mascheroni constant C ~ 0.577 215 . . ., and ^^{z) is the logarithmic derivative of 
the Gamma function [91, 92] , 

^(z) = ^lnr(z). (7.57) 
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With the help of the relation 



k=0 ^ ' 



the conjecture (7.54) can be rewritten as 



E i^^TTT-) ' - T ^ W - 4C C(3) + C(2) . (7.59) 



V fc + 1 / 4 

fc=0 ^ ' 

We proceed to calculate numerically, to high precision, the infinite sum 



^a(fc), (7.60) 



fc=0 



using the CNC transformation. 

In order to establish the rate of convergence of (7.60), we investigate the asymptotic behaviour of the 
a(k) as k oo. The logarithm of the Gamma function can be expanded into an asymptotic series [see 
Eq. (4.03) on p. 294 of 

Inr(z) = i^z--j\nz - z + - ln(2^) + 2.(2. -'l)z^-^ + ^^'^^^ 

where 

Here, Jx]] is the integral part of x, i.e., the largest integer m satisfying m < x, Bk{x) is a Bernoulli 
polynomial defined by the generating function [see Eq. (1.06) on p. 281 of Ref. [91]]: 



^ = Brnix)—:, \t\ < 2n , (7.63) 
1 — 1 ml 

and 



exp(i) — ^ 

' m=0 



B,n = S™(0) (7.64) 
is a Bernoulli number (p. 281 of Ref. [91]). The following asymptotic relation for il]{z) follows: 

*W.ln,-l_'f ||^+o(_L)^ (7,65) 

The leading asymptotics of the remainder of the sum (7.60) after adding A^ — 1 terms can thus be derived 
easily. We have for large A:, 

ln(A: + 2)2 ln(/c + 2) ln(fc + 2) ( ^ \ , 

^ lfcTTF"(fc + l)Mfc + 2)"6(fc + l)Mfc + 2)2+^lFj ' ^^-^^^ 

Based on these formulas, the remainder of the sum (7.60), for large A^, can be written as 

T.m-^ + ^ + j^+0^^y (7.67) 

k=N ^ ^ 

Here, the Euler-Maclaurin formula [Eqs. (2.01) and (2.02) on p. 285 of Ref. [!)]]] has been used in order 
to convert the sum over the a{k) in the asymptotic regime of large k [see Eq. (7.66)] into an integral 
plus correction terms. In order to calculate (7.59) to an accuracy of 200 decimals, Eq. (7.67) says that 
we would be required to add on the order of 10^"^ terms. Without the use of convergence acceleration 
methods, this would represent a formidable computational task. 
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Using the CNCT, it is easy to calculate the sum (7.59) to 200 digits, based on multiprecision arith- 
metic [235] and a Linux personal computer, within a few hours. We obtain for the 246th and the 247th 
CNC transform defined according to Eq. (7.41), 

rcNc(246) = 2.37254 51620 38445 67035 68130 69148 85258 25756 18499 54254 
97013 57806 20011 72404 62937 46020 32218 23862 67095 00004 
69194 36541 28946 10390 15116 52595 90270 23975 58737 74256 
23420 48480 95165 00802 19816 35378 76591 98589 60393 32103, (7.68) 

and 

rcNc(247) ^ 2.37254 51620 38445 67035 68130 69148 85258 25756 18499 54254 
97013 57806 20011 72404 62937 46020 32218 23862 67095 00004 
69194 36541 28946 10390 15116 52595 90270 23975 58737 74256 

23420 48480 95165 00802 19816 35378 76591 98589 60393 32112. (7.69) 

The apparent convergence to 200 decimals can be verified against the right-hand side of Eq. (7.59). Of 
course, the right-hand side of Eq. (7.59), which involves only rationals, zeta functions and the Euler- 
Mascheroni constant, 

^C(4)-4CC(3)-f C2C(2), 

can be easily evaluated to 200 decimals using known algorithms which are included in computer algebra 
systems [90]. 

The evaluation of the terms a{k) proceeds as follows. For small index k, it is easy to write a recursion 
relation relating a{k) and a(fc -I- 1) based on the (trivial) recursion for the ip function, 

^{k + 1) = iik) + y . (7.70) 
k 

For large index k, the asymptotic formula (7.65) can be used in order to calculate the ip function to high 
precision. The point at which index one may switch from the recursion to the asymptotic method depends 
on how many explicit values for Bernoulli numbers are available to the machine. We use the first values 
for the first 60 Bernoulli numbers, to 250 decimals, for our calculation. We switch from one method to 
the other when the index k of a{k) has reached a value of 500. 

With 84 308 a{k) terms evaluated (out of which 1364 by recursion and 82944 by the asymptotic method), 
we evaluate the first 247 transforms with the results presented above in Eqs. (7.68) and (7.69). If the 
terms of the series (7.54) were added on a term-by-term basis, then about 10^°^ would be required for an 
accuracy of 200 decimals in the final result. The reduction of this number to roughly 84 000 corresponds 
to an acceleration of the convergence by roughly 200 orders of magnitude. We have also carried out, 
using enhanced precision, a calculation to 430 decimals, involving about 500 CNC transformations and 
arithmetic with 600 decimal figures. These evaluations not only confirm the conjecture (7.59) to high 
precision, but they also represent one of the most accurate experimental verifications of the convergence 
properties of the delta transformation (7.31) today. Moreover, it is observed that the rate of convergence 
of the CNC transforms, which is approximately one decimal figure per transformation, remains constant 
over a wide range of transformation orders. This observed behaviour is consistent with linear convergence 
[p « 0.1 in Eq. (7.6)] in higher transformation orders. By contrast, the series (7.59) is only logarithmically 
convergent. This corresponds to convergence acceleration according to the definition (7.9). 
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7.2.3 Other Applications of the CNCT 



We would like to mention the existing applications in the domain of quantum electrodynamic bound-state 
calculations (see [9, 12, 43] and Ch. 1 of this Thesis). Another existing application concerns the quantum 
electrodynamic effective action (see Ref. [2U]). 

J. Becher [239] has investigated the applicability of the combined nonlinear-condensation transformation 
to series of the form 

°° ^2fe+l 

RJx) = y — ^ — , (7.71) 

^ ^ 2fc + 1 p ^ ^ 



fe=0 



as well as 



and 



oo 



rr. , ,X X ■ 1 C0Sll(2fc+l)x 



Up(x,h) = V : ,vn^-T);7 : (7-73) 



{2k + l)P cosh(2fc + 1)5 ' 



1 cosh(2A: -I- l)x 



k=0 



{2k + I)?' sinh(2fc + 1)6 



Series of this type occur naturally in the context of plate contact problems with periodic boundary 
conditions [240]. The arguments p, x and 6 are real and positive for cases of practical relevance. For x ^ b 
and p « 1, the series Tp and Up are very slowly convergent. In App. A. 2 of [239], it is demonstrated by 
way of numerical example that the CNCT is able to efficiently accelerate the convergence of these series 
in problematic parameter regions. 

In the numerical calculations, it is necessary to evaluate terms with large index k. This can lead to 
numerical overflow because of the large arguments of the hyperbolic functions. Clearly, representations 
like 

g(2fc+l)(i;-6) _(_ g-2x(2fe+l) 

(2A: + 1)P 1 + e-26(2fc+i) 



Tpi^^b) = -l^^TTT^ . , :-2.m^n (7-74) 



A:=0 

provide a solution for this problem. 

Let us recall that considerable effort has been invested in the development of efficient numerical methods 
for the evaluation of the series (7.71) - (7.73) [241, 242, 243, 244, 24-5]. These alternative methods make 
intensive use of special properties of the series. They involve integral transformations and infinite series 
over numerical integral [244], and they make use of special properties of Legendre's chi-function [243] 
which is related to the functions (7.71) - (7.73). 

A further application of convergence acceleration methods discussed by J. Becher [239] is the calculation 
of the incomplete Gamma function r(0,x) whose asymptotic expansion for large argument reads 



r(0,^) - ^ E("ir^' a;-.oo.. (7.75) 



It is observed (see App. A. 3 ibid.) that the divergent asymptotic series of this function for large positive 
argument can be resummcd effectively using the delta transformation, without any recourse to the CNCT. 



7.3 Conclusions and Outlook for Convergence Acceleration 

We have discussed several applications of the convergence acceleration methods introduced in Sec. 7.1: in 
statistical physics (Sec. 7.2.1), in experimental mathematics (Sec. 7.2.2), and other applications, mainly 
in the evaluation of special functions (Sec. 7.2.3). Specifically, it is observed that a combination of the 
methods introduced in Sees. 7.1.5 and 7.1.6 leads to an efficient calculational scheme for the Lerch 
transcendent (f> defined in Eq. (7.51). This special function provides a generalization of several kinds of 
probability density functions which are of significance for the statistical analysis of DNA sequences (see 
Table 7.1). 

Other applications can be envisaged: For example, the evaluation of generalized hypergeometric functions 
~ notably of two-variable hypergeometric functions (Appell functions) Fi and F2 [24G] - has been histor- 
ically problematic. In the recent investigation [247], the evaluation of i^i(a,/3i,/32,7;a;,?/) is considered: 
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The algorithm consists in two steps: (i) the sophisticated use of analytic continuations to map the point 
(x, y) into the convergence region of the defining hypergeometric series of the Appell function, and (ii) 
the computation of the Appell function within the region of convergence. The analytic continuations give 
rise - in addition to Fi functions of transformed argument ~ to the G2 function, which can in turn be 
written in terms of Appell's F2 function. This implies that numerical methods also have to be devised for 
the F2 function. The second step of the algorithm described in [247] is the actual numerical calculation of 
the Appell function(s) (Fi and F2) within the region of convergence. This calculation proceeds by inves- 
tigating the criterion to = max(|a::|, \y\) < 0.5. If the criterion is fulfilled, then both the Fi and the F2 are 
evaluated by single- index series expansions in terms of 2^"i's, for example employing the formula [248, 249] 



If the criterion is not fulfilled, the authors resort to the solution of a third-order differential equation 
for the evaluation of the Appell function [starting from the point {x,y) = (0,0) where Fi vanishes]. It 
appears possible that an application of convergence acceleration algorithms to the series (7.76) could lead 
to an efficient calculational scheme in the problematic case to > 0.5 where the authors of [247] had to 
resort to an integration of a third-order differential equation. 

Possible improvements over existing algorithms concern "multi-stage transformations" . For example, it 
was observed in [20] that the rate of convergence of the output data resulting from Borel summation 
(see Sec. 8.2.2) of a specific divergent series describing the quantum electrodynamic effective action could 
be further accelerated by using two iterations of the epsilon algorithm [defined in Eq. 7.14)], i.e. by 
evaluating Pade approximants to the Pade approximants of the Borel transforms (in this sequence). 
Iterated forms of a number of convergence accelerators have been discussed in [156]. In general, the 
iteration of suitable convergence accelerators depends on a thorough analysis of the remainder structure 
resulting from previous "stages" of the transformation, and it is not a priori guaranteed that a "multi-stage 
transformation" with only one iterated algorithm will help. In principle, the use of iterated algorithms can 
only be helpful when the evaluation of the sequence transforms or Pade approximants is computationally 
cheap in comparison to the evaluation of further terms of the input sequence which could potentially be 
used to calculate higher-order terms in a single-stage transformation process. 




X 2Fi{a + r, I3i + r,j + 2r,x) 2Fi{a + r, (32 + r^'y + 2r,y) . 



(7.76) 
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Chapter 8 



Divergent Series 



8.1 Introduction to Divergent Series in Physics 

The observation that many perturbation series encountered in quantum mechanics and quantum field 
theory are divergent, sparked the development of large-order perturbation theory [-50]. This subject is 
inherently connected to the resummation of divergent series, and with the development of appropriate 
numerical methods which exploit all available information contained in a necessarily finite number of 
perturbative coefficients. Of prime importance in this context is the Borel method to be discussed in 
Sec. 8.2.2. 

As already stressed in Ch. 8, there are connections between the convergence acceleration methods dis- 
cussed in Ch. 7 and the resummation methods discussed in the current Chapter. Specifically, the methods 
introduced in Sees. 7.1.4 and 7.1.5 will be used below in Sec. 8.3.2 in the context of the resummation of 
the divergent perturbation series for the QED effective action. 

Essentially, we will consider four physical applications of resummation methods: (i) the energy displace- 
ment of a hydrogen atom in a background electric field, including the autoionization width (Sec. 8.2), 
(ii) zero-dimensional field theories (Sec. 8.3.1), (iii) the QED effective Lagrangian for background mag- 
netic and electric fields (Sec. 8.3.2), and (iv) the energy levels of the double-well potential which receive 
nonperturbative contributions from multi-instanton effects (Sec. 8.3.3). 

Several important aspects and basic concepts connected with divergent series and the large-order behavior 
of perturbation theory can be illustrated by investigating the (divergent) perturbation series of the Stark 
effect [251]. This effect, which is also known as the LoSurdo-Stark effect in view of the existence of the 
investigation [252]), describes the energy displacement of a hydrogen atom in a static external electric 
field and can be expressed as a perturbation series (power series) in the electric field strength. 

The first nonvanishing perturbation (in atomic units) is the second-order effect 



where the sum over m runs over the entire spectrum, including the continuum, but excluding the IS ground 
state, and Em is the nonrelativistic (Schrodinger) energy of the mth state. A well known, but perhaps 
surprising result says that the coefficients of the terms of order F^, F^, F^, . . . grow so rapidly that the 
series in F ultimately diverges, irrespective of how small the field strength is. The convergence radius of 
the factorially divergent perturbation series is zero. The perturbation series for the Stark effect [251, 252] 
can be formulated to arbitrarily high order [253] . The perturbative coefficients grow factorially in absolute 
magnitude [254]. The perturbation series is a divergent, asymptotic expansion in the electric field strength 
F, i.e. about zero electric field. This means that the perturbative terms at small coupling first decrease in 
absolute magnitude up to some minimal term. After passing through the minimal term, the perturbative 
terms increase again in magnitude, and the series ultimately diverges. 

By the use of a resummation method, it is possible to assign a finite value to an otherwise divergent 
series, and various applications of resummation methods in mathematics and physics have been given. 
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e.g., in [24, 156, 165, 179, 255]. When a divergent series is resummcd (cf. Ch. 6), the superficial precision 
hmit set by the minimal term can be overcome, and more accurate results can be obtained as compared 
to the optimal truncation of the perturbation series. However, the resummation of the relevant divergent 
series is problematic in the considered case, because the Borel transform, from which the physically 
correct, finite result is obtained by evaluating a generalized Borcl integral ~ see Eq. (8.11) in Sec. 8.2.2 
below -, exhibits a rich singularity structure in the complex plane. The perturbation series is formally 
not Borel summable. 

The Stark efi^ect and its associated divergent perturbative expansion including the nonperturbative, non- 
analytic imaginary contributions have attracted considerable attention, both theoretically and experi- 
mentally [157, 253, 254, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, 
273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 283, 284, 285, 286]. Experiments have been performed in 
field strengths up to a couple of MV/cm [287, 288, 289, 290]. 

Rather mathematically motivated investigations regarding the Borel summability of the divergent per- 
turbation series for the Stark effect were performed in [265] and [159], and it was established that the 
perturbation series of the Stark effect is Borel summable in the distributional sense (for the definition of 
"distributional Borel summability" we refer to [158]). 

In comparison to the previous investigation [157], we use here an optimized resummation method that 
combines an analytic continuation of the Borel transform with Fade approximants in the conformal 
variable and the explicit use of the "leading renormalon pole". All of these methods will be described 
in Sees. 8.2.2 and 8.2.3. These optimized Borel transforms represent an alternative to so-called order- 
dependent mappings which have been discussed intensively in the literature [291, 292, 293, 294]. 

After the discussion of the Stark effect (Sees. 8.2.2 - 8.2.4), we turn our attention to the resummation 
of the divergent perturbation series for zero-dimensional field theories (Sec. 8.3.1), to the QED effective 
action (Sec. 8.3.2), and to the double-well potential (Sec. 8.3.3). 



8.2 The Stark Effect: A Paradigmatic Example 
8.2.1 Perturbation Series for the Stark Effect 

In the presence of an electric field, the SO (4) symmetry of the hydrogen atom is broken, and parabolic 
quantum numbers ni, 712 and m are used for the classification of the atomic states [KiS]. For the Stark 
effect, the perturbative expansion of the energy eigenvalue E{ni,n2,m, F) reads [see Eq. (59) of [253]], 

00 

E{n,,n2, m,F)^Y. ^n%m , (8-1) 

where F is the electric field strength. For N ^ 00, the leading large-order factorial asymptotics of the 
perturbative coefficients have been derived in [269] as 

where A[!^Ji-m is given as an asymptotic series, 

00 

A^^Xm-K{n,,n„m,N) ^ a^'"^™ (2 + 771 + iV - /c)! . (8.3) 

The quantities q^'^j"* ^re constants. The if-coefficients in Eq. (8.3) arc given by 

K(n,,nj,m,N) = ~ [27rn3nj! (n^- + m)!] exp {3 (71, - n^)} 6^"^+™+! (3nV2)^ . (8.4) 

Here, the principal quantum number n as a function of the parabolic quantum numbers ni, 712 and m is 
given by [see Eq. (65) in [253]] 

77, = 77,1 + 772 -I- |?77| + 1 . (8.5) 

According to Eq. (8.2), the perturbative coefficients E^n^m^ for large order iV — > oo of perturbation 
theory, can be written as a sum of a nonaltcrnating factorially divergent series [first term in Eq. (8.2)] 
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and of an alternating factorially divergent series [second term in Eq. (8.2)]. Because the i 



inEq. (8.3) 



are multiplied by the factorial (2ni + m + N — k)l, we infer that for large perturbation theory order TV, 
the term related to the a^' ^ coefRcient (k = 0) dominates. Terms with fc > 1 are suppressed in relation 
to the leading term by a relative factor of XjN^ according to the asymptotics 

{2nj + m + /V - A;)! 



1 



l + O 



(2nj + m + iV)! 

for — > oo. The leading (fc — 0) -co efficient has been evaluated in [2.j4] as 



an = 1 ■ 



(8.6) 



(8.7) 



According to Eqs. (8.2), (8.3) and (8.7), for states with ni < 712, the nonalternating component of the 
perturbation series dominates in large order of perturbation theory, whereas for states with ni > 712, 
the alternating component is dominant as iV ^ 00. For states with ni — 712, the odd-A^ pcrturbative 
coefficients vanish [2(1' l], and the even-A^ coefficients necessarily have the same sign in large order [see 
Eq. (8.2)]. According to Eq. (8.2), there are subleading divergent nonalternating contributions for states 
with ni > 77,2, and there exist subleading divergent alternating contributions for states with ni < 712. 
This complicates the resummation of the perturbation series. 
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Figure 8.f: Integration contour C+i for the evaluation of the generalized Borel integral 
defined in Eq. (8.f f). Poles displaced from the real axis are evaluated as full poles, whereas 
those poles which lie on the real axis are treated as half poles. 



8.2.2 Borel— Pade Resummation 

The resummation of the perturbation series (8.f ) by a combination of the Borel transformation and Fade 
approximants proceeds as follows. First we define the parameter 

A = 2 max(77i, 7t2) + 7n, + f . (8.8) 
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The large-order growth of the perturbative coefficients [see Eqs. (8.2) and (8.3)] suggests the definition 
of the (generahzed) Borel transform [see Eq. (4) in [2()]] 



Eb{z) = £'b(?t.i, ^2, m, z) 



E 



(TV) 

■nin2m ]\[ 



(8.9) 



where we consider the argument z of E-b{z) as a complex variable and A is defined in Eq. (8.8). The 
additive constant (in this case A) in the argument of the Gamma function is chosen in accordance with 
the notion of an "asymptotically improved" resummation (see also [2'<]). It is observed that the additive 
constant A can be shifted by a small integer without affecting the convergence of the Borel rcsummed 
series. Because the perturbative coefficients E^n^m diverge factorially in absolute magnitude, the Borel 
transform E^^z) has a finite radius of convergence about the origin. The evaluation of the (generalized) 
Laplace-Borel integral [see Eq. (8.11) below] therefore requires an analytic continuation of E-giz) beyond 
the radius of convergence. The "classical" Borel integral is performed in the z-range z S (0, oo), i.e. along 
the positive real axis [see e.g. Eqs. (8.2.3) and (8.2.4) of [ICo]]. It has been suggested in [295] that the 
analytic continuation of (8.9) into regions where F retains a nonvanishing, albeit infinitesimal, imaginary 
part can be achieved by evaluating Fade approximants. Using the first M+l terms in the power expansion 
of the Borel transform E-b{z), we construct the Fade approximant [we follow the notation of [179], see 
also Eq. (7.12)] 



Vm{z) = 



p.//2l/p/ + l)/2l 



(z) , (8.10) 



where \x\ denotes the largest positive integer smaller than x. We then evaluate the (modified) Borel 
integral along the integration contour C+i shown in Fig. 8.1 in order to construct the transform TEm{F) 
where 

TEMiF) = [ dtt^-' eM-t)VMiFt). (8.11) 

The successive evaluation of transforms TE]\i{F) in increasing transformation order M is performed, and 
the apparent convergence of the transforms is examined. This procedure is illustrated in Tables I and II 
of [2.5]. 

The contour C+i is supposed to encircle all poles at t — Zi in the upper right quadrant of the complex 
plane with arg z^ < 7r/4 in the mathematically negative sense. That is to say, the contribution of all poles 
Zi with Rezi > 0, Imz; > and Rez,; > Imz^, 



-2 7ri V Rest^"^ exp(-t)pM(Ft) , 

^ ^ t — Zi 

i 

is added to the principal value (P.V.) of the integral (8.11) carried out in the range t € (0, oo). We have, 

TEMiF) = (F.V.) / dtt^-^ expi-t)rMiFt) 
Jo 

-27ri V Rest^^i exp(-t) ^^/(i^t) ■ (8.12) 

^ — ^ t—Zi 
i 

The principal- value prescription [first term in Eq. (8.12)] for the evaluation of the Laplace-Borel integral 
has been recommended in [295, 296]. This prescription leads to a real (rather than complex) result 
for the energy shift and cannot account for the width of the quasistationary state. The additional pole 
contributions [second term in Eq. (8.12)] are responsible for complex- valued (imaginary) corrections which 
lead, in particular, to the decay width. 
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By contour integration (Cauchy Theorem) and Jordan's Lemma, one can show that the resuh obtained 
along C-fi is equivalent to an integration along the straight line with argz = 7r/4, 

TEm{F)^c^ I dtt^-^ cxp{-ct)VM{Fct) (8.13) 
Jo 

where c — exp(i7r/4). This contour has been used in [157] (see also p. 815 in [2')7]). The exponential 
factor exp(— ci) and the asymptotic behavior of the Pade approximant Vm{FcA) as i ^ go together 
ensure that the integrand falls off sufficiently rapidly so that the Cauchy Theorem and Jordan's Lemma 
can be applied to show the equivalence of the representations (8.12) and (8.13). 

The representation (8.13) illustrates the fact that the integration in the complex plane along C+i analyt- 
ically continues the resummed result in those cases where the evaluation of the standard Laplace-Borel 
integral is not feasible due to poles on the real axis. The representations (8.11) and (8.12) serve to clarify 
the role of the additional terms which have to be added to the result obtained by the principal-value 
prescription in order to obtain the full physical result, including the nonperturbative, nonanalytic contri- 
butions. Note that, as stressed in [25], the pole contributions in general do not only modify the imaginary, 
but also the real part of the resummed value for the perturbation series. 

Formally, the limit of the sequence of the TEm{F) as M s- oo, provided it exists, yields the nonpertur- 
bative result inferred from the perturbative expansion (8.1), 

lim TEm(F) = E(F) = E(ni,n2,m,F) . (8.14) 

M^oo 

Because the contour C+i shown in Fig. 8.1 extends into the complex plane, the transforms TEm{F) 
acquire an imaginary part even though the perturbative coefficients in Eq. (8.1) are real. 

In the context of numerical analysis, the concept of incredulity [298] may be used for the analysis of the 
convergence of the transforms TEm{F) of increasing order M. If a certain number of subsequent trans- 
forms exhibit apparent numerical convergence within a specified relative accuracy, then the calculation 
of transforms is stopped, and the result of the last calculated transformation is taken as the numerical 
limit of the series under investigation. It has been observed in [25, 295] that for a number of physically 
relevant perturbation series, the apparent numerical convergence of the transforms (8.11), with increasing 
transformation order, leads to the physically correct results. 

It is observed that the rate of convergence of the transforms (8.11) can be enhanced if instead of the 
unmodified Pade approximants (8.10) leading "renormalon" poles are explicit used for the construction 
of modified approximants. For the ground state, this entails the following replacement in Eq. (8.11): 

Vm{z)^V'm{z), 

where 



V'm{z) ^ 



1 



,M-t-4, /,M-3, 



(z) , (8.15) 



where £'b(C) = (1 ^ C^) £'b(C)- For the excited state with quantum numbers rti = 3, 712 = and m = 1, 
we replace 

Vm{z)^V"m{z), 

where 



V"m{z) = ^ 



1 



,M + 2^ LM - 1, 



(z) , (8.16) 

E'' 



where £'b(C) = (1 ~ C^) -E'b(C)- The resummation method by a combination of Borel and Pade techniques 
- current Section ~ will be referred to as "method I" in the sequel. 



8.2.3 Doubly-Cut Borel Plane 

According to Eqs. (8.2) and (8.3), the perturbative coefficient En^n^m, for large N, can be written as 
the sum of an alternating and of a nonalternating divergent series. In view of Eqs. (8.4) and (8.7), we 
conclude that the series defined in Eq. (8.9), 

(TV) 



Eb[z) = > — : 7Z . 

^„r(7V + A) 
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has a radius of convergence 

about the origin, where n is the principal quantum number [see Eq. (8.5)]. Therefore, the function 

£bH = 1^ r(iv + A) ^ ' ^^-^^^ 

has a unit radius of convergence about the origin. It is not a priori obvious if the points w ~ —I and 
w = +1 represent isolated singularities or branch points. The asymptotic properties (8.2) and (8.3) 
together with Eq. (8.6) suggest that the points w = —I and w = +1 do not constitute poles of finite 
order. We observe that the leading factorial growth of the perturbative coefficients in large perturbation 
order N is divided out in the Borel transform (8.18), which is a sum over N. The perturbative coefficient 
En^n^m Can be written as an asymptotic series over k [see Eq. (8.3)]. We interchange the order of the 
summations over N and fc, we use Eq. (8.6) and take advantage of the identity 

E^=Li.H. (8.19) 

Af=0 

The Borel transform £b{w) can then be written as a sum over terms of the form Tk{w) where for k — > oo, 

Tk{w) ^ Cin„ n„m) a^^'"'"" Likiw) . (8.20) 
The coefficient C {ui , rij , m) is given by 

C{n„nj,m) ^ - [27rn3nj! (n^- + m)!] cxp {3 (n, - n^)} 62"^+™+i. (8.21) 

These considerations suggest that the points w ~ — \ and u; = +1 represent essential singularities (in this 
case, branch points) of the Borel transform £b.{w) defined in Eq. (8.18). For the analytic continuation of 
£b(w) by conformal mapping, we write w as 

w ^ (8.22) 
1 + y2 

(this conformal mapping preserves the origin of the complex plane). Here, we refer to w as the Borel 
variable, and we call y the conformal variable. We then express the Mth partial sum of (8.18) as 



iV=0 ^ ' 



M( \ \ ^ J^nin^m ^ ]\[ 



M 



N=0 



where the coefficients Cn are uniquely determined [see, e.g., Eqs. (36) and (37) of [299]]. We define the 
partial sum of the Borel transform, re-expanded in terms of the conformal variable y, as 



M 



JV=0 



We then evaluate (lower-diagonal) Fade approximants to the function S'^ (y) 



[M/21 / KM + l)/2] 



(y) ■ (8.25) 

c/M 



We define the following transforms. 



T'EuiF) = s^ dw w^-^ cxp{~w) £"b (yiw)) . (8.26) 
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At increasing M , the limit as M —t oo, provided it exists, is then again assumed to represent the complete, 
physically relevant solution, 

E{F)= lim T"Em{F). (8.27) 

We do not consider the question of the existence of this limit here (for an outline of questions related to 
these issues we refer to [■:)()()]; potential problems at excessively strong coupling are discussed in Sec. IIC 
of [301]). 

Inverting Eq. (8.22) yields [see Eq. (8.26)] 



VI + - \/l - w 
^/l + IV + \/l — w 

The conformal mapping given by Eqs. (8.22) and (8.28) maps the doubly cut w-plane with cuts running 
from u" = 1 to w = oo and w = —1 to w — —oo unto the unit circle in the complex y-plane (i.e., it is a 
surjective mapping). The cuts themselves are mapped to the edge of the unit circle in the y-plane. 

In comparison to the investigations [299] and [300] , we use here a different conformal mapping defined in 
Eqs. (8.22) and (8.28) which reflects the different singularity structure in the complex plane [cf. Eq. (27) 
in [299]]. We also mention the application of Pade approximants for the numerical improvement of the 
conformal mapping performed according to Eq. (8.25). In comparison to [27], where the additional Pade- 
improvement in the conformal variable is also used, we perform here the analytic continuation by a 
mapping whose structure reflects the double cuts suggested by the asymptotic properties of the pertur- 
bative coefficients given in Eqs. (8.2), (8.3) and (8.6) [cf. Eq. (5) in [27]]. 

The method introduced in this Section will be referred to as "method 11". One goal of the current 
investigation is to contrast and compare the two methods I and II. A comparison of different approaches 
to the rcsummation problem for series with both alternating and nonaltcrnating divergent components 
appears useful, in part because the conformal mapping (without further Pade improvement) has been 
recommended for the resummation of quantum chromodynamic perturbation series [299, 300]. 



8.2.4 Numerical Calculations for the Stark Effect 

In this section, the numerical results based on the resummation methods introduced in Sees. 8.2.2 and 8.2.3 
are presented. Before we describe the calculation in detail, we should note that rclativistic corrections 
to both the real and the imaginary part of the energy contribute at a relative order of (Za)"^ compared 
to the leading nonrclativistic effect which is treated here (and in the previous work on the subject, see 
e.g. [157, 269]). Therefore, the theoretical uncertainty due to rclativistic effects can be estimated to be, at 
best, 1 part in 10^ (for an outline of the rclativistic and quantum electrodynamic corrections in hydrogen 
see [1, 2, 48, 49, 50, 51, 302]). Measurements in very high fields are difficult [287]. At the achievable 
field strengths to date (less than 0.001 a.u. or about 5MV/cm), the accuracy of the theoretical prediction 
exceeds the experimental precision, and rclativistic effects do not need to be taken into account. 

The perturbative coefficients En^n^m. defined in Eq. (8.1) for the energy shift can be inferred, to arbitrarily 
high order, from the Eqs. (9), (13-15), (28-33), (59-67) and (73) in [253]. The atomic unit system is used 
in the sequel, as is customary for this type of calculation [253, 259, 262, 264]. The unit of energy is 
rrie (? = 27.211 eV where a is the fine structure constant, and the unit of the electric field is the field 
strength felt by an electron at a distance of one Bohr radius OBohr to a nucleus of elementary charge, 
which is 1/(4 TTEo) (e/aiohr) ~ 5.142 x lO^MV/cm (here, eg is the permittivity of the vacuum). 
We consider the resummation of the divergent perturbative expansion (8.1) for two states of atomic 
hydrogen. These are the ground state (rii = 712 = m = 0) and an excited state with parabolic quantum 
numbers rii = 3, 712 = 0, m = 1. We list here the first few perturbative coefficients for the states under 
investigation. For the ground state, we have (in atomic units), 

F (F\- ^ 9 2 3 555 4 2512779 e 13 012 777 803 . 
E,o,{F)-----F -^^F 16 384 ^ ^^'^^^ 

The perturbation series for the state ni = 3, 712 = 0, m = 1 is alternating, but has a subleading 
nonaltcrnating component [see Eq. (8.2)]. The first perturbative terms read 

F (F\- 31875 2 I 54140 625 3 715 751953125 , 

i530i(F)--- + yF-^-F + F F +... (8.30) 
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Note that for F = 0, the unperturbed nonrelativistic energy is recovered, which is — l/(2n^) in 
atomic units. In contrast to the real perturbative coefficients, the energy pseudoeigenvahie (resonance) 
E{ni,n2,m, F) has a real and an imaginary component, 

E{ni,n2,m,F) = Rc En^no.miF) - ■^Tmn27niF) , (8.31) 
where Tn-^n2m{F) is the autoionization width. 

Using a computer algebra system [90], the first 50 nonvanishing perturbative coefficients are evaluated 
for the ground state, and for the state with parabolic quantum numbers tii = 3, 712 = 0, jti = 1, we 
evaluate the first 70 nonvanishing perturbative coefficients. The apparent convergence of the transforms 
defined in Eqs. (8.11) and (8.26) in higher order is examined. In the case of the Borel-Pade transforms 
defined in Eq. (8.11), use is made of the replacements in Eqs. (8.15) and (8.16) ["leading renormalon poles 
arc being put in by hand"]. This procedure leads to the numerical results listed in Tables 8.1 and 8.2. 
The numerical error of our results is estimated on the basis of the highest and lowest value of the four 
highest-order transforms. 

An important result of the comparison of the methods introduced in Sees. 8.2.2 and 8.2.3 is the following: 
Both methods appear to accomplish a rcsummation of the perturbation series to the physically correct 
result. Method I (Borel+Pade with leading renormalon poles, see Sec. 8.2.2) and method II (Borcl+Pade- 
improved conformal mapping, see Sec. 8.2.3) appear to lead to results of comparable accuracy. 

To date, a rigorous theory of the performance of the resummation methods for divergent series of the type 
discussed in this work (with alternating and nonalternating components) does not exist. The logarithmic 
singularities introduced by the branch points of higher-order polylogarithms [see the index k in Eq. (8.19)] 
are difficult to approximate with the rational functions employed in the construction of Fade approxi- 
mants. A solution to the problem of approximating the logarithmic singularities, based on finite number 
of perturbative coefficients, would probably lead to further optimizimation of the rate of convergence of 
the transformed series. Within the current scheme of evaluation, the problematic logarithmic singularities 
may be responsible, at least in part, for certain numerical instabilities at higher transformation order, 
e.g. in the resuh for T"Ejn{F = 2.1393 x lO"'*) in Eq. (8.32) below. 

For the atomic state with quantum numbers ni = 3, 712 = and m = 1, the evaluation of the transforms 
TEm{F) defined in Eq. (8.11) (method I) and of the transforms T"Em{F) defined in Eq. (8.26) (method 
II) in transformation order M = 67,68,69,70 for a field strength of = 2.1393 x 10"**. Method I leads 
to the following results, 

0.015 860 468 199 2 - i 0.529 048 x 10"^ , 
0.015 860 468 200 9 - i 0.529 047 x 10"^ , 
0.015 860 468 198 9 - i 0.529 048 x 10"^ , 
0.015 860 468 194 5 - i 0.529 015 x 10"^ . (8.32) 



-0.015 860 468 200 4 - i 0.529 047 x 10"^ , 
-0.015 860 468 200 3 - i 0.529 047 x 10"^ , 
-0.015 860 468 200 4 - i 0.529 047 x 10"^ , 
-0.015 860 468 203 3 - i 0.529 046 x 10"*^ . (8.33) 



TEe^{F = 2.1393 x 10""*) = 

TEdfiiF = 2.1393 x lO"'') = 

T£:69(F = 2.1393 x lO"'') = 

TEjo{F = 2.1393 x lO"**) = 
Method II yields the following data, 

T"Eqj{F = 2.1393 X 10"^) = 

T"EesiF = 2.1393 x 10"^) = 

T"EeQ{F = 2.1393 x 10"**) = 

r"En{F = 2.1393 X 10"**) = 
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Numerical results obtained by resummation are presented in Tables 8.1 and 8.2 for a variety of field 
strengths and for the two atomic states under investigation here. 



t (a.u.) 




i 000 ) 


0.04 


-0.503 771 591 013 654 2(5) 


3.892 699 990(1) x 10 


0.06 


-0.509 203 451 088(2) 


5.150 775 0(5) x 10 ^ 


0.08 


-0.517 560 50(5) 


4.539 63(5) x 10~^ 


0.10 


-0.527 419 3(5) 


1.453 8(5) X 10 


0.12 


-0.537 334(5) 


2.992 7(5) X 10 


0.16 


-0.555 24(5) 


7.131(5) X 10 


0.20 


—0.570 3(5) 


1.212(5) X 10 


0.24 


—0.582 6(1) 


1.767(5) X 10 


0.28 


-0.591 7(5) 


2.32(3) X 10-1 


0.32 


-0.600(5) 


2.92(3) X 10-1 


0.36 


-0.604(5) 


3.46(3) X 10-1 


0.40 


-0.608(5) 


4.00(5) X 10-1 



Table 8.1: Real and imaginary part of the energy pseudoeigenvalue Eooo{F) for the ground 
state of atomic hydrogen (parabolic quantum numbers ni = 0, 712 = 0, m = 0). 



F(a.u.) Rc£;30i(^^) r3oi(F) 

1.5560 X 10-* -0.016 855 237 140 761 7(5) 0.421 683(5) x 10-^ 

1.9448 X 10-4 -0.016 179 388 257 0(5) 0.143 773(5) x 10"^ 

2.1393 X lO"'' -0.015 860 468 20(1) 0.105 09(5) x 10-*^ 

2.5282 X lO-'' -0.015 269 293(1) 0.176 39(5) x 10-* 

2.9172 X 10-4 -0.014 742 60(3) 0.999 96(9) x 10-'' 

3.3061 X 10-4 -0.014 260 2(3) 0.295 4(2) x lO'^ 

Table 8.2: Real part and imaginary part of the energy pseudoeigenvalue _E3oi(i^) for the excited state 
with parabolic quantum numbers ni = 3, 712 = 0, m = 1. The field strength F is given in atomic units. 
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925 
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QSfi 
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fi9n 989 '=i'=iQ '=iQ'=i 
u.uzu zoz ooy oyo 


68 


6.206 


411 


X 


10175 


0.620 


004 


547 


622 


620 282 559 595 


69 


-6.752 


937 


X 


10178 


0.620 


004 


547 


622 


0.620 282 559 595 


70 


7.455 


631 


X 


10181 


0.619 


986 


566 


577 


0.620 282 559 595 



exact 0.620 282 559 595 0.620 282 559 595 0.620 282 559 595 



Table 8.3: Evaluation of the perturbation series for Z($) [zero-dimensional 0^-tlieory, see 
Eq. (8.34)] for ,9=1. The s„ are the partial sums of the divergent asymptotic series (8.35), 
the Fade approximants [v — \v /2^/\v /2Y\ are calculated according to Eq. (7.13), and the 
delta transforms are evaluated according to Eq. (7.31). 



8.3 Further Applications of Resummation Methods 
8.3.1 Zero— Dimensional Theories with Degenerate Minima 

We consider the generating functional in a zero-dimensional theory (in this case, the usual path integral 
reduces to an ordinary integral). First, we briefly consider the <i>'*-theory in zero dimensions [see Eq. (9- 
177) ff. in [34]]; the generating functional reads 



d$ 
V2^ 



exp 



2 ^ 



5.34) 



The strictly alternating divergent asymptotic expansion in powers of g for g ^ Q reads, 

°° 4^F(27V-M/2) 

JV=0 



+ 1) 



(-.9) 



N 



.35) 



On using the known asymptotics valid for N oo, which in this case yield the "large-order" asymptotics 
of the perturbative coefficients. 



F(27V+ 1/2) 
F(iV + l) 



4 AT 



r(«) 



?.36) 
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it is easy to explicitly establish the factorial divergence of the series (see also p. 888 of [-')7]). The gen- 
erating functional in zero dimensions has been proposed as a paradigmatic example for the divergence of 
perturbation theory in higher order. It can be resummed easily to the nonperturbative result; in particular 
it is manifestly Borel summable, and no singularities are present on the positive real axis. Specifically, 
in Table 8.2.4, we consider the resummation of the series (8.35) by Pade approximants (Sec. 7.1.4) and 
delta transformations (Sec. 7.1.5). The rapid rate of convergence due to the delta transformations in 
comparison to the Pade approximants is obvious; this finding is consistent with the results of Ref. [24]. 

Complications arc introduced by degenerate minima. As a second example, we consider the modified 
generating hmctional [compare with Eq. (2.6) on p. 15 of [2^.0] and with Eq. (40.1) on p. 854 of [297]]: 



Z'($) 



d$ 



exp 



-Jl>'(l-x/5*) 



d$ 



exp 



.37) 



The expansion of the exponential in powers of the coupling g leads to a divergent asymptotic series, 

N 



N=o _"L ^-^^ ^ 



- - 5 $ 



N=0 



d$ 



N 



2tt 



-1/2*^ ^ 



3=0 



r{2N-j + i) 



4\ J 



N=0 



r(7V- 1/2) 



^ 8^r(2Af + 1/2) ^ 



N=Q 



(8.38) 



where C^(a;) denotes a Gegenbauer (ultraspherical) polynomial. Note that terms of half- integer power 
of g entail an odd power of the field and vanish after integration. The first few terms of the asymptotic 
expansion read, 

Z'($) = 1 4- 6 g + 210 g2 + 13860 ff^ 



-^1351350 / + 174594420 g^ 



4-28109701620 / + 5421156741000 5^ 



-hl218404977539750 5® + 0{g^) . 



5.39) 
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For the perturbative coefficients 

8^r(2Ar + l/2) 

V^r(iv + i) ' (^■^'^^ 

we estabhsh the following asymptotics, 

Cn ^7V-^32^r(Af + 1). (8.41) 

Try 2 

Due to the nonalternating character of the expansion (8.38), it is not Borel summable in the ordinary 
sense. Rather, it is Borel summable in the distributional sense [158, 159]. Here, we present numerical 
evidence supporting the summability of the divergent expansion (8.39) based on a finite number of 
perturbative coefficients. The final integration is carried out along the contour Co introduced in [25] [see 
also Eq. (8.44) below]. The same contour has also been used for the resummation of divergent perturbation 
series describing renormalization group (anomalous dimension) 7 functions [27]. As explained in [25], the 
integration along Co, which is based on the mean value of the results obtained above and below the real 
axis, leads to a real final result if all perturbative coefficients are real. 



M 


partial sum 


TZm 


(.9 = 0.01) 


2 


1.081 000 


1 


.102 


326 


3 


1.094 860 


1 


.096 


141 


4 


1.108 373 


1 


.089 


875 


5 


1.125 832 


1 


.090 


695 


6 


1.153 942 


1 


.092 


000 


7 


1.208 154 


1 


.091 


596 


8 


1.329 994 


1 


.091 


389 


9 


1.642 718 


1 


.091 


553 


10 


2.545 239 


1 


.091 


545 


11 


5.438 230 


1 


.091 


503 


12 


1.5 X 10^ 


1 


.091 


525 


13 


5.5 X 10^ 


1 


.091 


527 


14 


2.2 X 10^ 


1 


.091 


519 


15 


9.5 X 10^ 


1 


.091 


523 


16 


4.5 X 10^ 


1 


.091 


523 


17 


2.2 X 10"* 


1 


.091 


521 


18 


1.2 X 10^ 


1 


.091 


522 


19 


6.9 X 10^ 


1 


.091 


522 


20 


4.1 X 10^ 


1 


.091 


522 


exact 


1.091 522 


1.091 522 



Table 8.4: Resummation of the asymptotic series for the generating functional of a zero- 
dimensional theory with degenerate minima given in Eqs. (8.38) and (8.39). We have g = 
0.01. Results in the third column are obtained by the method indicated in Eq. (8.44) along 
the integration contour Cq (see [25]). The partial sums in the second column are obtained 
from the asymptotic series (8.38). 



In particular, the resummation of the divergent expansion (8.39) is accomplished as follows. We first 
define the Borel transform of the generating functional by [see Eq. (4) in [20] and the discussion after 
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Eq. (8.9)] 



E 



N=0 

Padc approximants to this Borcl transform arc evaluated 



[Af/2l/l(Af + 1)/21 



(z) , (8.43) 



where [x] denotes the largest positive integer smaller than x. We then evaluate the (modified) Borel 
integral along the integration contour Cq introduced in [2"]; specifically we define the transform TZ^ig) 

TZM{g)^ I dteM-t)V'M{gt). (8.44) 

J Co 

In this case, poles above and below the real axis must be considered, and the final result involves no imag- 
inary part. The particular case of g = 0.01 is considered. Values for the partial sums of the perturbation 
series (8.39) and the transforms defined in Eq. (8.44) are shown in Tab. 8.4. The transforms exhibit ap- 
parent convergence to 6 decimal places in 20th order, whereas the partial sums of the perturbation series 
diverge. Between the second and forth term of the perturbation series, (the forth term constitutes the 
minimal term), the partial sums provide approximations to the exact result. It might seem surprising that 
the minimal term in the perturbative expansion is reached already in forth order, although the coupling 
assumes the small value g = 0.01. This behavior immediately follows from the large geometric factor in 
Eq. (8.41) which leads to a "resultative coupling strength parameter" of g^cs = 0.32. "Nonperturbative 
effects" of the order of exp(— l/^ios) provide a fundamental limit to the accuracy obtainable by optimal 
truncation of the perturbation series; this is consistent with the numerical data in Table 8.4. 

We have also investigated the resummation of the divergent series (8.39) via a combination of a conformal 
mapping and Fade approximants in the conformal variable. The situation is analogous to the Stark effect: 
Results are consistent than those presented in Table 8.4 obtained by the "pure" Borel-Pade and in this 
case slightly more accurate. The radius of convergence of the Borel transform Z^{z) defined in Eq. (8.42) 
is s = 1/32 [cf. Eq. (8.17) for the Stark effect], and the appropriate conformal mapping in this case reads 



(8.45) 



(1 + y) 



[cf. Eq. (8.22)]. The inverse reads 



y{w) 



1 - VT~^ 



[cf. Eq. (8.28)]. The conformal mapping (8.45) maps the complex w-planc with a cut along (l,oo) unto 
the unit circle in the complex y-plane. While the zero-dimensional model example given in Eq. (8.37) does 
not exhibit all problematic features of degenerate anharmonic double- well oscillators, certain analogies 
can be established; these comprise in particular the need to evaluate the mean value of Borel transforms 
above and below the real axis (see also [27]). 



8.3.2 The Effective Action as a Divergent Series 

Maxwell's equations receive corrections from virtual excitations of the charged quantum fields (notably 
electrons and positrons). This leads to interesting effects [303, 304]: hght-by-hght scattering, photon 
splitting, modification of the speed of light in the presence of strong electromagnetic fields, and - last, 
but not least - pair production. 

When the heavy degrees of freedom are integrated out (in this case, the "heavy particles" are the electrons 
and positrons), an effective theory results. The corrections can be described by an effective interaction, the 
so-called quantum electrodynamic (QED) effective Lagrangian. The dominant effect for electromagnetic 
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fields that vary slowly with respect to the Compton wavelength (frequencies uj <C 2mc^/h) is described 
by the one-loop quantum electrodynamic effective (so-called "Heisenberg-Euler" ) Lagrangian which is 
known to all orders in the electromagnetic field [304, 305, 306, 307, 308, 309]. 

The Heisenberg-Euler Lagrangian A£, which constitutes a quantum correction to the Maxwell La- 
grangian, is usually expressed as a one-dimensional proper-time integral [see e.g. Eq. (3.43) in [304], 



AC 



lim 

c.?)^0+ 



oo+i„ 



17-1 



ab coth(eas) cot(e5s) — 



-62 



1 



.46) 



To clarify the notation, we introduce the well-known Lorcntz invariants J- and Q which arc given by 



V2 {a' 



.47) 



G = V4^^p.(*-F)^" = E B = ±ab, (8.48) 

where E and B are the electric and magnetic field strengths, Ff^^, is the field strength tensor, and 
(*F)'"' denotes the dual field strength tensor {*F)'^'^ = (1/2) e'^'^P"' Fp^. By a and h we denote the secular 
invariants^ 

b = VV^^+P^~F. (8.49) 

These Lorentz invariants are referred to as secular invariants because they emerge naturally as eigenvalues 
of the field strength tensor; these eigenvalues are conserved under proper Lorentz transformations of the 
field strength tensor. There are connections between the different representations [304]: If the relativis- 
tic invariant Q is positive, then it is possible to transform to a Lorentz frame in which E and B are 
antiparallel. In the case G < 0, it is possible to choose a Lorentz frame in which E and B are parallel. 
Irrespective of the sign of G we have in the specified frame 

a = \B\ and b — \E\ if and only if B is (anti-)parallel to £J . (8.50) 

In any case, because a and b are positive definite, we have 

ab ~ \E ■ B\ > for any Lorentz frame and G , (8.51) 

which clarifies the sign ambiguity in (8.48). We give in (8.47) and (8.48) seemingly redundant definitions, 
but it will soon become apparent that each of the alternative "points of view" has its applications. The 
Maxwell Lagrangian is given by 

Cci = -T= - V4 F^, F^- = V2 - B^) = V2 {b' - a^) . (8.52) 

As it is obvious from Eq. (8.46), the correction AC to the Maxwell Lagrangian is conveniently written in 
terms of the secular invariants a and b. 

It was observed as early as 1956 [310] that the quantum correction (8.46), when expressed as a per- 
turbation series in the usual QED perturbation theory parameter a = e'^/{ATr), constitutes a divergent 
series. This divergent character of the QED perturbative expansion is supported by Lipatov's argument 
based on a saddle-point expansion of the generating functional [311, 312, 313], and explicit estimates for 
the large-order (factorially divergent) behaviour of the QED perturbative coefficients have been obtained 
in [314, 315]. Recently, the divergent character of the large-order behaviour of perturbative expansions 
in quantum field theory has found an exquisite confirmation in explicit 30-loop calculations of renormal- 
ization group 7 functions in a six-dimensional (j)^ theory, and in a Yukawa theory [316, 317]. 

Here, we will study the divergent series generated by expanding the effective action (8.46) in powers of e^. 
We will distinguish the cases of a magnetic and an electric background fields. We start with a magnetic 
field of strength B in which case the effective action reads 

SB . ACiE . „, B, . J - i - 1} oxp (-^ ,) (8,53) 
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[put a = B, b ~ in Eq. (8.46)]. The nonpcrturbative results in Eq. (8.53) can be expanded in powers of 
the effective coupUng 

9B = . (8.54) 

This results in the divergent asymptotic series 



2e"i?^ ^ (-l)"+i4" |62„+4| 



7r2 ^ (2?i + 4) (271 + 3) (2n + 2) 



3b 



262^2 
2~9B 



720 1260 ^ 630 



1 2 

9B- -^gB 



+ ...-2.33X 10^"'\g|" + . . .] (8.55) 

for ge —>■ 0. Here, |;S2n+4| denotes the modulus of the (2 n + 4)th Bernoulli number. The Bernoulli numbers 
alternate in sign. 

sign(62„+4)-(-l)"+^ (8.56) 

The expansion coefficients 

(-1)»+1 4" \B2n+4\ 

(2n + 4)(27i + 3)(2n + 2) ^ ' 

obviously display an alternating sign pattern and grow factorially in absolute magnitude, 

»-ta:ii£eiL±a(,,o(2--«.,) (a.») 

for n —^ oo. 

The series (8.55) is an alternating series. In [163, 255, 310], it was shown that the series (8.55) is Borel 
summable. A further expansion of the perturbative coefficients c„ for large n is carried out in [255] in 
order to perform the Borel summation; this analysis follows the asymptotic expansion in Eq. (8.58) above, 
and higher-order terms in the asymptotics of the coefficients for large n are also used. Here we consider 
the resummation of the divergent series (8.55) with the help of rational approximants which use as input 
data only a finite number of perturbation theory coefficients in numerical form. Specifically, we use Pade 
approximants defined in Eq. (7.12) and discussed in Sec. 7.1.4 and the nonlinear sequence transformations 
d and S defined in Eqs. (7.30) and (7.31), respectively. The nonlinear transformations were discussed in 
Sec. 7.1.5. 

The transformations are applied to the divergent perturbation series Eq. (8.55), the transforms are calcu- 
lated with the help of the recurrence relations (7.27) and (7.28). The evaluations presented in Tables 8.5 
and 8.6 are carried out for values gs = 1/10, = 1 and = 10 of the expansion parameter defined 
in Eq. (8.54). The first column in the tables contains the index n, the second column contains the partial 
sums of the input series and in the last three columns are the results of the Fade Approximation, the 
d-transformation and the delta-transformation. 

The data in the tables is presented to a numerical accuracy of 10^^. The convergence of the transforms 
is indicated by underlining those decimal places which appear to have converged to the nonpcrturbative 
result. The numerical data in the Tables 8.5 and 8.6 indicate that the most favourable numerical results are 
obtained with the Weniger delta-transformation. For the numerical calculations, exact rational arithmetics 
was used [!iii]. Numerical experiments for g^ larger than 1 (not shown) and data presented in [24] indicate 
that the Weniger delta transformation appears to resum the divergent series (8.55) at least up to ^b = 200, 
and to arbitrary precision. The other resummation prescriptions - Fade, and the Levin d-transformation 
- fail in this domain of strong coupling. 

In the case of a background electric field, Eq. (8.46) yields the following expression 

. AA., i» ^ 0, . r $ {coU, . - 1 - i } CP (-i .) . (8.58, 
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— 91 1 393 673 
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10 


^.308 


952 


549 


-0.211 393 


982 


-0.211 393 670 


-0.211 393 671 


exact 


-0.211 393 671 


-0.211 393 671 


-0.211 393 671 


-0.211 393 671 



Table 8.5: Evaluation of the perturbation series for 5*3 given in Eq. (8.53) for ge = 1/10. Results are 
given in terms of the scaled function S^b given by Sb = 10~^ x (— i?^)/(8 tt^) x S^b- The performance of 
three different resummation methods is compared: in the second column, Padc approximants are used. In 
the third column, we list results obtained using the Levin transformation defined in Eq. (7.30), whereas 
in the last column, the delta transformation (7.31) is employed. The apparent convergence is indicated 
by underlining the decimal figures that have stabilized in increasing transformation order. 



In analogy to the expansion parameter for the magnetic case from Eq. (8.54), we can define the effective 
coupling as 

5E = ■ 8.60 

The expansion of the effective action in powers of gE leads to the formal power series 



9E J2 
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7r2 (2n + 4)(2n + 3)(2n + 2) 



9i 



—9E 



1 



1 



720 1260^ 630 



5E 



:9e 



+ ... + 2.33 X 10i"^.g|" + . . .] . (8.61) 

The perturbation series is nonalternating. Because the remainder estimate tOn = On+i [see Eq. (7.29)] 
is in general not valid for nonalternating series, we cannot assume that any of the previously discussed 
resummation methods (Pade, or the delta or d-transformation) are able to resum the perturbation series 
Eq. (8.61). This is confirmed by the explicit numerical data presented in Table 8.3.2. 

There exists a "nonperturbativc" imaginary part for the uniform background electric field which is a 
priori not contained in the perturbation scries whose coefficients are all real and positive [see Eq. (8.61)]. 
The imaginary part can be directly inferred from Eq. (8.59) by residue calculus, the result is [34, 255] 

i-^--^E;;i-ph^ ■ (8-62) 



The characteristic factor 



\ e E 

n—l ^ 



exp ( -— ) (8.63) 
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25 
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exact 


-0.164 598 939 
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Table 8.6: Evaluation of the perturbation series for 5b given in Eq. (8.53) for = 1- Results are given in 
terms of the scaled function 5b given by S'b = 10~^ x (— i?^)/(8 tt^) x 5b. The apparent convergence 
is indicated by underlining the decimal figures that have stabilized in increasing transformation order. 



is nonpcrturbative since Eq. (8.62) has an essential singularity in the limit — > 0, and it is reminiscent 
of characteristic expressions occurring in the description of quantum mechanical tunneling processes (see 
p. 195 in [34]). In our context, the imaginary part (8.62) describes the pair production in an electric 
background field. As is well known, the vacuum in very strong electric fields becomes unstable with 
regard to particle-antiparticle pair production (see also the elucidating discussion in [303]). 

Spontaneous pair production becomes considerable only in strong external electric fields. Ordinary field 
strengths are much smaller than the critical field i^crit for which 

= 1 . (8.64) 

mi 

The critical field strength £^crit can be written (including h and c) as the ratio of the rest mass of the 
electron to the characteristic length scale which is set by the Compton wavelength of the electron: 

i?cHt = ^ = ^ - 10-^ . (8.65) 

eft/TOeC eft cm 

This critical field strengths are comparable to those felt by an electron in a hydrogenlike atom whose 
nucleus has the charge Z w 137. 

Nonalternating divergent perturbation series require considerable effort in their resummation. It has been 
shown in [25] that a generalized Borel-Pade method can be used to reconstruct the imaginary part of 
the full nonpcrturbative result. This method depends on the special integration contours in the complex 
plane, introduced in Sees. 8.2.2 and 8.2.3. In Table 8.8, we present results for the generalized Borel 
transforms constructed according to Eqs. (1) - (6) of [25]). 
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Table 8.7: Failure of the evaluation of the perturbation series for 5'e (electric background field) given 
in Eq. (8.59) for g-^ = 1. Results are given in terms of the scaled function Se given by Se = 10^^ x 
(— E'^)/{8 tt^) x Se- Observe the apparent lack of convergence. 



As for the complex resonances of the Stark effect (see the numerical results in Sec. 8.2.4), it is possible to 
reconstruct the nonpcrturbativc imaginary part (pair production rate) from the purely real pcrturbative 
coefficients in Eq. (8.61). 

8.3.3 The Double-Well Problem 

We discuss in the current section specific results, recently obtained in [2!)], for the quantum-mechanical 
double-well problem and the related multi-instanton expansion. The energy levels of the double-well po- 
tential receive, beyond perturbation theory, contributions which are non-analytic in the coupling strength; 
these arc related to instanton effects. For example, the separation between the energies of odd- and 
even-parity states is given at leading order by the one-instanton contribution, which, in the path integral 
formalism, corresponds to a tunneling of the particle from one minimum of the potential to the other min- 
imum. In order to determine the energies (of individual levels as well as the separations) more accurately, 
multi-instanton configurations have also to be taken into account. The multi-instanton configurations 
describe a particle that tunnels between the minima more than once. For even instantons, the particle 
returns to the minimum from which it started the motion. 

For the double-well problem, the energy eigenvalues of the states at nonvanishing coupling .9 7^ cannot 
in principle be obtained by analytic continuation from the unperturbed situation at vanishing coupling 
g ~ because a potential with degenerate minima introduces a degeneracy in the spectrum: for any 
one unperturbed state, two states emerge when the perturbation is switched on. These two states are 
separated by an energy shift which is nonpcrturbativc and nonanalytic in the coupling, i.e. vanishing 
to any order in perturbation theory. Therefore, the two states are described by the same perturbation 
series and yet differ in their energy by instanton contributions. Specifically, we consider the case of the 
double-well potential with the hamiltonian 

ff = -f^ + ^n'7), Viq) = lj'il-q)\ (8.66) 

Of course, the alternative formulation H = —d^/dq^ + (1 ^ y/9l)'^ gives rise to the same energy 
levels, but the form (8.66) illustrates that the coupling plays the formal role of h. The hamiltonian 
(8.66) can be obtained from H = —d^/dq"^ + ^q'^ {1 — \/gq)^ by the scaling q q/^- It has been 
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exact 0.001 147 154 0.001 147 154 + 12.52 x 10-' 



Table 8.8: Resummation of the asymptotic series for the QED effective action (8.61) in a constant back- 
ground electric field for = 0.05. Results in the third column are obtained by the method discussed 
in Sees. 8.2.2 and 8.2.3 along the integration contour C+i (see Fig. 8.1). The partial sums in the second 
column are obtained from the asymptotic series (8.61). 



conjectured [171, 172, 173, 174] that an asymptotic expansion for the energy eigenvalue can be obtained 
by finding a solution to the equation 

1 \ / 2\^(^'^) 

D{E,g)][—] exph^(£;,.g)/2] =±i, (8.67) 



271 \2 J \ 9 

which can be understood as a modified Bohr-Sommerfeld quantization condition. The plus and minus signs 
apply to even- and odd-parity states, respectively The conjecture (8.67), whose validity has been proven 
in [17G], has found a natural explanation in the framework of Ecalle's theory of resurgent functions [175, 
177, 318]. The functions D{E,g), A{E,g) constitute power series in both variables. The function D{E,g) 
describes the perturbative expansion; its evaluation is discussed in [171, 174]. The first terms read 

D{E, g)^E + g (^e'^ + ^+ g^ (ih E^ + ^E^+ 0{g^) . (8.68) 

The ground and the first excited state are both described by the same perturbation series which can 
be found by inverting the equation D(E,g) = 1/2. When the energy is expressed in terms of the naive 
perturbation series in g, the function [D{E,g) — N — 1/2] then vanishes in any order of perturbation 
theory, i.e. in all orders in g. Here, N is the quantum number of the unperturbed state which is a 
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harmonic oscillator eigenstate. The function A{E,g) essentially describes instanton contributions [174]; 
its first terms read 

A{E, .9) = ^ + 5 (l7 + + .9' (227 E^ + ^E^+ 0{g^) . (8.69) 

A solution to the equation (8.67) can be found by systematically expanding the energy eigenvalue E{g) 
in powers of g and in the two quantities 

m-lnU) and ^{g)^'-^^^m. (8.70) 
V 9 J 



Terms of order ^(g)" belong to the n-instanton contribution. The energy eigenvalue for nonvanishing 
perturbation 5 7^ can be described by two quantum numbers: the unperturbed quantum number N and 
the positive or negative parity of the state. We have (the upper index denotes the instanton order) 

i?^,±(5)-E<±(5) (8-71) 

n=0 

where the perturbation series (zero-instanton contribution) is given as 

<±(5)=E<^5^' (8.72) 



where the right-hand side is parity independent. For 71 > 0, the instanton contribution reads 

Nil n-1 



= (-] agr E ^(5)' E ^^i^' 9' . (8.73) 



fc=0 (=0 



The lower indices n, k and I of the e coefficients denote the instanton order, the power of the logarithm 
and the power of g, respectively. Some of the results that will be used in the sequel read, 
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,(0,+) _ (0,-) _ (0,+) _ (0.-) _ _^ _ ^ ^ 

"^200 ~ *'200 ~ ' ' '^201 ~ '^201 ~ 9 6 



(0'+) ("'-) 1277 

£202 = 4o2 ^ ' (8-^4) 



where 7 = 0.57221 ... is Euler's constant. Odd-instanton contributions have opposite sign for opposite- 
parity states and are responsible, in particular, for the energy difference of the ground state with quantum 
numbers (0, -I-) and the first excited state with quantum numbers (0,—). The dominant contribution to 
the separation of the two lowest energy levels is given by the one-instanton contribution: 

/ 71 6999 \ 
i?o,-(.g)-£;o,+ (.9)^2e(5) (^1 - - g - ^ 5' + ^(s') j + ^(^(5)') • (8.75) 

By contrast, even-instanton contributions have like sign for opposite-parity states and are responsible, in 
particular, for the displacement of the mean value (1/2) [i?o._ (g) -|-£'o,+ (3)] from the value of the general- 
ized Borel sum of the perturbation series B (^X]?=o ff^) (^°^ evaluation of the generalized Borel 
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sum of a nonaltcrnating divergent series, see Sec. 8.3.1). The dominant contribution to the displacement 
comes from the two-instanton effect, and we have 



\K=0 




53^^/2eT\ 23 



1277^^ /2e-'\ 13 



-0(5'ln(5)) 



+0(?(5)') 
The function [172] 



A(.9) = 4 



[Eo.-ig) - EoAg)]' 



■"(¥) 



.76) 



.77) 



relates the multi-instanton contributions to the energy eigenvalues, which can be evaluated numerically, 
and to the (generalized) Borel sum of the perturbation series which is evaluated by analytic continuation 
of the integration path into the complex plane (see [297]). The calculation of A{g) at small coupling is 
problematic because of severe numerical cancellations. From the equations (8.75), (8.76) and (8.77), we 
obtain the following asymptotics for A{g), 
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If we additionally perform an expansion in inverse powers of ln(2/g) and keep only the first few terms in 
{l/ln(2/g)} in each term in the g-expansion, the result reads 
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The higher-order corrections, which arc only logarithmically suppressed with respect to the leading terms 
1 + 3^, change the numerical values quite significantly, even at small coupling. In Table 8.9 we present 
numerical results for the function lS.{g) at small coupling; these are in agreement with the first few 
asymptotic terms listed in equation (8.78) up to numerical accuracy. Of course, for strong coupling, 
significant deviations from the leading asymptotics must be expected due to higher-order effects; these 
are indeed observed. For example, at g = 0.1 the numerically determined value reads A(O.l) = 0.87684(1) 
whereas the first asymptotic terms given in equation (8.78) sum up to a numerical value of 0.86029. 

The higher-order corrections to the two-instanton effect are related to the corrections to the leading 
factorial growth of the perturbative coefficients. This can be seen by expressing that the imaginary part 
of the perturbation series, when continued analytically from negative to positive coupling, has to cancel 
with the imaginary part of the two-instanton contribution which is generated by the logarithms ln(— 2/^). 
The corrections of order g h\{—2/g) and g^ h\{—2/g) yield the \/K~ and 1/i^^-corrections to the leading 
factorial growth of the perturbative coefficients. From the results for £21',^'' ~ 0, 1, 2) given in equation 
(8.74), we obtain 
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Tabic 8.9: Comparison of numerical values for the function A(g) defined in equation (8.77) in the 
region of small coupling to values obtained by calculating the first few terms in its asymptotic 
expansion given in (8.78). 



coupling g 


0.005 


0.006 


0.007 


0.008 


0.009 


A((7) num. 
A(g) asymp. 


1.0063(5) 
1.00640 


1.0075(5) 
1.00739 


1.00832(5) 
1.00832 


1.00919(5) 
1.00919 


1.00998(5) 
1.01001 



The analytic results should be checked against explicit values of the perturbative coefficients. We have 
determined the first 200 perturbative coefficients e'^^^ {K = 0, . . . , 200) of the perturbation in the form 
of rational numbers, i.e. to formally infinite numerical accuracy. This allows to verify the l/K~ and 
l/if^-corrections to the leading factorial growth in equation (8.80) to high accuracy, for example by 
employing Richardson extrapolation [319]. Using the 160th through the 200th perturbation coefficient as 
input data for the Richardson algorithm, the coefficients of the leading, of the 1/if -subleading and of the 
l/K"^ suppressed corrections are found to be consistent with the analytic results given in equation (8.80) 
up to a relative numerical accuracy of 10~^^, 10"^'^ and 10~^°, respectively. For completeness, we give 
here the numerical values of the 198th through the 200th perturbative coefficients, to 30 decimals. These 
read: 

E^°lgs = -5.50117 76962 885879352775694 38632 X lO^^'^ , 
E^°lgg = -3.28445 39841 65780 00616 21912 32835 X 10''^^ 

^o°2oo = -1-97082 14193 09543 76979 53006 07410 x 10''^° . (8.81) 
Values for all 200 coefficients are available [234]. 

It is an interesting consequence of the expansion (8.71) that the energy difference (£^0,- — £'0,+), at small 
coupling, is described to high accuracy by the one-instanton contribution (n = 1 in equation (8.73)). For 
g = 0.001, we obtain to 180 decimals, 

£;o,+ (0.001) = 0. 49899 54548 6210917168 91308 3948192163 68209 47240 
20809 66532 93278 69722 01391 15135 28505 3829445798 
45759 95999 06739 55175 84722 67802 81306 96906 01325 
25943 77289 94365 88255 24440 17437 12789 27978 99793 , (8.82) 

0. 49899 54548 62109 17168 91308 39481 92163 68209 47240 
20809 66532 93278 69722 01391 29839 92959 55803 70812 
27749 92448 48259 36743 64757 68328 84835 35511 34663 

06309 82331 51885 23308 08622 84780 52722 10103 67282 . (8.83) 

Decimals which differ in the two energy levels are underlined. The results have been obtained by lattice 
extrapolation using a modified Richardson algorithm which is constructed according to ideas outlined 



whereas 
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in [156]. Calculations were performed on IBM RISC/6000 workstations while making extensive use of 
multiprecision libraries [235, 237, 238]. We define Vnig) as the Afth partial sum of the onc-instanton 
contribution £'g"^l(p) — s'^^^g), 



M 

VM{9) = 2^{g)Y,e^°,-^g^ . (8.84) 

Using exact rational expressions for the coefficients ef^^ (j < 141), we obtain 
T'moIO.OOI) X 10^1 = 

1.47046 44541 75092 50138 19899 64494 15198 15678 00350 05260 35283 

86053 33378 03660 50415 75193 50528 41826 73433 99328 21246 74888 , (8.85) 



75141(0.001) X 10^1 = 

1.47046 44541 75092 50138 19899 64494 15198 15678 00350 05260 35283 

86053 33378 03660 50415 75193 50528 41826 73433 99328 21246 74887 . (8.86) 

These values are in excellent agreement with the numerically determined energy difference (see the results 
presented above in equations (8.82) and (8.83)) 

[£;o.^(o.ooi) - £:o^+(o.ooi)] x lo^^ = 

1.47046 44541 75092 50138 19899 64494 15198 15678 00350 05260 35283 

86053 33378 03660 50415 75193 50528 41826 73433 99328 21246 74887 . (8.87) 

The first 70 decimals in equations (8.82) and (8.83) are the same because the one-instanton contribution 
is of the order of 1.4 x 10^^^. The accuracy to which the one-instanton contribution describes the energy 
difference £^0,- (0.001) — i?o,+ (0.001) is limited by the threc-instanton effect which for g = 0.001 is of 
the order of 8 x 10~^^^. Note that the two-instanton effect (which for g = 0.001 is of the order of 
4 X lO^^^'^) does not limit the accuracy to which the onc-instanton contribution describes the energy 
difference because it has the same sign for opposite-parity states. 

We have demonstrated that the behavior of the characteristic function A{g) defined in equation (8.77) 
at small coupling is consistent with higher-order corrections to the one- and two-instanton contributions, 
specifically with the instanton expansion of the energy levels governed by the equations (8.71) and (8.73), 
with the assumption that the instanton contributions given by equation (8.73) should be Borel summed, 
with the explicit results for the higher-order coefficients listed in (8.74) and the analytically derived 
asymptotics for the function A{g) given in equation (8.78). The corrections of relative order 1/K™ 
to the leading factorial growth of the perturbative coefficients - see equation (8.80) - are consistent 
with the analytically evaluated ln(— 2/5)-corrections to the two-instanton effect and with the explicit 
values for the first 200 terms in the perturbation series (8.72). The nonperturbative energy difference 
Eo^-{g) — i?o,+ (3) at small coupling g is described, to high accuracy, by the one-instanton contribution 
only. 

8.4 Divergent Series: Some Conclusions 

A priori, it may seem rather unattractive to assume that the quantum electrodynamic perturbation series 
may be divergent even after the regularization and the rcnormalization. However, as shown by explicit 
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nontrivial 30-loop calculations of rcnormalization group 7 functions in a six-dimcnsional (j)^ theory, and 
in a Yukawa theory (presented recently in [316, 317]), we believe that the ultimate divergence of the 
perturbative expansion can be regarded as a matter-of-fact, clearly demonstrated by explicit high-order 
calculations. Therefore, it appears meaningful to explore the physical implications of this divergence. 

We have discussed four physical applications of resummation methods for divergent series: (i) the energy 
displacement of a hydrogen atom in a background electric field (Sec. 8.2.4), (ii) zero-dimensional (model) 
field theories (Sec. 8.3.1), (iii) the QED effective Lagrangian (Sec. 8.3.2), and (iv) the energy levels of the 
double-well potential (Sec. 8.3.3). The mathematical structure of these problems can be characterized as 
follows, in the order of increasing complexity: 

• The QED effective Lagrangian for a background magnetic field, expressed as a perturbation series 
in a, is manifestly Borel summable [163, 255, 310], and the nonperturbative result can be inferred 
by a number of different resummation methods, as discussed in Sec. 8.3.2. 

• The QED effective Lagrangian for a background electric field, as well as the perturbation series for 
a hydrogen atom in an electric field, is not Borel summable. The same applies to the perturbation 
series for a zcro-dimcnsional field theory with degenerate minima. However, these three problems 
admit a treatment according to the concept of generalized (distributional) Borel summability [158]. 
Numerical results (see [25, 28] and Sees. 8.2.4 and 8.3.2) confirm the distributional Borel summabil- 
ity, and the use of additional asymptotic information about the perturbative coefficients ["leading 
renormalon poles", see Eqs. (8.15) and (8.16) in Sec. 8.2.2] accelerates the convergence of the Borel- 
transformed perturbation series. 

• The double-well potential represent a very problematic case: the full nonperturbative solution to 
the eigenvalue problem cannot be obtained in principle from perturbation theory. The reason is the 
following: the perturbation introduces an additional degeneracy in the spectrum: each unperturbed 
level splits into two energy levels when the perturbation is nonvanishing. Both of these levels are 
described by one and the same naive perturbation series (see Sec. 8.3.3). The energy difference is 
nonperturbative and nonanalytic in the coupling strength, and finds an explanation in the theory 
of instantons. 

For the Stark effect and the QED effective Lagrangian (electric field), the existence of nonperturbative 
contributions is intimately linked with the failure of the Carleman criterion (see for example [163], The- 
orems Xn.17 and Xn.l8 and the definition on p. 43 in [l(i4], p. 410 in [1i'5], Ch. 6 of this Thesis or the 
elucidating discussion in Ref. [16()]). The Carleman criterion determines, roughly speaking, if nonanalytic 
contributions exist for a given effect which is described by a specified perturbation series. In this sense, the 
divergence of the perturbative expansion is physically important: for the QED effective action, it allows 
for the existence of nonanalytic, nonperturbative contributions like the pair-production amplitude which 
is not contained in the real perturbation series (8.61) and cannot be obtained on the basis of perturbation 
theory alone. 

Also, we would like to illustrate here the utility of resummation methods in those cases where perturbation 
theory breaks down at large coupling. As explained in Sees. 8.2.4, 8.3.1 and 8.3.2, even in situations where 
the perturbation series diverges strongly, it can still be used to obtain meaningful physical results if it is 
combined with a suitable resummation method. In a relatively weak field, it is possible to obtain more 
accurate numerical results by resummation than by optimal truncation of the perturbation series (see 
also [25]). In a strong field, it is possible to obtain physically correct results by resummation even though 
the perturbation series diverges strongly (see the discussion in Sec. 8.2.4 and the data in Tables 8.1, 8.2 
and 8.4). By resummation, the perturbation series which is inherently a weak-coupling expansion can be 
given a physical interpretation even in situations where the coupling is large. Returning to the analogy 
to quantum field theory, one might be tempted to suggest that physically complete results can in many 
cases be obtained after regularization, rcnormalization and resummation. 

Can the full nonperturbative result (for energy levels etc.) be inferred in all cases by a resummation of 
the divergent perturbation series? The answer is, unfortunately, no. This is demonstrated in Sec. 8.3.3 
by way of example. In theories with degenerate minima (such as the double- well problem), corrections 
to energy levels are caused by so-called instanton contributions that follow naturally by an expansion 
of the path integral around nontrivial saddle points. The n-instanton contribution is characterized by a 
nonperturbative factor ^(g)" where [see Eq. (8.70)] ^(p) = exp[— l/(6g)]/y^, and g is the couphng. The 
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expansion of ^(g) in powers of g vanishes in all orders of g. The nonperturbative energy splitting between 
even- and odd-parity states is investigated by considering the function A{g) defined in Eq. (8.77). This 
function would vanish if the generalized Borel sum of the perturbation series could reproduce energy 
levels exactly. However, it does not vanish, as demonstrated by the data in Tab. 8.9 (the numerical data 
can be interpreted naturally by considering so-called two-instanton eff'ects). At the same time, in a zero- 
dimensional theory with degenerate minima (considered in Sec. 8.3.1), resummation is successful despite 
the existence of a nontrivial saddle point. The reason for the success of the resummation in this model 
problem appears to lie in the fact that no degcneraces are introduced by the additional saddle point. The 
perturbation series determines the generating functional uniquely, which is not the case in the double- well 
potential where two distinct energy levels (of opposite parity) share the same perturbation series. This 
intriguing situation will warrant further investigation in the near future. 
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Chapter 9 

Conclusions 



Wc proceed to the interpretation of the results obtained in this Thesis, which has a dual subject: "quantum 
electrodynamic bound-state calculations and large-order perturbation theory" . 

The significance of quantum electrodynamics as one of the most appealing physical theories (a "jewel" of 
theoretical physics according to Richard Feynman) does not require any further explanation. The esthetic 
appeal of bound-state quantum electrodynamics stems from the accuracy of the experimental verifications, 
the significance of the theory for the determination of the fundamental constants, and the conceptual 
complexity of the calculations which derives from the apparent simplicity of the physical systems under 
study, when an accurate understanding is required in higher orders of perturbation theory. However, 
the divergence of the perturbative expansion in large orders, which persists even after the regularization 
and renormalization, raises fundamental questions regarding the internal consistency of the theoretical 
predictions. The problems associated with the divergence of the perturbative expansion do not only 
plague quantum electrodynamics, but even occur in ordinary quantum mechanical perturbation theory. 
This Thesis is divided into two Parts, which represent two attempts to advance our understanding of the 
higher-order corrections as well as of the conceptual questions raised by the divergence of the perturbative 
expansion. 

In the first Part, we discussed aspects of certain quantum electrodynamic bound-state calculations, which 
have been carried out over the past years. Not all of these calculations [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 
14, 15, 16, 17, 18] could be described in full detail. Both the experimental accuracy as well as the accuracy 
of the theoretical calculations have improved dramatically in recent years [42]. On the theoretical side, the 
calculations at low nuclear charge number may be based on either of two methods: the analytic approach 
that is based on the Za-expansion and the numerical approach in which all electron propagators are 
kept in exact relativistic form. Progress has been achieved due to advances in the analytic approach that 
profit from an adequate formulation of the problem, and due to the development of numerical techniques 
which have led to the highly accurate numerical evaluation of QED corrections in the realm of low nuclear 
charge number. This has led to the favourable situation where analytic and numerical calculations can 
be checked against each other. In view of the complexity of the calculations and their importance for 
the determination of the fundamental constants [s!)], the existence of independent cross-checks for the 
calculations is highly desirable. 

Numerical calculations at low Z have been made possible by convergence acceleration methods described 
in the second Part of this Thesis (notably by the combined nonlinear-condensation transformation) . These 
techniques have led to the solution of the severe numerical difficulties associated with the singularity of 
the propagators for equal radial arguments, and to a reduction in computing time by three orders of 
magnitude. Results obtained for the one-photon self energy are several orders of magnitude more accurate 
than previous calculations of the effect. They represent the first direct evaluations of the one-loop self 
energy at the nuclear charge numbers Z — \ and Z = 2, which are of crucial importance for precision 
spectroscopy (atomic hydrogen and hydrogenlike helium). 

The analytic calculations of the one-loop self energy for higher excited states will be complemented in the 
near future by numerical calculations based on the new techniques [DH], leading to a cross-check of the 
type mentioned above. Severe difficulties associated with the multitude of analytic terms due to the more 
complicated wavefunctions of the highly excited states and the associated angular momentum algebra 
have been a problem for self energy calculations in this area. 
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The calculations for the bound-state two-loop self energy report on the first evaluation of the highly 
problematic nonlogarithmic higher-order two-loop binding correction Bqq. Let us recall that analytic 
work on the corresponding Aqq correction for the less involved one-loop problem has extended over three 
decades [32, 68, 69, 70, 71]. The development of analytic methods for the evaluation of Bgo leads to 
an improved understanding of the scaling of the two-loop effect at low Z and will enable a detailed 
comparison of analytic and possibly available numerical results for the two-loop effect in the future. At 
the same time, these results lead to improved predictions for the hydrogen and helium fine structure. 

The spin-dependence of quantum electrodynamic corrections is a conceptually important issue. We focus 
on spinless particles and present a simplified derivation of the generalized Breit hamiltonian for a system 
of two spinless particles. Although scalar QED is a renormalizable theory, we found the issue of bound- 
state calculations within this theory insufficiently addressed in the literature. Our calculations for the 
leading-order self energy and the relativistic recoil correction are the first corrections of the "self energy 
type" to be evaluated within bound-state QED involving only scalar particles. 

As a further application, we would like to mention the recent investigation [18]. 

In the second Part of the Thesis, we investigate convergence acceleration and resummation methods 
for divergent series. The convergence acceleration techniques address the fundamental problem of the 
acceleration of a slowly convergent, nonalternating series. Severe numerical instabilities associated with 
the formation of higher-order weighted differences for nonalternating series have been a major obstacle 
for algorithms that try to accomplish this task. The methods find applications in areas as diverse as DNA 
sequence analysis and experimental mathematics. 

Recent highly nontrivial 30-loop calculations [316, 317] have convincingly demonstrated the divergence 
of the perturbative expansions in quantum field theory originally conjectured by Dyson [320]. We discuss 
nonlinear sequence transformations as an alternative to the "usual" Borel method for the resummation of 
divergent series that result from perturbative expansions in quantum field theory. These transformations 
have favourable asymptotic properties and lead in many cases to better numerical results than Pade 
approximants. We also discuss generalizations of the Borel method which lead to consistent results even 
if the full, nonpcrturbativc physical energy level or vacuum-to-vacuum amplitude acquires an imaginary 
part due to quantum mechanical tunneling. This imaginary part, which represents the autoionization 
decay width (Stark effect) and the electron-positron pair-production amplitude (in the case of the QED 
effective action), can be derived starting only from the real (not complex!) perturbative coefficients. 
In the double- well problem, the situation is different. Even- and odd-parity states acquire an energy 
separation which is nonperturbative and nonanalytic in the coupling strength and is an effect which in 
principle cannot be derived from perturbation theory alone. The energy separation is given by so-called 
multi-instantons which follow naturally from a path integral representation of the partition function 
and correspond to configurations where the particle tunnels repeatedly between the two minima of the 
double-well potential. The results lead to a better understanding of the energy shifts due to multi- 
instanton effects, by evaluating higher-order corrections to the two-instanton effect, and they also lead to 
an accurate verification of the instanton expansion via a comparison to numerically determined energy 
levels. 

The common theme of all investigations discussed in is to to explore the predictive limits of quantum 
theory. It has been the aim of this Thesis to present results for some of the essential QED corrections 
which influence the spectrum of bound systems and are of current experimental interest, as well as 
to work towards a solution of the questions regarding the predictive limits of field theories set by the 
ultimate divergence of the perturbative expansion. Certain methods developed along this endeavour have 
meanwhile found applications in other areas (applied biophysics and mathematics). 
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